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5.11 03JECTIVE FUNCTION 

If r s represents the cost of steel in rupees per 
unit weight of steel and r^, the cost of concrete in rupees 
per unit volume of concrete, the objective function, being 
the material cost per running metre of wall, can be expressed 
as 


f (X) 


st 


+ r V 
c c 


5.27 


5.12 DATA FOR THE PROBLEM 

5.12.1 Soil Properties 

Five different classes of soil are considered and 

they are; gravel, dense sand, loose sand, silt and clay. The 

first three types belong to the cohesion less variety while 

the other two to the cohesive one. Referring to IS t 1498 

(1970), Terzaghi (1960), Tschebotarioff (1951) and Huntington 

(1967), representative values are selected for the unit 

weight (w ), angle of internal friction (^) , permissible 
s 

beaming capacity (o ) and cohesion (c) for the soil retained * 

■ per 

These values are given in Table 5.1. 

5.12.2 Pro assigned parameters 

Three different values are selected for the top 
width b, namely 20 cm, 25 cm and 30 cm. The range of overall 
height h chosen is from 3.0 to 6.0 m. 
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= vector of generalized nodal displacements 

= modified displacement vector in a standard 
eigenvalue problem 

= subvectors of generalized displacements at 
nodes A, B and c of a member 


{q x > / {q Y > | 
{q Z },{q<P} I 


subvectors of generalized displacements 
= corresponding to flexure -in principal planes, 
axial stressing and flexural torsion respect- 
ively at a node in a member 



Following such a rational approach/ the optimal 
solution obtained is obviously a local minimum. For sane 
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cases, such optimal solutions were successfully checked for 
the Kuhn -Tucker conditions for the local minimum. However, 
afterwards it was felt unnecessary to check every optimal 
solution for these conditions and as such they are not 
incorporated in the programmes 0PT2, 0PT3 and 0PT4. The 
results obtained for the various cases in the programme 0PT2 
are shown in Tables 6.2 through 6.37. The results obtained 
from the programmes 0PT3 and 0PT4 are graphically represented 
m Figures 6.1 through 6.4. 

6.1.2 Programme OPTl 

A few optimal solutions obtained from this programmes 
were also successfully checked for the Kuhn-Tucker conditions 
for the local minimum. Later on, however, these conditions 
have not been included in the programme OPTl , as it was 
not found necessary. The results have been graphically 
represented in Figures 6.5 through 6.7. 

6.2 DISCUSSION OF RESULTS 

The results are discussed in the order in which 
the various optimisation programmes have been handled in 
Chapters 3 through 5. 
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R 

{R} 

2 

r 

r 

r 


S ± ,S 

SR, 


{s i> 


SR, 

4 

T 

T s 
T 


a . 


tf. 


^ lmax ' ^ lmin 


tf , 


2max ' ^ 2min 


tw 


tw _ / tw . 
max mxn 


v 

I V 


V, 


= resultant foundation pressure on the base of a 
retaining wall 

= vector of fixed end reactions for a member 
= cross-sectional parameter 
= penalty parameter 
= unit cost of concrete '■ 

= .unit cost of steel 
= shearing strength of soil 
= direction vector 

= slenderness ratio of compression flange of an 
I -section 

= slenderness ratio of a column 

= total twisting moment applied at a section 

= St'. ( Variant's torque at a section 

= warping-., torque at a section 

= thickness of top flange of an I -section 

= maximum and minimum values of top flange 
thickness of an I -section 

= thickness of bottom flange of an I -section 

= maximum and minimum values of bottom flange 
thickness of an I -section 

= thickness of web of an I -section 

= maximum and minimum values of web thickness 
of an I -section 

= potential energy of a system 

= summation of vertical forces on the base of a 
cantilever retaining wall 

= volume of concrete 
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For the top width values of 0.25 and|0.30 m, the 
variation in the optimum cost is from Rs. 502 . jL 8-1 5 92 .4 2 and 
Rs. 563.29-1656,82 respectively, the optimal design being 
governed by the lower bound on the stem root thickness. 

6. 2. 4. 2 Dense sand 

The lower bound on the stem root thickness continues 
to be the active constraint for heights from 3,0 to 4.5 m, 
the top width being 0.20 m. The corresponding variation in 
the optimum cost is from Rs, 464.95 to Rs. 893.75. For heights 
above 4.5 m, the active constraint is on the safety against 
failure of the wall by sliding on the foundation. The 
optimum cost variation is from Rs. 1112.62 to Rs. 1619.25 as 
the height changes from 5.0 to 6.0 m. 

When the top width has values of 0,25 and 0.30 m, 
the optimum cost ranges from Rs. 521.56-1670.26 and Rs. 569.13- 
1694.29 respectively as the wall height varies from 3.0 to 

6.0 m. In most of these cases, the lower bound on the stem 
root thickness governs the optimal design. 

5. 2. 4. 3 Loose sand 

For a top width of 0.20 m, the optimum cost ranges 
from Rs. 544.19 to Rs. 2156.26 as the wall height varies from 

3.0 to 6.0 m. The lower bound on the stem root thickness 
behaves as an active constraint for wall heights upto 4.0 m, 
after which the active constraint is on the safety against 
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For larger heights/ a change over from the fixed 
to simple supports will substantially increase the cost but 
it may not do so for that from simple supports to a canti- 
lever. 

6.3.3 Optimization of Thin-Walled sections under Dynamic 
Constraints 

For a span of 3.0 m, as the support conditions are 
changed from fixed to simple/ there is no increase in the 
optimum weight/ when the bounds on the first natural vibra- 
tional frequency have a broad range; in case of a narrow range, 
the increase in the optimum weight is by 16.54%. For a span 
of 6.0 m, these increases are 74.90% and 23.22% respectively. 

Hence, fixed beams are significantly economical 
compared to simple beams, for larger spans and subject to a 
broad range of bounds on the natural vibrational frequency; 
in case of smaller spans, the economy achieved may be 
marginal . 

On the other hand, if the prescribed bounds have a 
narrow range, a small economy may be expected by choosing a 
fixed beam for small or large spans. 

6.3.4 Optimization of a cantilever Retaining Wall 

6. 3. 4.1 Effect of top width 

When the earth fill is gravel, for a height of 
3.0 m, the optimum cost increases substantially by 23.95% as 
the top width changes from 0.20 to 0.30 m. However, this 
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v X '' V y = transverse shear forces in the principal planes 

of a member 

V sdl /V s(:a2/ V sd3 = volume of distribution steel in the stem, toe 

and heel respectively of a cantilever retaining 
wall 


V sml'' V sm2 #V sm3 “ vo ^ ume °f main steel in the stem, toe and heel 

respectively of a cantilever retailing wall 

= volume of temperature reinforcement in a 
cantilever retaining wall 

V ss' V sh' V st = total volume of steel in the stem, heel and toe 

respectively of a cantilever retaining wall 

W = weight of cantilever retaining wall 

W 1 /W 2' * * * /W 4 = wei 9 hts °f different portions of a cantilever 

retaining wall 



X = {X > 
x imax' x imin 

y 1 -y 2 


= weight of steel in a cantilever retaining wall 
= unit weight of reinforced concrete' ■ 

= displacement component in z-direction 
= unit weight of soil 
= unit weight of steel 
= design vector 

= a design point 


★ 

= {x^} = constrained optimal point 

= maximum and minimum values of a design 
variable x^ 

= spacing of main and distribution steel 
respectively, in the stem. 
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Y 3' Y 4 

Y 5' Y 6 


z 




Px' .Py 

Yi ( Z ) 

A 

6 

6 , 6 ’ 


ii 

K 


= spacing of main and distribution steel 
respectively in the toe ,, 

= spacing of main and distribution steel 
respectively in the heel 
of a cantilever retaining wall 

= eccentricity of the foundation pressure 
. respltant with respect to the base centre 

= relative weightages assigned to various 
functions in a multiobjective function 

= cross-sectional parameters , 

= displacement modes 

= maximum difference of f friction values in the 
Powell's, method , , 

* i • i * . 

= thickness of a cross -section, .. ;; 

= obliquity of soil pressure resultants P and P 
respectively 

= displacement component in y~direction 

= kl, a constant for a thin -walled member 


X 


* * 

X , X. 

* min 
£ 

^per 

P 

a 

„(U rr (2) 

cr / cf 

max 
o . 


= coefficient used to define a convex function 
= eigenvalue 
= Lagrange multipliers 
= optimum step length 
= smallest eigenvalue 

= displacement component in x -direction 
= permissible value of deflection in a beam 
= mass per unit volume 
= normal stress at a point 

= normal stress at points (1), (2) etc. in a 
cross-section 

= maximum value of normal stress at the ith 
node 
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= critical value of normal stress in a member 

= permissible value of normal stress 

= critical value of buckling stress in a column 

= critical value of compressive stress in a 
column at working loads 

= permissible value of compressive stress in a 
column 

= permissible compressive stress in bending for 
concrete 

= permissible stress in tension for steel 
= shear stress at a point 

= shear stress at points (1)/ (2), (3) etc. in a 
cross-section 

= maximum value of shear stress at ith node 
= permissible shear stress in concrete 
= critical value of shear stress in a member 
= permissible value of shear stress 
= penalty function 

= angle .-of internal friction for soil 

= rotation of member about the longitudinal axis 

= normalized sectorial co-ordinate 

= natural frequency of transverse vibration 

= value of natural frequency correspondir^g to 
the first mode of transverse vibrations! 

= upper and lower bounds on the natural frequency 
of transverse vibration 
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Table 6.26 

Optimal sectional Dimensions 

Loading; uniformly distributed 

Intensity; 1000 kg/m ; Lateral Eccentricity; 0.1 bf^ 

(from Table 6.3) 


Span 
in m 



cant: lever 




Optimal 

dimensions 

in mm 


7 

» 

! Optimum 

1 

1 

J Active 


bf 1 

1 

t f i ; 

M 2 ; tf 2 i 

1 !- 

i 

hw ; 

i 

tw 

J weight 

! in kg/m 

1 

| constra- 
J int /s 

t 

1 .0 

51.0 

6.4 

51.0 6.4 

101 .0 

4.0 

8.31 

none 

1.5 

79.0 

6.7 

79.0 6.7 

145.0 

4.7 

13.73 

* 

*13 

2.0 

97.6 

7.2 

97.6 7.2 

194 .0 

5.3 

19.17 

913 


^Maximum normal stress in the beam 
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SYNOPSIS 


Engineering is the activity through which designs 
for material objects are produced, in the field of Civil 
Engineering, a structural engineer is confronted with the 
problem of designing a wide range of structures, such as 
massive dams, retaining walls etc. (thick structures) at 
one end of the spectrum and thin -walled component members 
like rolled steel joists, columns etc. at the other end. 

It is well-known that the primary function of a 
structure is to resist and/or transmit forces through solid 
matter. To increase the efficacy 0 £ a structure, the 
structural designer should not only work out a satisfactory 
design but also should endeavour to see that it is an 
optimal one. When a range of designs exists within a 
selected design concept, the optimal design is chosen in 
such a way that it not cnly satisfies all the limitations 
and restrictions placed on it, but it is best in sane sense, 
when compared to other possible designs in the set. 

Although, the basic principles of structural 
analysis are same for all structures, the approaches for 
thick and thin structures are at variance. Even the design 
objectives in -these two classes of structures are different. 
In other words, the set of constraints to be considered in 
the design process are to be selectively chosen. However, 
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Table 6.36 

Optimal sectional Dimensions 


Loading: uniformly distributed 


intensity: 

3000 kg/m ; 

Lateral Eccentricity: 0. 

(from 

1 bf x 

Table 6.7 





Cantilever 



Span 









in m 

Optimal dimensions 

in mm 


Optimum 

Active 








weight 

constra- 


bfi * 

tff 

j bf 2 j 

tf 

Z 2 

hw | 

tw 

in kg/m 

int/s 









* 

1.0 

90.0 

6.9 

90.0 

6.9 

175.0 

5.1 

16.76 

g 13 

1.5 

116.0 

8.7 

116.0 

8.7 

241.0 

6.0 

27.12 

g 13 

2.0 

151.8 

9.5 

151.8 

9.5 

303.0 

6.7 

38.55 

g 13 


* Maximum normal stress in the beam, 
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i 

the tools available for the structural optimization, such as 
mathematical programming, optimal synthesis etc. can be 
applied with equal ease to both these categories of structures. 

A reinforced concrete cantilever retaining wall is 
chosen to represent the thick-walled class. In such structures, 
the material cost being the main objective, optimal designs 
are obtained for minimum cost for five types of soil commonly 
met with in practice. 

As rolled steel sections (I -shape) are quite 
commonly used in Civil Engineering construction, structural 
I -sections are chosen to represent the thin -walled category. 

In such a case, weight being the criterion, optimal designs 
of beams and columns of I -section are obtained subject to 
various types of loading and boundary conditions, under the 
static, the stability and the dynamic constraints. 

Under the optimization subject to static constraints, 
the loading cases considered are all of practical signifi- 
cance. The case of concentrated loads acting at selected 
points cn a beam is considered, for example. Such type of 
loads are quite comm on in the form of reactions transferred 
on the main girder through a cross-girder system. A similar 
other case of practical importance is of a laterally eccentric 
distributed load that produces significant torsion in the 
beam. Normally a section with unequal flange widths and 
thicknesses is expected to perform better in such a loading 
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situation? however , it is interesting to note that the actual 
optimal section possess equal values of these dimensions. 

In both the problems of optimization of thick- and 
thin -walled structures, the constraints used are those 
specified by the Indian codes of Practice. In the latter 
problem, the significant feature is that, the analysis is 
carried out by the discrete element technique incorporating 
Vlasov's thin-walled beam theory. The consistent matrix 
formulation, thus generated, is quite versatile so that 
thin -walled beams, having any shape of cross-section formed 
by an assemblage of plate elements, can be analysed by this 
technique. By suitably imposing the boundary conditions at 
the intermediate nodes of the discretized structure, the 
optimization problem can be extended to include continuous 
beams as well. 

Lastly, the' discrete element analysis considers 
the effect of warping of the section, an attempt never made 
so far, for obtaining optimal designs of thin-walled sections. 
The results thus obtained are, therefore, more rational and 
could be utilized in practice. 
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1 . INTRODUCTION 


1.1 GENERAL 

A structure is an element or an assemblage of 
elements put together into a shape, to resist and/or transmit 
forces. It may also be looked upon as a system consisting of 
subsystems which in turn comprise many structural elements. 

A structure is referred to as a planar or a space 
structure depending on whether or not the centre line of the 
elements forming the structure lie in a single plane. 

Depending on the shape (structural form) and nature of external 
and/or internal influences (forces), structures are called 
rods t beams, trusses , frames , -plates shells e tc. All 
structures are invariably three dimensional, but depending on 
the relative ratios of sides, they are treated as one - , two- 
and three - dimensional continua; rods, beams, trusses, frames 
belong to the first category while plates shells etc. belong 
to the second. Whenever it is possible, three dimensional 
structures are reduced to lower dimensional category for 
purposes of analysis and design. 

1.2 INFLUENCES (FORCES) 

The environment in which a structure is situated, 
is usually idealized as a system or systems of forces. Forces 
caused by interactions between a structure and its environ- 
ment are referred to as external forces and forces between 
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different parts of an isolated system are known as internal 
forces. The forces acting on a structure may be static or 
dynamic, internal or external, fixed (deterministic) or random 
(probabilistic) etc. 

1.3 STRUCTURAL ANALYSIS 

The purpose of structural analysis is to study or 
predict the response or behaviour of a loaded structure. It 
is based on four requirements namely, the equilibrium t the 
condition of compatibility , the force- deformation relationship 
of the material used and the known boundary conditions. 

If all the foregoing conditions are satisfied, the solution is 
exact. However, a solution may yield useful information 
insofar as to give a bound on the exact solution, if a few of 
the foregoing requirements are satisfied. 

Depending on the range of mechanical behaviour 
considered, fcn analysis may be Ibir&fzdly classified as elastic 
or plastic. Elastic analysis predicts the behaviour at 
working or design loads while the plastic analysis does so at 
collapse. In this thesis, the elastic analysis is dealt with, 
exclusively. 

1.4 STRUCTURAL DESIGN 

The primary purpose of structural design is force 
transmission by solid matter with certain design objectives 
such as structural safety, rigidity , serviceability , or 
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geometric, functional and aesthetic requirements. Tradi- 
tionally , a trial design is an educated human guess based on 
intuition, common sense, professional judgement and experience. 
Certain assumptions are made, a trial design worked out, tested 
by detailed examination (analysis), modified in the light of 
the results obtained and re-examined. This iterative process 

s, 

is continued until all the prescribed conditions (constraints ) 
are satisfied. Thus, the design is treated as a process of 
trial and error, except when formulated as either a direct 
design or an optimal design problem. 

1.5 STRUCTURAL DESIGN IN CIVIL ENGINEERING PROBLEMS 

In the field of civil engineering, a structural 
engineer is confronted with the problem of designing a wide 

range of structures, such as massive dams, retaining walls etc. 

% 

(thick structures ) at one end of the spectrum and thin -walled 
component members like rolled steel joists, columns etc. 

(thin structures ) at the other end. 

Although the basic principles of structural analysis 
are same for all structures, the approaches for thick and thin 
structures are at variance. In structures like darns, retaining 
walls etc. there exists a significant interaction between the 
surrounding medium and the structure. As such, the internal 
forces between the idealized environment and the structure 
have to be considered in its design. In thin-walled structures. 
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the surrounding medium has no influence insofar as its struc- 
tural design is concerned . 

The term 'failure' of a structure may mean fracture , 
fatigue , instability or total collapse. The failure may occur 
at a localized portion in the structure (local failure) or 
the gross structure may fail (overall failure) . m a thick- 
walled structure, it is sufficient to analyse it for the overall 
failure. On the other hand, due to the 'size effect', one has 
to investigate for local failure as well in a thin -walled 
structure, which makes its analysis more involved. 

A designer should not only prepare a working design 
of a structure, but should also endeavour to see that it is 
an optimum one, keeping a certain objective in mind. Even the 
design ob 3 ectives in these classes of structures are 
different, usually. Thick-walled structures in civil engineering 
are invariably of masonry or plain or reinforced concrete. 

These are usually cast-in -situ wherein the material cost is 
bulk of the total cost. The optimal design of such structures 
would be aimed at the minimal cost. 

Thin-walled structures are usually erected from 
prefabricated metallic component members that are assembled 
at site. Transportation and handling being an important 
operation involved in such a case, one would prefer to cut 
down the weight of the structural component members to a 
minimum. Weight is, therefore, chosen as the design objective 
in such a class of structures. 
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1.6 mechanical behaviour of engineering structures 

Structures respond in different ways when subjected 
to external influences depending on their geometric charac- 
teristics and on the mechanical and other properties of the 
materials used. The mechanical properties of a material may 
be represented by determining a functional relation between 
the stresses and strains or more generally between generalized 
stresses and generalized strains. 

The actual behaviour of a structural material is 
idealized by making certain reasonable assumptions. A typical 
relation between generalized stress and generalized strain is 
shown in Figure 1.1(a). 

For most materials encountered in engineering 
practice/ it is practically impossible to describe the entire 
stress-strain curve with a single mathematical expression. 
Hence the stress-strain curve is idealized as a piecewise- 
continuous one. Depending upon the aspects that are most 
important in a given problem, only that part of the curve 
pertaining to the range of material behaviour under consid- 
eration is idealized. 

Steel and concrete are the most common materials 
used in civil engineering structures . Structural steel 
exhibits a marked plastic behaviour from the elastic one, 
the yield point separating the two. The stress-strain curve 
for such an elastic-perfectly plastic material can be 
idealized as shown in Figure 1.1(b). 




Elastic- perfectly plastic Linearly elastic Non linearly elastic Visco -elastic 

Fig-l.t Stress strain- curves* 
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The stress-strain relationship for a material in 
the elastic range may be either linear [Figure 1.1(c)] or 
non-linear [Figure 1.1(d)] and the material is then referred 
to as linearly- or non-linearly-elastic. Structural steel 
can be reasonably approximated to follow a linear elastic 
behaviour/ although strictly speaking, it is a non-linear one. 

Concrete exhibits a time -dependent behaviour and as 
such, its stress -strain characteristic changes according to 
the applied strain rate. Such behaviour is referred to as 
visco-elastic typically, represented in Figure 1.1(e). 

Concrete is heterogeneous in composition. The 
constituent crystals in the micro-structure of the structural 
effect are randomly oriented. However, for all practical 
purposes in structural analysis, both concrete and steel are 
assumed to be homogeneous , isotropic and linearly ■ elastic . 

In addition to the structural material properties, 
a structural designer has to take cognizance of character- 
istics of the surrounding medium. For example, in a structure 
like a retaining wall, this medium being the soil retained, 
such physical properties as the unit weight, porosity , 
cohesion, shear strength, angle of internal friction, 
angle of repose, permeability , bearing capacity etc. are 
to be taken into account. 

1.7 FORCE -DEFORMATION RELATIONSHIP IN A STRUCTURE 

A structure subjected to the action of external 
forces gives rise to the internal forces (stresses ) that 
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deform the structure and eventually maintain its equilibrium. 

The deformations that usually include translational and 
rotational displacements are generally very small. Hence 
the analysis of internal forces may be based on the dimensions 
of the unde formed shape without appreciable error. This 
approach to the solution of structural systems is called the 
small deflection theory e which is used in the present 
analysis. Within the limitation of small deflection/ struc- 
tural systems are usually stressed below the elastic limit of 
the material. Hence the small deflection theory is also 
known as the elastic theory . 

If the stress-strain relation of the material is 
linearly elastic# the equations representing the force 
deformation relations are also ■P%T l& Ss? functions Q f internal 
forces. Such structural systems are called linear / whose 
important characteristic is that the principle of superposition 
can be used for their solutions# as is true in the present 
case. 

1.8 STRUCTURAL OPTIMIZATION 

1.8.1 General 

Engineering is the activity through which designs 
for material objects are produced. In many design problems# 
there are several possible alternative design concepts ; for 
example# a girder for a highway bridge can be concrete or 
steel# and once material is chosen# several approaches to 
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using it are available. The goal of the design process may- 
range from providing a practicable solution to a problem where 
none is previously known/ to improving on or replacing an 
existing design. Adopting this point of view that a range of 
design exists within a selected design concept, the concept 
of the optimum design is initiated because of the presumption 
that these values are to be chosen in such a way that the 
design will be the one that satisfies all the limitations and 
restrictions placed on it and is best in sane sense. 

1.8.2 Optimization Problem 

An optimization problem is usually concerned with 
the extremization of a function (functional) of several 
variables (functions) within prescribed constraints . since 
any maximization technique can be transformed into a minimi- 
zation technique and vice versa ( duality ) * without loss of 
generality, optimization can be taken to mean minimization » 
There is no single method available for solving all optimization 
problems efficiently. Hence a number of optimization methods 
have been developed for solving different types of optimization 
problems. The optimum seeking methods are also known as 
mathematical programming techniques and are generally 
studied as a part of the operations research branch in 
mathematics, that deals with the general problem of 
optimality. 
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1.9 MATHEMATICAL PROGRAMMING TECHNIQUES 

Table 1.1 gives the various mathematical programming 
techniques along with the other well-defined areas of 
operations research. 

The mathematical programming techniques are useful 
in finding the minimum or maximum of a function of several 
variables subject to limitations (constraints ) that are 
expressed as equalities or inequalities „ The stochastic 
process techniques can be used to analyse problems which are 
described by a set of random variables having known 
■probability distribution . The statistical methods enable one 

to analyse the experimental data and build empirical models to 
obtain the most accurate representation of the physical 
situation . 

In the procedure characterized as mathematical 
programming formulations, the representation of inequality 
constraints is of critical importance since this permits the 
design to be identified as one in which not all structural 
elements are subject to limiting conditions under specified 
loads. Also, in such formulations, the orientation towards 
many variable problems fits quite well with the recent trend 
in analysis towards discrete element representations that 
require large order systems. Hence, the scope of the present 
thesis is confined to the application of one of the mathema- 
tical programming techniques using the discrete element 
analysis of the idealized structure. 
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1.10 HISTORICAL DEVELOPMENT 

The existence of optimization methods can be traced 
back to the days of Newton, Lagrange and Cauchy. The 
development of differential calculus methods of optimization 
was possible because of the contributions of Newton and 
Leibnitz to calculus. The foundations of the calculus of 
variations were laid by Bernoulli, Euler, Lagrange and 
Weirstrass . The method of optimization for constrained 
minimization problems which involves the addition of unknown 
multipliers was invented by Lagrange. Cauchy made the first 
application of the steepest descent method to solve 
unconstrained minimization problems. In spite of these early 
contributions, very little progress was made until the middle 
of the twentieth century, when the high speed digital computers 
made the implementation of the optimization procedures possible 
and stimulated further research on new methods. 

In 1S>47, Dantzig developed the simplex method for 
linear programming problems. 

The principle of optimality was established by 
Bellman in 1957, for dynamic programming problems which paved 
the way for development of the methods of constrained opti“ 
mization . 

Duff in, Zener and Peterson, in 1967, developed 
geometric programming. Avriel and Williams, in 1970 > extended 
this method to include any rational function of posynomial 
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terms and called the method as complementary geometric 
programming. 

Gomory, in I960, did pioneering work in integer 
programming , which is one of the most exciting and rapidly 
developing areas of optimization. The reason for this is 
that most of the real world applications fall under this 
category of problems. 

Dantzig, in 1955, and charnes and Cooper, in 1959, 
developed stochastic programming techniques and solved 
problems by assuming design parameters to be independent and 
normally distributed. 

Network, analysis methods are essentially management 
control techniques and were developed during the years 1957 
and 1958. 

The foundations of game theory were developed by 
von Newmann in 1928 and since then it has been applied to 
solve several mathematical, economics and military problems. 
Only during the last few years, game theory has been applied 
to solve some of the engineering design problems. 

1.11 NON-LINEAR PROGRAMMING 

If the optimization problem involves any of the 
objective or constraint functions as non-linear, the problem 
is called a non-linear programming problem. This is the 
most general programming problem and all other problems can 
be considered as special cases of this. Depending upon 
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whether the constraints exist or not in the problem/ it is 
classified as a cons trained or an unconstrained one. 

Mathematical programming was first applied to 
structural optimization in the late 1950's. According to 
early contributions by Livesley (1956) and Pearson (1958)/ 
limit design was treated as a linear programming problem/ 
while Schmit (1960) cast the elastic design as the more 
general non-linear programming problem. 

1.11.1 Unconstrained Optimization 

Major developments in the area of numerical methods 
of unconstrained optimization were made in the United Kingdom 
•in the early 1960’s. 

Under one such category of methods called direct 
search methods / unconstrained minimization problems were 
initially handled by the random search and later by the 
univariate methods. 

The simplex methods were introduced by Spendley, 
Hext and Himmsworth (1962). The original simplex triangle 
was modified and more effectively used by Nelder and Mead 
(1965) and by Box, Davies and Swann (1969). 

Hooke and Jeeves (1961) initiated the pattern search 
methods. Rosenbrock (1960) independently introduced the 
method of rotating co-ordinates t which can be considered 
as the development of the former one. Powell (1964) extended 
the basic pattern search method. It is a very powerful and 
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the most widely accepted direct search method. Powell has 
suggested some modifications to facilitate the convergence of 
this method when applied to non -quadratic objective functions , 
which have made it more versatile. 

Among other category of the unconstrained minimi- 
zation methods/ come the descent methods. One such type is 
the steepest descent method pioneered by Cauchy in the middle 
of nineteenth century. The method was used with moderate 
success on a wide variety of problems through the middle 1950's. 

Fletcher and Reeves (1964) greatly improved the 
convergence characteristics by modifying it into a conjugate 
gradient method. 

The oldest method for solving a set of non-linear 
equations is the Newton's method. Based on this, a class of 
methods called quasi- New ton methods was developed, since 
these can be regarded as approximations to the Newton's 
method in some sense. One such method is th q- variable metric 
method introduced first by Davidon (1959) and developed later 
on by Fletcher and Powell (1963). 

1.11.2 Constrained Optimization 

Many techniques have been developed for the 
solution of a constrained non-linear programming problem. 

All these techniques can be classified into two broad cate- 
gories namely, the direct methods and the indirect methods. 

In the direct methods, the constraints are handled in an 
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explicit manner whereas in most of the indirect methods, the 
constrained problem is solved as a sequence of unconstrained 
minimization pr obi em s . 

Among the direct methods , one such approach is by 
approximating the constraints 0 Kelly (1960) used the cutting 
plans method for solving convex programming problems. 

Another- method developed by Zoutendijk (1960) makes 
a successive use of feasible directions in obtaining the 
constrained minimum. 

The most significant under the indirect methods are 
the penalty function methods. These are further categorized 
into the interior penalty and the exterior penalty function 
methods. They were successfully developed by carrol (1961), 
Zangwill (1967) and Fiacco and McCormick (1963, 1964, 1968). 

1.12 CIVIL ENGINEERING APPLICATIONS OF OPTIMIZATION 

Optimization, in its broadest -sense, can be applied 
to solve any engineering problem. In the field of civil 
engineering, optimization studies are useful to 

(a) develop a fully automated design with the aid of high 
speed digital computers of structures like frames, 
foundations, bridges, towers, dams, retaining walls etc. 
for minimum cost. 

(b) Obtain minimum weight design of thin -walled structures 
for earthquake, wind and other types of random loading. 
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(c) develop optimal plastic design of structures. 

(d) select a suitable configuration and topography for a 
water resources system and to design it for the maximum 
efficiency. 


1.13 STATEMENT OF AN OPTIMIZATION PROBLEM 


An optimization or a mathematical programming 
problem can be stated as follows 


f 


*1 


find X 


| X 2 \ 

j . v which minimizes f(X) 

V f 

n j 

4 


subject to the constraints 


1.1 


(X) <0/ j = 1, 2, ..., m 

and, lj (X) = O, j = m+l, m+2 , ..., p . 

where X is an n-dimensional vector called the design vector _ 
f(X) is called the objective function and g^(X) and 1^ (X) 
are respectively the inequality and the equality constraints. 

The number of variables 'n' and the number of 
constraints 'p' need not be related in any way. Thus, in a 
given mathematical programming problem, p could be less than, 
equal to or greater than n. In some problems, the value of p 
might be zero which means that there are no constraints in 



the problem. Such type of problems are called unconstrained 
optimization problems. Those problems for which p is not 
equal to zero are known as constrained optimization problems. 


1.13.1 Design vector 

Any engineering system or component is described by 
a set of quantities , some of which are treated as variables 
during the design process. In general, certain quantities are 
usually fixed at the outset and these are called p reassigned 
parameters . All the other quantities are looked upon as 
variables in the design process and are called design or 
decision variables x j _, i =1, 2, ..., n. The design variables 
are collectively represented as a design vector 


( X 1 

! 


X 


{x ± } 


nj 


The design variables of an optimum structural design 


problem may consist of the member sizes, parameters tnat 
describe the structural configuration and the mechanical or 
physical properties of the material as well as other quanti=» 
fiable aspects of design. The topology (i.e. the pattern of 
connection of members in the mathematical model) of the 
complete structure is difficult to take into account although 
it is treated to a limited extent when the algorithm employed 
permits members to reach zero size. 
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If an n~dimensional cartesian space with each co- 
ordinate axis representing a design variable x^, i = 1/ 2, 

. .., n, is considered, the space is called the design variable 
space. Each point in the n -dimensional design space is called 
a design point and it represents either a possible or an 
impossible solution to the design problem. 

1.13.2 Design Constraints 

In many practical problems, the design variables can-“ 
not be chosen arbitrarily; rather they have to satisfy certain 
specified functional and other requirements. The prescribed 
restrictions that must be satisfied in order to produce an 
acceptable design are collectively called design constraints . 

A constraint may take the form of a limitation imposed 
directly on a variable or group of variables (explicit 
constraint ) or may represent a limitation on quantities whose 
dependence on the design variables cannot be stated directly 
(implicit constraint ). 

The constraints which represent limitations on the 
behaviour or performance of the system are termed as b ehavioural 
or functional constraints. The constraints which represent 
physical limitations on the design variables like availability, 
f abricability and transportability are known as geometric or 
side constraints. The behavioural constraints in structural 
design are usually limitations on stresses or displacements. 

But they may also take the form of restrictions on such factors 
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as vibrational frequency or buckling strength. Behaviour 
constraints are generally implicit, although typically 
explicit behaviour constraints are given by formulae presented 
in design specifications. 

In general, the given set of constraints may include 
the equality as well as inequality constraints. In theory, 
each equality constraint is an opportunity to remove a design 
variable from the optimization process and thereby reduce the 
number of dimensions of the problem. However, as the elimi- 
nation procedure may be awkward and algebraically complicated, 
this approach is not always adopted. 

The idea of an inequality constraint is of major 
importance in optimum structural design. If equality const- 
raints only were stipulated in a design limited by stresses 
alone, all procedures would lead to fully stressed design . 

It has been proved that fully stressed designs are not 
necessarily minimum weight designs and for optimality, it 
is essential to permit designs in which not all stress cons- 
traints are satisfied identically i.e. inequality constraints. 

In general, constraint equations are non-linear 
functions of the design variables when the behaviour of the 
structure is elastic and statically indeterminate although 
the analysis technology is linear . 
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1.13.3 constraint Surface 

Each constraint condition appears in design space 
as a surface representing the locus of design points which 
cause the constraint to be satisfied as an equality constraint. 

For illustration, consider an optimization problem 
with only inequality constraints g^(X) ^ 0. The set of values 
of X that satisfy the equation g^(X) = 0 forms a hypersurf aoe 
in the design space and is called a constraint surface. 

It may be noted that this is an (n-1 ) “dimensional 
subspace, n being the number of design variables. The const- 
raint surface divides the design space into two regions; one 
in which g^ (X) <0 and the other in which g^ (X) > 0. Thus 

the points lying on the hypersurface will satisfy the const- 
raints gj(X) critically ■, whereas the points lying in the 
region where g^Cx) > 0 are infeasible or unacceptable s and 
the points lying in the region where g^(X) < 0 are feasible 

or acceptable . The collection of all the constraint surfaces 
gj(X) =0, j =1, 2, . m which separates the acceptable 
region is called the composite constraint surface. For 
continuous design variables, the constraint surface is in 
general continuous s and is curved for the elastic structural 
design of an indeterminate structure. 

Figure 1.2 shows a hypothetical two-dimensional 
design space where the infeasible region i* indicated by 
hatched lines. A design point which lies on one or more than 
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one constraint surface is called a bound point and the asso- 
ciated constraint is called an active constraint. The design 
points which do not lie on any constraint surface are known 
as free points. Depending upon whether a particular design 
point belongs to the acceptable or unacceptable region/ it 
can be identified as one of the following four types: 

(a) free and acceptable point/ 

(b) free and unacceptable point, 

(c) bound and acceptable point, and 

(d) bound and unacceptable point. 

All these four types of points are shown in Figure 

1 . 2 . 

1.13.4 Objective Function 

The conventional, design procedures aim at finding 
an acceptable or adequate design, which merely satisfies the 
functional and other requirements of the problem. In general, 
there will be more than one acceptable design and the purpose 
of optimization is to choose the best one out of the many 
acceptable designs available. Thus a criterion has to be 
chosen that constitutes a basis for the selection of one of 
several alternative acceptable designs. The criterion with 
respect to which the design is optimized, i.e. whose least 
(or greatest) value is sought in an optimization procedure, 
when expressed as a scalar function of the design variables, 
is termed the cost function 9 merit function or objective 
function. 
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.The objective function represents the most 
important single property of a design and its choice is 
governed by the nature of a problem. The objective function 
for minimization is generally taken as weight in aircraft; and 
aerospace structural design problems, in civil engineering 
structural designs, it is usually the minimisation of cost. 

The maximization of mechanical efficiency would be the 
obvious choice in mechanical engineering systems design 
problems . 

In some situations, there may be more than one 
criterion to be satisfied simultaneously. One way to handle 
such a problem is to take the actual objective function as a 
weighted sum of the multiple oh jective functions . Thus if 
f r (X) and f 2 (X) are the two objective functions possible, the 
new objective function f(X) may be constructed for optimization 
as 

f(X) = G 1 f 1 (X) + o 2 f 2 (X) ... 1.2 

Of, toeing the relative weightages assigned. 

1.13.5 Objective Function Surfaces 

The locus of all points satisfying f (X-) = C = 
constant, forms a hypersurface in the design space and for 
each value of 'C' there corresponds a different member of a 
family of surfaces. These surfaces called the objective 
function surfaces are shown in a hypothetical two-dimensional 
design space in Figure 1.3. 
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Once the objective function surfaces are drawn along 
with the constraint surfaces/ the optimum point can be deter- 
mined without much difficulty. However/ as the number of 
design variables exceeds two or three/ the constraint and 
objective function surfaces become complex even for visualiz- 
ation and the problem has to be solved purely as a mathematical 
problem . 


1.13.6 Gradient 

An important concept employed continually in the 
chapters that follow is that of the gradient of the objective 
function/ vf . The gradient is a vector composed of the deri- 
vatives of f with respect to each of the design variables, i.e. 


Vf = 


8f 

"3X^ 

8 f 
9 x 0 
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f 


... 1.3(a) 


n J 


Similarly the gradient g^ for the jth constraint surface is 
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... 1.3(b) 
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1.13.7 Local and Global Minima 

A structural designer, while applying the mathema- 
tical programming methods , should be able to distinguish 
between local minima and the global minimum [Figure 1.4(a)] . 
If, for example, the optimization procedure leads from the 
starting point A to the point B , analytical tests applied at 
point B will indicate that no further moves are possible 
without a violation of constraints. Thus a local minimum 
will have been reached. The absolute or global minimum, is 
seen to exist a point C, however. 

A local minimum coincides with the global minimum 
when the constraint surface is convex. Mathematically, a 
function f(x) is said to be convex for any pair of points 



and all X, 0 ,< X < 1 , 

f [X.X 2 + (1 - X)^] < Xf(X 2 ) + (1 - X) fCXj.) ... 1.4 

Figure 1.4(b) illustrates a convex function in one 
dimension. It can be easily seen that a convex function is 
always bending upwards and hence it is apparent that the 
local minimum of a convex function is also a global minimum. 



Side constraint 



(a) Local and global minima (b) convex function . 
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Thus, the property of convexity is of value in ascertaining the 
uniqueness of a solution for an optimum design point. 

Unf ortunately, constraint surfaces in structural 
design problems are not always convex. Kavlie and Moe (1971), 
Adidam S.R. (1972), and Moses and Onada (1969) have given 
examples of structural design problems with concave constraint 
surfaces . 


1.13.8 Kuhn-Tucker Optimality Condition 

The analytical test for a local minimum alluded to 

above is the Kuhn-Tucker optimality condition after the men 

who developed them in 1951. The mathematical statement of 

this condition is as follows. 

In a constrained optimization problem, consider the 

division of the constraints into two subsets and J where 

indicates the indices of those constraints which are 

active at the optimum point and includes the indices of 

the remaining (inactive) constraints. The Kuhn-Tucker opti- 

mality conditions to be satisfied at a constrained minimum 
★ 

point X can be expressed as 


af , r , 

+ i Aj * 

3X i j J 9X 


0, i— 1, 2, 


and 


A . > 0. 

J 


n 


1 5 

• • • J. ^ w» 


In other words, these are the necessary conditions for a 
local optimum design point, but in general not sufficient. 
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However, there is a class of problems, called convex 
programming problems, wherein both the objective and the 
constraint functions are convex, the Kuhn -Tucker conditions 
are necessary and sufficient for a global minimum. 

Rewriting Equation 1.5 in the form 

vf(x ) = - . i x. ?g.(x ) ... 1.6 

j £ J 

will also mean that in order to satisfy the Kuhn-Tucker 
optimality conditions, the gradient of the objective function 
must be expressible as the negative of a linear combination 
of the active constraint gradients . The non -negative multi- 
pliers Xj (j eJ^) are also termed the Lagrange multipliers 
since the development of Equation 1.5 may be accomplished by 
the application of the Lagrange -multiplier approach to the 
constrained minimization . 

1.14 SCOPE OF THE THESIS 

The thesis is primarily concerned with the appli- 
cation of a suitable mathematical programming technique to 
obtain optimal designs of two different classes of civil 
engineering structures. A reinforced concrete cantilever 
retaining wall is chosen to represent the thick-walled 
category. In such structures, the material cost being the 
main objective, optimal designs are obtained for minimum cost 
for five types of soil commonly met with in practice. As 
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rolled steel section (I -shape) are quite commonly used in 
construction, structural I -sections are selected to represent 
the thin-walled class, in such a case, weight being the 
criterion, optimal designs of beams and columns of I-section 
are obtained for various types of loading under different sets 
of constraints viz. the static, the stability and the dynamic. 

Both these are essentially constrained optimization 
problems involving non-linear programming. Two different 
non-linear programming techniques are employed in the respec- 
tive problems — in the former, Powell's method is used while 
in the latter, the solution is attempted using Davidon- 
Fletcher-Powell * s method. In both the problems, the const- 
raints used are those specified by the Indian Codes of Practice. 

The thesis is divided into six chapters, chapter 1 
introduces the fundamental concepts of structural analysis 
and structural design which can be found in any standard 
textbook on structural theory. The critical aspects that 
govern categorization of the civil engineering structures 
into the thick- and thin -walled class are then discussed. A 
brief historical background of different mathematical 
programming methods is followed by a selective literature 
survey of the various non-linear programming techniques. 
Summarizing, in chapter 1, it is intended to develop a 
proper background for the formulation and subsequent solution 
of the optimization problems alluded before. 
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in the optimization of the thin -walled class, the 
important feature is that, the analysis is carried out by the 
discrete element technique incorporating Vlasov's thin-walled 
beam theory. The analysis, therefore, considers the effect of 
warping of the section, an attempt never made so far, for 
obtaining optimal designs of thin -walled sections. Chapter 2, 
therefore, discusses the consistent matrix formulation for the 
discrete element technique for linear and eigenvalue problems 
of structures assembled from thin-walled segments with open 
cross-section. The formulation is exact for linear static 
problems and a close upper bound for eigenvalue problems. 

The technique presented may be considered as an extension of 
the procedure used for solid beam structures. 

In optimum structural design, the behavioural 
constraints are usually the prescribed limitations on stresses 
and deflections. For such a set of static constraints, the 
formulation of the optimization problem for thin-walled 
sections is discussed in chapter 3. 

Thin-walled component members are widely used in 
aerospace structures wherein the payload plays a significant 
part in creating the structural and/or the dynamic instabi- 
lity in the members. In the optimal design of such sections, 
the behavioural constraints may be in the form of the 
prescribed buckling strength or the specified range of 
natural vibrational frequency respectively. Optimal design 
of thin -walled sections under such sets of the stability or 
the dynamic constraints are dealt with in Chapter 4. 
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An important reason to group these two sets together 
is that the problem of stability and the problem of free 
transverse vibrations, both come under the category of the 
eigenvalue problems. Mathematically, an eigenvalue problem 
in the general form consists of determining values of A 
for which the ' n' homogeneous linear equations in *n' 
unknowns 


[A] X = A . [B] X 


... 1.7 


have a non- trivial solution X such that 


X 



and ^ 0 for all i, 


The (n x n) matrices IaI and [B] in Equation 1.7 take the 
form of the stiffness ■ and the geometric matrices respect - 
ively in the stability problem, while in the free transverse 
vibrations problem, matrix [B] is replaced by the mass matrix. 

Chapter 5 discusses the formulation for obtaining 
optimal designs of cantilever retaining walls for different 
classes of soil retained. The behavioural constraints take 
into account different modes of failure of the wall as well 
as the specifications governing the reinforced concrete 
design. 
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In chapter 6, the results obtained in the preceding 
chapters have been critically evaluated and conclusions 
drawn. The thesis concludes with the presentation of scope 
for further studies. 


- o - 
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2. DISCRETE ELEMENT ANALYSIS OF THIN -WALLED BEAMS 

2 . 1 GENERAL 

According to their spatial character/ structural 
elements can be classified into four classes viz. the massive 
bodies „ plates and shells 3 solid beams snd thin -walled beams 

Massive bodies are those whose three dimensions are 
comparable e.g. a cube, a sphere, an ellipsoid etc. These 
may be looked upon as an elastic continuum filling all space 
so that the problem of determination of stresses and defor- 
mations can be dealt by the mathematical theory of elasticity. 

Bodies having one dimension (thickness) small 
compared with the other two (length and breadth) which are of 
the same order, can be treated as plates and shells. Thin 
slabs -T-orouod, rectangular, trapezoidal etc. in shape belong 
to the former category while thin shells of various shapes — 
cylindrical, conical, spherical etc. belong to the latter one. 
The general theory of such bodies is based on geometrical 
hypotheses valid for thin deformable bodies within a certain 
degree of approximation . 

2.1.1 Solid Beams 

Solid beams are bodies characterized such that two 
of their dimensions, usually the cross-sectional ones, are 
small compared with the third dimension, the length. Because 
of this particular form assumed, the analysis of solid beams 
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involves a number of geometrical hypotheses to simplify the 
calculations. This may be explained with reference to the 
following two kinds of basic problems one comes across in the 
strength of materials . 

In the problem of transverse bending , out of the 
six components of the general strain tensor * only one corres" 
ponding to the direction parallel to beam axis is preserved 
and other five components are taken to vanish, such an 
assumption is possible because of the law of -plane sections 
which means that any cross-section of beam which is initially 
plane, remains plane after deformation. This geometrical 
hypothesis, postulated by Navier, forms the basis of the 
elementary theory of bending of beams. 

The problem of torsion follows St. venant's theory of 
pure torsion i.e. no extension or shearing strain is assumed 
to exist in the plane of cross-section. In other words, a 
plane section remains plane i.e. free from any warpage or 
distortion , ■ after deformation . This theory allows to deter- 
mine the tangential stresses in the cross-section of beam. 

In general, the class of 'solid beams’ includes all 
beam systems ~ plane or three-dimensional, statically deter- 
minate or indeterminate whose elements are beams following 
the laws mentioned above. 

2.1.2 Thin -Walled Beams 


These are bodies characterized by the fact that 
their three dimensions are all of different magnitude. They 
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can be looked upon as essentially having the form of a long 
prismatic shell such that the thickness of cross-section (6) 

<< the characteristic dimension of cross-section (d) << the 
length of the shell (1). 

Rolled, welded or rivetted structural steel sections 
used as columns or beams, elements of girders and frames etc. 
wherein the ratios 6/d and d/1 may be of the or dev of 0.1 or 
less, can be categorized under thin-walled beams for all 
practical purposes, certain engineering structures do have 
such proportions that they can be considered as thin-walled 
structures e.g. certain types of girder and arch bridges with 
sufficiently rigid cross-section, suspension bridges of a 
roadway with a channel or 1-section girder system, long rein- 
forced concrete ribbed structures, cylindrical or prismatic 
shell arches, bunkers, pipe lines etc. 

2.2 DISTINCTIVE FEATURES OF THIN -WALLED BEAMS 

When subjected to torsion, thin-walled beams undergo 
longitudinal extension as a result of relative warping of the 
section, which does not remain plane. Hence, complementary 
longitudinal normal stresses proportional to these strains are 
created which have to maintain the internal equilibrium of 
longitudinal forces in each cross-section. These are not 
examined in the theory of pure torsion of solid beams , but 
can attain large values in thin-walled beams with open (rigid 
or flexible) cross-section and beams with closed (flexible) 
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cross-section. Thin-walled beams with closed (rigid) cross- 
section may be considered as solid beams so far as their 
behaviour under combined flexure and torsion is considered. 
Complementary longitudinal stresses that occur in such beams have 
a local character. By St. Venant's principle, they fall off 
rapidly along the beam. 

2.3 VLASOV'S THEORY . 

In the theory of thin -walled beams of open section, 
Vlasov (1959) has postulated the following two assumptions: 

(a) A thin-walled beam of open section can be considered 
as a shell of rigid (i.e. undef ormable) section, 
and, (b) Shearing deformation of middle surface of the beam 
cross-section can be assumed to vanish. 

These can be so generalized that a loading state at 
a section can be considered as superposition of two independent 
states and the total deformation can be obtained by super- 
position of deformation in each case. 

This can be illustrated referring to Figure 2.1(a) 
wherein the original loading is viewed as a statical equivalent 
of two independent states shown in Figures 2.1(b) and (c). 
Obviously, Figure 2.1(b) is a state of flexure that can be 
explained by the law of plane sections . Figure 2.1(c) can be 
looked upon as a typical of restrained torsion i.e. torsion 
in which independent longitudinal elements of the beam undergo 
bending in addition to torsion. This is referred to as 
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'flexural torsion that suggests the existence of normal 
stresses due to flexure in a cross-section in addition to 
tangential stresses due to torsion* As a result , the cross- 
section does not remain plane — it is warped and hence this 
kind of torsion is also called the warping torsion . 

In the case of a thin-walled beam subjected to 
torsion, the subsequent presence of the longitudinal stresses 
infers that, part of the work done by the applied twisting 
moment (T) is used in developing the warping torsion (Tq), 
while the remainder will develop shear stresses associated 
with the pure or the St. Venant torsion (t ) . Thus 

T = T n + T ... 2.1 

ft s 

When the warping torsion produces normal stresses 
at each section of a thin-walled beam, their stress resultant 
is developed as a generalized force that has the form of a 
system of self -balancing longitudinal tensions . Such resulting 
pair of bending moments are referred to as longitudinal 
bimoment of the beam. The bimoment has the dimension of 
force multiplied by area. 

2.4 ANALYSIS OF THICK- AND THIN -WALLED BEAMS 

TO calculate the stresses and strains at a given 
section of a thick-walled beam, whose behaviour is governed 
by the law of plane sections s one has to know an internal 
axial force and bending moment together with the area and 
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moments of inertia of the cross-section. To calculate the 
deflection of such a beam, one must also know the distribution 
of external forces along the length of the beam. 

All these values are not sufficient for the analysis 
of a thin -walled beam, because the hypothesis of plane sections 
is no more valid . 

Hence Vlasov introduced two additional , new types 
of forces termed the flexural twist and the bimoment . He 
also introduced additional functions of the sectorial 
properties of a section, similar to first and second moments 
of area. The introduction of these new forces and properties 
of a section enables a complete analysis for determining the 
stresses and distortions of thin -walled beams. 

2.5 CONSISTENT MATRIX FORMULATION OF THE PROBLEM 

In the problem of optimal structural design of 
rolled steel i-shaped sections, the behavioural constraints 
have been chosen to be limitations on stresses and deflections 
(’static 9 constraints) as well as prescribed buckling 
strength (stability constraints) and specified range of 
natural transverse vibrational frequency (dynamic constraints) . 
This, being a multivariate problem, fits quite well with the 
recent trend in analysis towards dvscrete element represen- 
tations. The structure under optimization can be looked upon 
as an assemblage of thin-walled discrete elements with open 
cross-section. Hence, it is essential to establish a 
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consistent matrix formulation for static problems, elastic 
stability and dynamic response of a structure assembled from 
thin-walled members of any open cross-section in general. An 
I -shaped section being a particular case of an open section, 
the foregoing results can be readily used. 

Hitherto the theory of flexure and torsion of thin- 
walled members, having as a distinctive feature the occurrence 
of normal stresses as a result of torsion (due to warping) has 
been elaborated. However, as far as the structural systems 
assembled from such members are concerned, very little has 
been done to treat the problem on the lines of modem tech- 
niques. Since the single thin-walled member is by itself 
statically indeterminate regardless of boundary conditions, 
the number of redundant forces is considerably higher than 
for a similar structure assembled from solid beams. This 
undoubtedly emphasises the need for the development of a method 
suitable for the application of high speed digital computers . 

Krahula (1967) derived the stiffness matrix for a 
thin -walled member. One has to mention also the paper by 
Renton (1967) considering one of the aspects of stability. 

Both these papers however lacked generality and consistency 
of corresponding procedures developed for solid beam structures 
by Archer (1963, 1965) and Hartz (1965). Kollbrunner (1968) 
developed a displacement method and a new calculation model 
for the thin -walled bars. Krajcinovic (1969) developed a 
consistent and general matrix formulation of thin -walled beams 
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applicable to all linear and linearized structural problems. 
Recently, Nitssche (1976) developed stiffness matrices for 
thin-walled beam members of closed section subjected to 

flexure and torsion, considering the effects of warping torsion. 

* 

2.6 FORMULATION OF THE PROBLEM 

A thin -walled member, as a basic element of a 
structure articulated from many such components, is considered 
as a spatial system composed of plates undergoing both 
bending and axial stressing in plane. The general theory of 
thin-walled members of open cross-section, as developed in 
its final form by Vlasov (1959), is essentially used in the 
formulation. The two basic assumptions on which the theory 
is based are enumerated in Article 2.3 and these reduce the 
number of degrees of freedom at each cross-section to four. 

Thin -walled beams treated according to such a theory are 
distinguished from solid beams by experiencing longitudinal 
strains as a result of torsion (due to warping). In other 
words, instead of Navier's hypothesis of plane sections, a 
more general hypothesis governing the kinematics of the 
member is introduced . 

Consider a member having a thin -walled, open cross- 
section, as shown in Figure 2.2. The section being rigid, 
the displacement of a point on the middle surface of the 
member can be described by three component displacements K , 

T) and w Q in the direction of the principal axes x, y and z 
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where d ± and 6 are the breadth and thickness of the i th 
plate; I is the sectorial moment of inertia defined by 

I nn = / Q 2 dA ... 2.3(c) 

where is the normalized sectional co-ordinate; and 

T = T n + T ...2.1 

it s 


The total torque 'T' is a vector sum of the St. venant's (T ) 

s 

and the warping torque (Tj as already indicated in Equation 

2.1. The primes denote differentiation with respect to z. 

It should be noted that Equations 2.2 imply the 

orthogonality of co-ordinates x, y and fi or that in other 

words, the normal force N and beading moments M and M are 

z yc y 

reduced to the centre of gravity G, while the torque T and 

transverse forces V and v are reduced at the shear centre S 

x y 

(Figure 2.2). It may also be noted that the transverse forces 

V and v cannot be defined in terms of deformations as a 
x y 

result of imposed assumptions about the deformation. 

The equations of equilibrium in terms of componental 
deformations for the linear static problem are due to the 
orthogonality of x , y and & uncoupled. 


z -direction, 
x-direction. 


EAw 


EI xx? 


= p + m 

1 X X 


EI yy 11 


= p + m 

p y y 


... 2.4(a) 
... 2.4(b) 


y-direction 


2.4(c) 
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Moment equilibrium 

about z-direction, EI QQ ^ G K 9 


I 

+ m fi ... 2.4(d) 


where p , p and p are distributed loads the directions 
X Y 2 

of the x, y and z axes respectively while m , m , m and 

!X y is 

are distributed bending moments/ torque and bimoment respect- 
ively acting upon the member. 

Equation 2.4(a) governs the problem of an axially 
stressed member while Equations 2.4(b) and (c) govern the 
problem of flexure in the two principal planes . Equation 2.4(d) 
governs the problem of warping torsion of a thin-walled member. 

Considering the order of differential equations in 
relations 2.4(a) through (d), a set of fourteen boundary 
conditions (seven for each terminal cross-section of the 
member) should be attached to this system in order to define 
the boundary value problem. Such seven possible displacements, 
called the nodal degrees of freedom for each terminal cross- 
section of the element are; K, T), w / the displacements 

I I 

along axes x, y and z; 5 , T) and <p, the rotations about 

I 

these axes, and 9 , the warping of the section. The solution 
of the problem is attempted by assuming a displacement mode. 


2.7 DISPLACEMENT FUNCTIONS 

Define the displacement field q(z, t) as 
q(z, t) = 1 q ± (t) Y-j/z) 


q l Y l + q 2 T 2 + 


• « « 


+ Vn 
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where Y^(z) are some assumed displacement modes; q(t) are 
the unknown amplitudes at the structural nodes i.e. the 
nodal dis -placements and n is the total number of generalized 
nodal displacements, so chosen that it equals the nodal 
degrees of freedom for an element. 


Moreover, the vector {q^} = 


! q l I 


/ 

\ 


1 


can be 


Y 9 " 


l 

J 


subdivided into four subvectors 



z 

The first two generalized displacements (subvector q ), 

characterize the axial stressing of the member. These, in 

sequel, be neglected as the axial deformations are expected 

to be of the small order. Modes y^ 2 ) corresponding to next 

x y 

eight nodal displacements q^ (subvectors q and q ) charac- 
terize the flexure of beam in the two principal planes, 
while the remaining four (subvector q^ ) characterize its 
flexural torsion. The corresponding twelve nodal displace- 
ments characterizing the twelve displacement modes, are 
shown in Figures 2.3 and 2.4. • 
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2.7.1 f lexur e 


For flexure in the principal plane x-z, the first 

four displacement functions can be obtained with reference to 

% 

Figure 2.5. The displacement mode Y^(z) is the equation of 
the elastic curve corresponding to q 1 =1 and all other 
displacements remaining zero [Figure 2.5(a)] . The governing 
Bernoulli »Euler equation for flexure is given by 


d£ 

El — V = 0 

** az 4 


? 7 

• • * • / 


Successive integration leads to 

z 3 z 2 

£ = c x ~ + c 2 ~ + c 3 z + c 4 ... 2.8 

where c-j_ / c c 3 and c^ are the constants of integration. 
Using the boundary conditions 

£ = 1, £ 1 = 0, at z = 0 

2 9 

• » • ^ 

i 

and £ = £1=0, at z = 1 


the displacement function for the first mode is obtained as 


Y x (z) 


(1 




... 2.10(a) 


Similarly the second ^ third and fourth modes in flexure are 


Y 2 (z) 



... 2.10(b) 


Y 3 (z) 



UjfcjfK' # 10(c) 

02234 » 

Jjp% 


5 ,5 S¥ ^ 




Fig . 2 .5 Displacement 
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and Y 4 (z) = \ (1 “ f> * ... 2.10(d) 

The foregoing Y's are the set of displacement modes for 
flexure in the principal plane x=z. 

Proceeding on similar lines, corresponding set of 
displacement modes for flexure in the principal plane y-=z is 

Y 5 (z) 

V 2) 

Y 7 (z) 

and Yg( z ) 

2.7.2 Flexural Torsion 

Starting from Equation 2.4(d), the solution of the 
homogeneous part is 

<p(z, t) = CgChkz + c^shkz + c^kz + Cg ... 2.12 

where k 2 = GK/EI^q ... 2.13 

Again, Yj_ i- s the solution to the homogeneous equation 

for g. =1 and all other q. =0 (i = 9, ..., 12; j = 1, 2,..., 
i J 

12; i ^ j)/ so that 


(1 + ~f)(l ~ f) , 


- (1 ™ -) 

1 V V ' 


... 2.11 


z_ (3 „ 

l 2 ( l} ' 


— (1 - -) 
2 V l } 

1 
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Y 9 ( 2 ) = "Jg" [ ( 1 - ch k ) ch kz + sh < sh kz 

= kz sh k +1 ~ ch k + < sh <] 

Y 10 (z) = [(k ch k ~ sh k) ch kz + ( ch k - 1 

k1 D 

- k sh k ) sh kz + kz ( ch k « 1 ) + sh k < ch<] 

... 2.14 

Yu (z) = ~~ [(ch k - 1 ) ch kz - sh k sh kz 

+ kz sh k » 1 = ch <_ ] 

and 

1 

Y 10 (z) = — ~-f~- [(sh k - k) ch kz + (1 - ch K ) sh kz 

K 1 Z D 

+ kz + kz ( ch k =■ 1 ) + k •= sh k ] 

where 

D = 2(1 = ch k) + k sh k ...2.15 

and < = kl ...2.16 

with ' 1' denoting the length of the member. The symbols 

'ch* and 'sh 1 refer to the hyperbolic cosine and sine functions 

respect ively. 

For the generalized nodal displacements and the 
corresponding generalized nodal forces shewn in Figures 2.3 
and 2.4, for the element under consideration. 
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the nodal displacement vector {q^ } = [q^, q 2 , .../ q^ 2 ] 


^ 5 a' ^a 1, k b' ^b 1, ^a' 11 a 1 ' ^b' ^b 1 ' ^a 1 ' ? a i2/ ^b 1 ' ^b 3-2 ^ T 


.. 2.17 


and the generalized nodal force vector {Q i > = [ Q 1 , Q 2 ,...,Q 12 ] 

= Ka' V*' V xB' M yB A ' V' M XA 71 ' V' M XB 7 ^ V 1 ' 




2nT 


B ' flB' 


... 2.18 


The potential energy may then be written as a scalar 


product 


V = (Q.. } ... 2.19 

Some displacements and forces are multiplied by the length of 
the beam 1 to various powers in order to achieve dimensional 
homogeneity. The nodal displacements q^ have the dimension 
of length while the generalized nodal forces have the 
dimension of a force. 


2.8 NODAL FORCE VECTOR 

The nodal force vector {q^} has to be determined at 
a number of nodes (structural co-ordinates ) that characterize 
the division of the structure ( discretization ) into a finite 
number of elements , and then it is assembled for the entire 
structure. The load applied to the structure can, in general, 
be either distributed or concentrated at certain locations 
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(not necessarily coinciding with the nodes) or both. It is 
therefore necessary to establish a relation between the nodal 
force vector {Q^} and external loads, which will be consistent 
with the structural idealization . in order to achieve this, 
one can obtain the equivalent nodal loads vector {F^} by 
considering the strain energy equivalence of the deformed 
structure with the original external loads. The nodal force 
vector {q.} can then be determined as the statical equivalent 
to the vector {F^} . 

For a flexure problem, one can obtain the nodal 
loads vector {F^} for any type of loading by using necessary 
relations given in Martin (1966) or Gere (1969) and thus can 
obtain {Q^}. 

Recalling the basic principles of the theory of 
structures and in particular the Mueller-Breslau theorem, 
one can define the displacement modes as the influence functions 
for the torsional modes ( Yg and Yu) and bimoment (Y j_q and 
Y 12 ) for a beam with both ends fully clamped (in sense of 
vanishing deplanation and twist). Having recognized this, 
one can write the relation for the generalized nodal force 
vector in the case of flexural torsion problems as 

1 1 , 

{Q.} = / Y. (z) m (z) dz - / Y,-(z) m (z) dz ... 2.20 

1 o 1 o n 

In case of concentrated forces, one can use Dirac delta 
functions to evaluate integrals in Equation 2.20. 
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The same relation has been obtained by Archer (1965) 
using virtual work argument and Betti's principle. 

2.9 MATRIX FORMULATION OF THE PROBLEM 

The general non- linear dynamic problem of elastic 
thin-walled structures presents considerable mathematical 
difficulties. However, a majority of problems of practical 
significance may be treated in a much simpler form assuming 
the steady state motion of the structure and constancy of 
parametric loads. This linearization essentially transforms 
the non-linear problem into an eigenvalue problem of much 
simpler form. 

In accord with the adopted model, it is further 
assumed that the loads are reduced to nodes. Although this 
assumption is not central to the derivation, it results in a 
much simpler form for all the governing matrices. 

The linearized differential equations of dynamic 
equilibrium defining the eigenvalue problem are (Figure 2.2): 

Transverse vibrations in x-z plane, 

El +P € + (M + a P) <p" + PA(i + a cp) = 0 ... 2.21(a) 

Transverse vibrations in y-z plane, 

El rf" + P T)' + (m - a p) 9 + P A(V) - a « ) = 

yy y -a. 


o 


... 2 . 21 (b) 
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Torsional vibrations about the longitudinal axis through the 
shear centre. 


.2 


EI na 9 + (r P + 2P v M^ r - 2p v M v - GK) 9 + (M v + a v P) l 


x y 


Y x 


(M y - a x P) n" + PA(a y 5 - a^Vj + r 2 9 ) = 0 


2.21(c) 


where terms in addition to those in Equations 2.4 are either 

due to the distortion of the beam or to the inertrci. External 

load is expressed through the resultant axial force P and 

resultant bending couples M and M . The co-ordinates of the 

"x y 

shear centre are a and a , and P is the mass per unit volume 

x y 

of the beam material, other parameters are 


2 ^xx I yy 2 2 

r = — LX- + a +■ a 

A x y 


P x = 2f~ /x(x2 + y 2 ) ^ - a x 


XX 


and (3, 


■X— f y(x 2 + y 2 ) dA - a 

yy 


... 2 . 22 (a) 


... 2.22(b) 

... 2.22(c) 


The differentiation with respect to time is indicated by dots. 

It may be noted that the left hand sides of Equations 
2.21 represent the resultant forces in direction of co-ordinate 
axes x and y and the resulting torque about the shear centre 
axis. Multiplying them by a set of virtual displacements 
6C, 6-r) and 69 and integrating over the whole domain one gets 
from Equations 2.21 
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1 

/ 

o 


{ [ EI xx 5 


+ [ei rj ' 

yy 


+ P£ 


+ P t) 


+ <M X 

+ (M y 


+ a p ) 

y 


a P ) 
x 


II 

cp 


9 


+ p k + a y 9 ) ] 6 K 
+ pA.(Vj - a 9 ) ] 6r) 

-K, 


+ f EI aa l f"” + (r3p + 2 Vy - 2 Py M x 


GK) 


9 + (M x + a y p ) 


II 

£ 


+ (My - a ^P ) T) + pA.(a 5 - + r 2 9 ) ] 69 } dz = O 

... 2.23 

Equation 2.23 simply states that the work of all -internal 
and external forces of an equilibrated structure vanishes for 
any arbitrary system of kinematically admissible virtual 
displacements , which is the well-known principle of virtual work. 

It may be recognized that Equation 2.23 may be inter- 
preted as the Galerkin formulation . 

Integration of Equation 2.23 by parts , with the 
interchange of derivation and variation leads to 


/ [El xx £ 6( £ ) + EI vv ti 6 ( r) ) +El fln 9 6 (9 ) + GK9 6 (9 )] dz 


yy 


+ / {p[ 56 ( 5 ') + T) ’ 6 (rj ' )} + 2{k + a p) <p' 6 c£ ) 


x y 


+ 2 (m - a x P) 9 6(rj ) + (tp + 2p x M - 2p M) 9 >6 (9 )} dz 


+ p / {A[ £6(0 + Y)6(r)) + f 2 g>6 (9)1 + * xx ‘i' 6 (?') 
o 

+ 1 yy Yi'd (t) ) + I nn $ 6 (9 ) + 2 A[ Sy *96 (9) - a x $6 <r) )] } dz 


= j o (o, 1 ) 


... 2.24 
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where J (0, 1), usually called conjunct or concomitant is 
the sum of all integrated terms. For natural boundary 
conditions 3 to be treated in sequel, conjunct vanishes and 
the problem is said to be self-adjoint. 

Introducing for displacements q(z, t), Equation 2.5, 
Equation 2.24 may be written in the form 

{6q ± } ([k_] - fgr ± j ]) Cqj) + {6q i >[m i j] } = 0 ... 2.25 

where [k. . ], [g. .] and [m. .] are called the stiffness , the 
stability and the mass matrix respectively. [g\^] denotes 
the form of [g . . ] , multiplied by a scalar. 

Finally, making use of Castigliano ' s first theorem , 
one obtains 

(Ojj] - t g ij^ {q j } + t m ijJ {q j } = {Q i } 2 * 26 

where (q^) is the vector of generalized nodal forces defined 
previously in Equation 2.18. 

It may be noted that the most frequently used way 
of derivation starting from stresses and strains is not a 
straightforward one in this case. This is due to the basic 
assumption about the deformation. 

2.10 STIFFNESS, STABILITY AND MASS MATRICES 

2.10.1 Stiffness Matrix 

The stiffness matrix [Kyi is defined by the first 
integral on the left hand side of Equation 2.24. It is 
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apparently a quasidiagonal matrix of twelth order, of the 


form 




i k' 


xx 


[k. .] = 

13 


! 0 


0 

_yy 


o 


o 

o 

kW 


2.27 


The first two submatrices representing stiffnesses 
in two principal planes are well-known from standard references 
such as Przemieniecki (1968), Ziehkiewicz (1971) etc. When 
expressed in the non-dimensional form, these are 


r, xxt EI xx 

[k J = — T" 


and [k yy ] = 


El 


III 


12 

-6 

=12 

-6 

12 

6 

-12 

6 


=6 

4 

6 

2 

6 

4 

-6 

2 


»12 

6 

12 

6 

-12 

-6 

12 

-6 


=6 

2 

6 

4 

6 

2 

-6 

4 


The elements of matrix are 

k?? = El / (y'.Y! + k 2 y!y') dz 
ij nn o 1 J 13 

i, j = 9, . . . , 12 

so that the final elemental submatrix [k^ reads 


2 . 28 (a) 


2.28(b) 


2.29(a) 



62 




EI Qfi 

DlT 


K^shic K^(l-=ch k) 


3 . 

< sh k 


K (l-ch K) 


k ( k ch K =sh < ) 


K (ch< -1) K(sh< - <) 
2 


3 

K sh< 


Symm, 


K “(ch l< -1) 

K ( K ch K “Sh K ) 


... 2.29 

where D and k are given by relations 2.15 and 2.16 respectively, 
2.10.2 Stability Matrix 

The elemental stability matrix [g^.] as defined by 
the second integral in Equation 2.24 is of the form 




g 


k * 


xx 

0 

Cpx 


0 


g- 


yy 


X( P 


; y 9 


g 


cpy 


g 


<p<p 


... 2.30 


Again/ referring to Przemieniecki (1969), submatrices 

XX yy 

[g ] and [g ], corresponding to flexure in the two principal 
planes, can be represented in the non-dimensional form as 


- 


P 

1 


6 

1 

6 

1 

5 

10 

5 

10 

1 

2 

1 

1 

10 

15 

10 

30 

6 

1 

6 

1 

~ 5 

10 

5 

10 

1 

1 

1 

2 

10 

30 

10 

15 


... 2.31(a) 


and 
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[s Yy ] 


p 

1 


6 

1 


6 

1 

5 

"io 


5 

10 

1 

2 


1 

1 

10 

15 


10 

*30 

6 

1 


6 

1 

“ 5 

10 

+ 

5 

10 

1 

1 


1 

2 

io 

30 


10 

15 

being 

considered 

as 

conrpres 

sivi 


... 2.31(b) 


submatrices [g^ X ] and [g^ y ] are transposes of [g X<? ] and [g y<J> ] 
respectively, using the second integral of Equation 2.24 one 
needs to evaluate submatrices [g X< ^] , and [ g99] only. 

The elements of these symmetric matrices are; 


x 

g . 

ij 


9 = 


Y9 = 
i j 


2 / (a P + m ) Y .Y • dz 
o y x i j 

(i =i; 4 , j = 9, ..., 12) 

i . . 

- 2 / (a^P - M y ) YjYj dz 

( i = 5 , ..., 3, j =9, ..., 12) 


2.32(a) 


2.33 


and g99 = / (r 2 p + 28 M - 26 M ) y y dz 

^ Y Y j- J 

(i/ j = 9, ..., 12). 

2.10.2.1 Submatrices [g Xc ^] and [ g Y 9] 


.. 2.34(a) 


Making use of Equation 2.32(a), 
to generate the following coefficients 


it is sufficient 
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g 2 ,lO = g 4 , 12 = 7T1 r- 1 + fo (« ch k « sh k)] 

1 K 

and 94,10 = g 2 ,i 2 = rr^ r- 1 - s < K - sh 

J- K 

... 2.32(b) 

where the symbol ' th 1 refers to the hyperbolic tangent function. 

The submatrix fg X< ^] can be obtained by multiplying 

the coefficients in relations 2.32(b) by the term 2(a P + M ). 

y x 

The coefficients of the symmetric submatrix [g^] 
may be easily deduced from relations 2.32(b) by increasing the 
first index Igy four multiplying them by the term 

-2<a x P - M y ). 

2.10.2.2 Submatrix [g^] 

Using Equation 2.34(a), one needs to generate the 


following coefficients 
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g 


9,9 

g 9 ,11 

1 g ll,ll ~ 2. 



D . 

N 

o 

If 

" g lo,n 

- ~~ ! T~ y S [ ( 4 + 
D 1 


9 12 ,12 

— .... A. r ( r y. 

h - 4 

o 

h- 1 

o 

3 2 L(cn 
Kl D 


[(1 - ch <)(3sh k - k) + k sh 2 x ] 


2 2 

[ ( 4 + + ~~~ sh k ) ( 1 - ch k ) +• 29h 2 _ k ] 


^ — 2 t(ch K “ 1 ) ( sh k + k) + k sh k(k « 2sh k) 
2 

“ ( k - sh k ch k ) ] 


g iO,12 


[( sh K 

<1 D 


2 

3k )(1 - ch <) + ( k ch k - 3sh <)] 


... 2.34(b) 

and 9 11/12 = -g 9#12 = -s 9#10 • 

The submatrix [g^] is obtained by multiplying the 
coefficients given by relations 2.34(b) by the factor 
(Pr 2 + 2P x M y - 2p y M x ) . 


2.10.3 Mass Matrix 


The mass matrix defined by the third integral in 
Equation 2.24 is of the following form 

... 2.35 

r XX *1 

Referring to Przemieniecki (1968), submatrices [m J 
and [m^] , corresponding to transverse vibrations in the two 
principal planes, can be represented in the non-dimensional 




m 


m 


_xx 

0 

9x 


m 


m 


o 

yy 

cpy 


m 


xcp 


m 


m 


yep 


epep 


form as 
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r xx-i 
[m J 


p Al 
4 20 


and 


[m yy ] = 


PAl 
4 20 


156 

■=■22 

54 

13 

-22 

4 

“13 

“3 

54 

-13 

1 56 

22 

13 

-3 

22 

4 

156 

22 

54 

“13 

22 

4 

13 

“3 

54 

13 

156 

“22 

-13 

=3 

“22 

4 


.. 2.36(a) 


... 2.36(b) 


As the submatrices [m^ X ] and fm are transposes 
of [m X( ^] and [m^] respectively, using the third integral of 
Equation 2.24, one needs to evaluate submatrices [m X< ^] , 

[m YC P] and [m^] only. The elements of these symmetric 
submatrices are 


X(j)' 

mj. 


2 pA a / y i Y 1 dz 


(i =1, ..., 4, j =9, ..., 12) 


... 2.37(a) 


mY? = - 2p A a / v.y.dz 

ij 1 3 


(i = 5, ..., 8, j = 9, ..., 12) 
_ 1 


and m?T = P A r s Y, Y, dz 


13 


i 3 


... 2.38 


.... 2 . 39 ( a) 


(i, j =9/ •••/ 12 ) 
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2.10.3.1 Submatrices [m X ^j and [m^^] 

Making use of Equation 2.37(a), it is sufficient to 
generate the following coefficients; 

m l,9 = m 3,ll = - k 3 + Kr 5 )shK + (24 + ^-)(1 - ch k)] 

m 2,9 =_ -m 2 / ll = " [ {6k + fo )shtc + (12 + ^ + ^ 2 ) (1 - ch tc)] 

m 3,9 = m l,ll = ~4 ^ ( " 12 ^ + lc3 + §5 K 5 )sh^ + 

k 4 

(“ 24 + 2~)(1 - ch k)] 

m 4/9 "” m 4 ,11 = 6lc + ^o )shK + (- 12 - k 2 + ~) 

(1 =» ch k) ] 
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m 4,lO “ m 2 / 12 


1DK' 


3 

K 

12 - 


[ (k * T~o)sh K 


(4 


+ _i_ ) oh < 

30 ; 


+ 4 + 


4 

20 J 


The symmetric submatrix [m X< P] can be obtained by 
multiplying the coefficients in relations 2.37(b) by the 
term 2P A a . 

The coefficients of the symmetric submatrix 
may be easily deduced from relations 2.37(b) by increasing the 
first index by four and multiplying them by ~2 P A a . 


2.10.3.2 Submatrix 

Using Equation 2.39(a), one needs to generate the 
following coefficients 

m Q Q =m - = —*r ~~ [(- 5sh k+2< - k ch k + ^sh <)(1 - ch < ) 

9,9 H/ll 1 d < 

+ (“• - 2 <) sh 2 < ] 


9,10 = 777 Ksh K - f H1 - chK > + 2K2sh2 K 

3 

- ~~ sh k(1 + 2 ch <) ] 

6 


m 9 , 11 °* 5 “ m 9,9 


m 


9,12 1 2 d 2 |< 2 


1 [-( 8 + — shic + ~) (1 - ch k) - 4sh 2 k 


.2 . 2 


sh k ( 2 + ch <) - k sh < ] 


+ 



* 6 9 


m 


10,10 = m 12,12 = “~ 2^3 [(3K + 3sh K)(1 - ch K) + b (7 + 2ch K) 


k sh k( 2 + -| ch k) + (6k + sh 2 K 1 


... 2.39(b) 


It may be noted that after checking the expression 
for the coefficient as given by Krajcinovic (1969), 

the term in last parantheses in the original expression was 
seen to be incorrect. The form appearing above is the amended 
one, the correction being indicated by the underscoring. 
Continuing 


m 


10,11 


1 Z n 

~ [(2 - §-)(! - ch K) + K(2sh< - K ch K)J - 


2 2 
1D< Z 


m 


9,10 


m 


1°, 12 I^D 2 * 3 


[(5k + | k 3 - 3sh lc )(l - ch K ) 


ch < 


+ K 2 sh k(“ + ch k) - k 3 - (2k + ~~) sh 2 < ] 


and 

i k 2 

[- (2 + -r-)(l ~ ch <) + «(« - 2sh k)] - m 9 12 
11,12 < ' ' 

... 2.39(b) 


The submatrix is obtained by multiplying the 

2 

coefficients given in relations 2.39(b) by the factor P Ar • 
It may be noted that, in the final account, all 
three matrices (stiffness, stability and mass) are symmetric, 
reflecting the fact that the problem is self-adjoint. All 
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derived matrices are computed in local co-ordinate system. 

In order to form the global 'matrix., necessary topological 
(transformation) matrices should toe used in the same way as 
for solid beam assemblages. 

2.11 SOLUTION OF THE PROBLEM 

Equation 2.26 is general enough to briefly discuss 
a host of distinct problems — such as static analysis, elastic 
instability, free steady state vibrations etc. It is empha- 
sised that the displacement method is used as the basis for 
the analysis of any of these problems and Equation 2.26 
enables one to evaluate the nodal displacements { qJ , being 
considered as unknowns in the problem. 



Having generated vector (F^ from the known 
external loading and matrix [k^j] assembled for the structure, 
from the known properties of the component thin -walled members, 
one can determine the unknown displacements, using Equation 
2 .- 41 , for a given set of boundary conditions. 
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This enables one to work out the internal forces 
at any section i.e. the generalized nodal forces which are 
expressible in terms of the nodal displacements/ now known. 
The stresses and deformations at any section can then be 
worked out using the internal force-stress and internal 
force-deformation relationships which take into account the 
prop erti.es of the section, in other words, this approach 
essentially means the stiffness method of analysis. 

2.11.2 Structural instability 

Assuming that all external (conservative) forces 
are characterized by a single parameter say > , one can write 
for the homogeneous case, from Equation 2.25 as 

([k^j] — ^ [^ij ]) {*1^} =0 ... 2.42 

Equation 2.42 represents a well-known eigenvalue problem 
wherein the nontrivial solution for { qJ - is obtained if 1 is 
an eigenvalue of the matrix [g. .] ^ [k. .]. The buckling load 
calculated from the smallest eigenvalue should be an 

upperboybd to the exact solution as a consequence of the 
applied variational principle. 

2.11.3 Vibration Analysis 

The general vibration problem (neglecting damping) 
of small oscillations about the equilibrium position is 
governed by Equation 2.25. For the steady state solution to 
the problem, the solution can be sought in the form 
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qj(t) 


qj exp 


(i ut) 


... 2.43 


where 1 • 1 is the natural frequency of the system and * i 1 the 
imaginery unit. Hence , Equation 2.25 can be written for the 
homogeneous case as 






Kj] + rhd V * 


ij 




= o 


... 2.44 


the curl denoting that the premultiplying scalar is not 
considered . 

The natural frequencies are now the eigenvalues of 
-1 - ^ 

the matrix [m . . ] ([k. .] - [g. .])• For the same reason as 

in the previous case/ the obtained natural frequencies are 
upper bounds to the exact solution. 


2.12 A LIMITING CASE 

Although the analysis of an assemblage of thin -walled 
members according to the procedure presented/ does not involve 
special computational difficulties, it is possible to simplify 

I 

it further in certain cases depending on the member geometry. 
The differential Equation 2.4(d) governing the problem of the 
torsion differs only by its second term on the left hand 
side from Equations 2.4(b) and (c), governing the flexure of 
the beam. The magnitude of this term is, furthermore, 
dependent on the parameter k (or in non-dimensional form k) , 
being the square root of the ratio of tors'ional rigidity 
(GK) and warping rigidity (EI M ). It is, therefore, of 



73 


interest to analyze the influence of this term on relations 
derived heretofore. The analysis of a limiting case < -*■ 0 
will not only give simple formulae but will provide the 
possibility of a qualitative analysis of the problem, in 
the sense that a very simple relation between the torsion 
and bending can be established. 

When the torsional rigidity can be neglected in 
comparison with the warping rigidity (GK << EI nfi , i.e. kl 0), 
the second term in Equation 2.4(d) indicating the effect of 
the st. ven ant's torsion can be discarded. Thus Equation 
2.4(d) is simplified to 

‘P = m fc + m fl ... 2.45 


The solution of the homogeneous equation 2.45 is 
obviously a cubic polynomial. Moreover the displacement 
modes yf have to be identical with corresponding modes for 
bending Using an adequate number of terms of series 

expansions for hyperbolic trigonometric functions 'sh' and 
'ch', the displacement function Y g (z) may be written as 

Y g (z) = D-j/D ‘ ... 2.46(a) 


where 


D, 


4 

K 

12 


2 2 
k (kz) 


+ 


k (kz)' 


... 2.46(b) 


and D 


2(1 -ch k) + k sh < 


... 2.15 
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as k approaches zero. 

As a result of this, for k - 0, the torsional moment 
can be evaluated as a transverse force and the bimoment as a 
bending moment,if the external torsional load is considered 
as a transverse load and the bimoment as distributed couples. 

in a similar manner, it may also be proved that 
when K approaches zero 

lim r*W] = 4^- - 4**] - S ay. 

1 I 

K ■* 0 XX 

Similarly, 


lim 

[S x? ] - 

lim 

r XX 

= [g J 

K+ o 


«-*- o 


lim 

[m x< f] = 

lim [ m^] 

11 

t — i 

3 <? 

o 


k-*0 



and 
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This means that all matrices differ only by a multi— 
P^ ca " t " or having the same form as for solid beams. Hence, this 
limiting case, besides being a useful and simple approximation 
(when applicable), justifies entirely the procedure presented 
as being nothing but an extension of established procedures 
for the solid beams. 

2.13 CONCLUDING REMARKS 

The presented matrix formulation for the static and 
dynamic analysis of the structures assembled from thin-walled 
members is by itself quite general and consistent. A number 
of related linearized problems such as lateral buckling, 
second order theory, parametric resonance etc. may also be 
treated once the stiffness, stability and mass matrices for Uhe 
structure are assembled using the relations established 
earlier . 

However, one important point may be emphasized in 
this regard. In case of linear static problems, the solution 
presented herein is exact since the displacement modes are 
the exact solutions of the governing homogeneous Equations 
2.4(b) and (c). For the dynamic and stability problems, the 
story is somewhat different. The choice of displacement •* 
modes is a matter of assumptions and thus, subject to 
discussion. In the present case it is assumed that the 
static deformation modes are similar in shape to buckling 
and dynamic modes. Quite often this assumption was 
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successfully employed in conjunction with variational methods, 
long before the advent of the finite element method. However, 
as a consequence, the buckling load and the natural vibrational 
frequency obtained using this formulation are expected to be 
upperbounds to the exact solutions. 


- o - 



77 


3. OPTIMAL DESIGN OF THIN -WALLED SECTIONS 
SUBJECT TO STATIC CONSTRAINTS 

3.1 BRIEF LITERATURE SURVEY 

3.1.1 Analysis of Thin -Walled sections 

Before the publication of the classical treatise on 
thin-walled elastic beams by Vlasov (1959), attempts were made 
by Patterson (1952) and Nylander (1956) to study combined 
flexural torsion and lateral buckling of I -shaped beams of 
mono-symmetrical cross-section. 

El Darwish (1965) provided charts and empirical 
formulae for the accurate calculation of torsion constant for 
different structural shapes of open cross-section. Nethercot 
(1974) proposed a relation, based on experimental results, to 
determine the torsional rigidity of 1-sections, and concluded 
that the presence of residual stresses increases it signifi- 
cantly. Marshall (1971) presented a basis for determining 
displacements and stresses arising from torsion of structural 
rectangular hollow sections, with particular attention to 
stress concentrations at the re-entrant comers. 

Home (1971) studied the behaviour of thin -walled 
prismatic members of doubly symmetric open cross-section in 
which the second moment of area about the major principal 
axis is very much greater than that about the minor principal 
axis. The governing equations for the member were derived 
and the solution obtained. 
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Conway (1972) proposed a very simple formula for 
calculating the stress at the re-entrant comer of a twisted 
structural angle or channel. 

Chang (1975) derived a method of computing the 
warping shear stresses due to restraints in a thin-walled beam 
of open cross-section, using finite elements. 

Holden (1972) presented the numerical solution of 
three different types of problems of finite deflection of 
uniform thin beams using the Bemoulli-Euler relation. 

3.1.2 Optimization of Thin-Walled Sections 

Not much work appears to have been done on the 
optimal designs of thin -walled structures. 

Haaijer (1961) obtained optimum proportions that 
result in the maximum strength or stiffness for a given cross- 
sectional area in I -shaped homogeneous beams of high strength 
steel. Schilling (1974) extended this work on hybrid beams 
and gave equations and curves showing the variation of bending 
strength and stiffness, and other related properties with 
various geometric parameters for both hybrid and homogeneous 
be ams . 

Barnett (1961) developed design procedures for 
minimum weight, statically determinate beams using a deflection 
criterion. Optimum flange width and web thickness variations 
were described for I-beams designed with respect to bending 


and shear deflection 
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Dupuis (1971) pointed out that the optimal design of 
beams subject to deflection constraint alone might not be a 
well posed problem. As such, he' introduced an upper bound on 
the bending stress as an additional constraint and further 
showed that these two behavioural constraints were complementary. 

Felton (1971) and Nelson (1972) considered basic 
aspects of the design of indeterminate structures composed of 
thin^walled beam elements. The minimum weight design of such 
elastic framed structures was obtained, subject to constraints 
against two types of local instability as well against yield. 

The solutions were examined numerically by using a non-linear 
programming technique. 

Bronowicki (1975) pointed out that the optimal 
design formulation of the type above, although infeasible in 
practice, provided a lower bound on obtainable weight. After 
reformulation and utilizing a mathematical programming search 
technique, he obtained optimal designs of potentially more 
practical utility in the case of continuous thin -walled beams. 
The set of behavioural constraints was expanded to include 
effects of shear stress on web buckling, a factor which was 
not considered in the previous studies. 

However all these papers have not considered 
the effect of warping of the section and the stresses 
associated with it , although this is a chharacteristic feature 
of a thin-walled beam, in obtaining the optimal solutions for 
I -shaped sections, the transverse loading is considered to be 
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in the web-plane, in practice, this being seldom possible, 
the torsion created is likely to produce the warping stresses 
of such magnitude worthy to be considered in the design. 

Hence, in the present thesis, these foregoing stresses 
have been regarded in various behavioural constraints consi- 
dered for the optimal design of I -sections. In addition to the 
usual in-plane loading, the one having lateral eccentricity 
with respect to the web plane is also included to consider 
torsional effects. 

Consequently, the discrete element analysis employed 
considers four degrees of freedom at each node of the beam 
element, as against the usual two. The additional two degrees 
are the longitudinat twist and the loarping of the section, 
the internal forces corresponding to these being the warping 
torsion and the bimoment , which together represent the 
effect of flexural torsion produced in any thin-walled beam 
member,. The optimal solutions obtained incorporating such 
provisions in analysis and design are therefore expected to 
be more rational,. 

3.2 DESIGN VARIABLES 

As it is proposed to obtain the minimum weight 
design by keeping the section constant throughout the member 
span, the problem indirectly reduces to optimization of the • 
cross-sectional dimensions. Referring to Figure 3.1(a), the 
given section can be looked upon as an assemblage of three 
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thin plates: the top flange plate of size (bf^, tf^), the 
bottom flange plate of size (bf^/ tf^) and these two being 
connected by the web plate of size (hw, tv/). These 
variables representing cross-sectional dimensions can be 
viewed as to constitute the design vector 


X = { x ± } 




v tw J 


... 3.1 


The plate elements being thin -walled , the section 
can be 'idealized and considered for analysis in the form 
shown in Figure 3.1(b). It is aimed at maintaining the 
sectional symmetry with respect to the minor principal axis, 
so that the centroid G and the shear centre S will always lie 
on the x-axis. With these assumptions, various properties of 
the section, with G as the origin of co-ordinates x and y, 
can be written as [Figure 3.1(b)]. 


Sectional area, A = bf^.tf^ + bf 2 . tf 3 + hw.tw ... 3.2(a) 


Position of the centroid G, 


x 


hw.tw. + bf^.tf^.hw 
A 


3.2(b) 
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3.3 DISCRETE ELEMENT ANALYSIS 

Consider a thin-walled beam of I -section carrying 
any transverse loading and subjected to any boundary conditions 
[Fig. 3.2(a)], The loading is considered to lie in the web- 
plane [Fig. 3.2(b)] , although in some cases it can as well be 
eccentric with respect to the web plane [Fig. 3.2(c)] . 

Adopting a two-element discretization of the structure and 
noting that the degree of freedom at each node is four, the 
generalized unknown nodal displacements and the corresponding 
nodal forces are as shown in Figure 3.3. 

The displacement vector {q^} for the structure can 
be written down in the form of subvectors {q A ), {q B } and {q C } 
such that 

V 


where 
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. * . 3.5(b) 


... 3.5(c) 


The corresponding nodal force vector {Q^} for the 

A 

structure is similarly a combination of subvectors {Q }/ 
{Q B } and {Q C } such that 


{Q.} 


where (Q } = 


f A 
Q 


* 

< 


V 


(?) 


Q- 


W 



... 3.6 


... 3.7(a) 
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and 


{Q C > = 


{q b > 


f \ 

i 0 c 


Q- 




/ Q < 


Q 8 } 


\ 


Q 


10 


'11 


V = 


V 


It may be noted that the displacements {q } and 



3.7(b) 


3.7(c) 


B 


{q } are governed by the prescribed boundary conditions . 

3.3.1 Equivalent Nodal Loads 

The governing equilibrium equation for calculating 
the nodal displacements in a structure is 


{F i } 


= [k. .] 

ij 


{q. } 
tL 


... 2.40 


where is the vector of external equivalent nodal loads 

corresponding to the displacements / given by Equation 3.4. 

In general, the actual loads on a structure may not 


be necessarily acting at the nodes. Hence these are divided 
into two types: loads acting at the nodes and loads acting 
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on the elements of the discretized structure. Loads of the 

latter type are replaced by their statically equivalent loads 

obtained as the fixed, end. vecctions and written as the 

vector {r}. These reactions when reversed and added to the 

actual nodal loads, generate the vector (f^}. 

For the types of loads considered in the problem, 

one can readily make use of the relations given in Gere (1969) 

to obtain the fixed end reactions, provided the loads lie in 

the web plane, if they are eccentric [Figure 3.2(c)], these 

produce a varying torsion along the length of the beam. In 

order to determine the fixed end reactions namely, the 

twisting moment and the bimoment produced in such a case, 

attention is drawn to Article 2.12, with particular reference 

to the simplification achieved in evaluating these, when the 
2 

ratio k (warping rigidity/torsional rigidity) is very small. 

2 

In the present problem the value of k is seen to be in the 
=5 -5 

range of 2 x 10 to 4 x 10 , small enough to be neglected. 

Hence, the fixed end twisting moment and the bimoment can be 
evaluated as the transverse force and the bending moment 
respectively, considering the external torsional load as the 
transverse load and the bimoment as distributed couples. 

Thus, obtaining the unknown nodal displacements 
from the relation 

up = [^-r 1 tpp 


2.41 
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the vector of the internal forces or generalized nodal forces 
can be determined from 

{Q ± } = {R> + [k...] {q ± } ... 3.8 

Again, if the effect of the St.- Venant torsion (T ) 

s 

is neglected. Equation 2.1 yields T fi - T. 

3.3.2 Internal Nodal Stresses 

Having obtained the four internal forces, v^, M^/l, 

2 

T^/1 and /I at every node, one can obtain the internal 
stresses at a node. 

As mentioned in Article 2.3, a loading state at a 

section in a thin-walled beam can be considered as a super- 

position of two independent states — the flexure and the 

flexural torsion. As such, the total stress at any point can 

be summed up as due to the internal forces V , M (due to 

x y 

flexure) and Tq,M fi (due to flexural torsion) [Figure 3.4(a)]. 

The contribution to the normal stress and the 
shear stress at any point due to the flexural torsion is 
governed by the following two theorems of Vlasov; 

(a) The stress in a longitudinal fibre of a thin- 
walled beam due to a bimoment is equal to a product of this 
bimoment multiplied by the principal sectorial co-ordinate 
of this fibre divided by the principal sectorial moment of 
inertia of . its cross-section. 
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(b) A shear stress in a fibre of a thin-walled beam 
caused by a flexural twist is equal to the product of the 
flexural twist multiplied by the sectorial statical moment of 
this point divided by the wall thickness (at this point) and 
the principal sectorial moment of inertia. 

Hence, in the present case of an I -section, 
combining the effects of the flexure and the flexural torsion/ 
the expressions for the maximum normal stress and the maximum 
shear stress at a node can be written as follows [Figure 
3.4(b)] . 

3. 3. 2.1 Normal stress 


The maximum normal stress is expected to occur 
either at point (1) or (2) in the flange s in which case its 
values are 


bf 


( 1 ) _ y 

c - — j 


M„. (hw - b,) M n . 2 


x 


. ~ (hw - b + a ) 


x x 


XX 


“fifi 


... 3.9 


bf , 


( 2 ) 


M _.b . . --rr- (b„ - a v ) 

+ 


•y~~x . * 2 '“x x' 


I 


xx 


•fin 


3. 3. 2. 2 Shear stress 

The maximum shear stress is expected to occur 
either at point (3), ( 4 ) or ( 5 ) in the web and in the 
respective case, its value is given by 
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.( 1 ) 


( 2 ) 


.(3) 


V x . [bf ^ .tf - (hw 


tw 


bf, 

Vl + T oJ-§ 1 


tf ^ (hw 


b + a )] 
x x J 


XX 


tw • I ttQ, 

V x* r bf i* tf i(bw - b x ) + (hw - b^) 2 ] 


x 

tw . I 


XX 


bf: 


V bf 2- t£ 2- b x> . T n-f-s" • tf 2 (b x - a x >] 


tw . I 


+ 


XX 


tw . I 


£2fl 


3.10 


It may be emphasised that the second part in 
Equations 3.9 and 3.10 is the stress ■increase in comparison 
to the conventional solid beam theory s due to the effect of 

warping of the section which is' a characteristic feature of a 
thin -walled member. 



Thus, at 

any node i, the 

maximum normal and 

shearing 

stresses are 





max 

a i 

r d) 

= max [ 0 , 

(2) 
tf 3 

... 3.11(a) 


T max 

l 

r (1) 

= max [ % , 

(2) (3) 1 

T / T J 

... 3.11(b) 

and 

in the whole 

domain of the 

beam, the critical 

values of 

the 

normal and the shearing stresses, to be considered for 

the 

optimal design having a unif orm section are 



^cr 

_ r max-, 
max [ow ] 

(i = 1,...,3) 

... 3.12(a) 


T cr = 

r max-, 
max [t^ J 

* 

(i 

... 3.12(b) 


for the two-element solution. 
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As already discussed in this article, once the nodal 
displacements {q ± } are obtained using Equation 2.41, the 
critical deformations that can be considered in such an optimal 
design are 


^cr 

= max 

11*!' 


q 9 

< ^ ' 1 hr 

= max 

[ l<v 

q 6' 

o 
\ — 1 

D 1 

<9D cr 

= max 

[ |q 3 / 

q 7' 

q ll 

(9 ' l 2 ) 

T cr 

= max 

Uq 4 / 

q 8' 

q 12 


... 3.13 


3.4 CONSTRAINTS IN THE OPTIMIZATION PROBLEM 

Out of the sixteen constraints prescribed in the 
problem, twelve are the side constraints, the remaining four 
being the behaviour constraints. Recalling that the design 
vector is 



the constraint relations can be obtained as follows. 
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3.4.1 Side Constraints 

3. 4. 1.1 Lower and upper bounds on the design variables 

Using subscripts 'min' and 'max' to mean the minimum 
and maximum values respectively of the design variables bf^/ 
tf^/ ... etc., the constraint relations in the normalized 
£ orm are 


g x (x) 

= 

or . 

1mm ^ 

x i 

... 3.14(a) 

g 2 (x) 

= 

t: ^lmin ^ 

x 2 

.... 3.14(b) 

g 3 ( X) 

= 

b£ 2min ^ 

X 3 

... 3.14(c) 

g 4 (X) 

= 

tf „ . 

2mm ^ 

x 4 

... 3.14(d) 

g 5 (x) 

= 

hw . 

mm 4 

x 5 

... 3.14(e) 

g 6 (x) 

= 

tw . 

mm . 

x 6 

... 3.14(f) 

g 7 (x) 


*1 - 1 

... 3.14(g) 


b£ lmax 

g 8 (x) 

= 

X 2 

1 - 1 
tf lmax 

... 3.14(h) 

g 9 (x) 

= 

X 3 

J - 1 
b£ 2m ax 

... 3 . 14 ( i ) 
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g 10 (x> 

X 4 

tf _ ^ 

2m ax 

... 3.14(3) 


g lx (x) 

X 5 

= - — ±_ „ i 

hw 

max 

... 3.14(k) 


g 12 (x) 

= _ 6 • - 1 
tw 
max 

... 3.14(1) 

3.4.2 Behaviour 

Constraints 


3. 4. 2.1 

Maximum 

normal and shear stresses in the beam 


Denoting the permissible normal and shear stresses 

as a 

per 

and t 

per 

respective ly , 



gi 3 ( x) 

a 

- 0^ 1 
per 

... 3.14 (m) 


gi4 ( x) 

T cr 1 

T 

per 

... 3 .14 (n ) 

3. 4. 2. 2 

Maximum 

deflection in the 

beam 


If 5p er -'- s the permissible def lection in the beam 


9 15<X) 

^cr _ ^ 

^per 

3.14 ( p) 

3. 4. 2. 3 

Maximum 

slenderness ratio 

of the compression flange 


If SR 1 

is the slenderness 

ratio of the compression 


flange with respect to the appropriate radius of gyration , as 
per the Clause 21.3 of IS : 800 (1962), this value should not 


exceed 300. So 
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SR 

g l6 ( * } = 305 ~ 1 *** 3 ' 14(q) 

3.5 DATA FOR THE ANALYSIS 

3.5.1 Number of Elements 

As has been established by Krajcinovic (1969), a 
two-element solution gives quite satisfactory results for the 
static analysis problems. Hence in the present case of opti- 
mization involving static constraints, the analysis using a 
two-element idealization [Figure 3.3(a)] is adopted. 

3.5.2 Boundary conditions in the Problem 

The support conditions chosen are: fixed, simple 
and one end fixed with the other end free (cantilever). 
Referring to the generalized nodal displacements shown in 
Figure 3.3(b), it is understood that a fixed boundary will 
restrain all the four displacements at the support node, a 
free boundary will allow these to occur while a simple 
boundary restrains only the latera-1 displacement and the 
longitudinal twist, allowing the transverse rotation and the 
warping to occur. Depending upon whether a particular 
displacement at a support is restrained or not, it is 
assigned a code number unity or zero. The generalized nodal 
displacements in a two-element discretized model have been 
numbered [Figure 3.3(b)] and using the coding system mentioned 

of boundary condition can be generated . 


above, any type 
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3.5.3 Elastic constants 

The beam material being structural steel/ in accor- 
dance with the clause F-l.l (Appendix f), IS s 800 (1962)/ the 
values of the elastic modulus 'E* and the modulus of rigidity 
'G 1 have been taken as 2.047 x 10^ kg/cm^ and 0.4 E respect- 
ively. 

3.5.4 Span of the Beam 

The range of span L chosen is from 3.0 to 6.0 m for 
the fixed and the simple supports/ and from 1.0 to 2.0 m for 
the cantilever. 

3.5.5 Types and Range of Loading 

Both the nodal and the non--nodal (applied on the 
elements) loads of the following types have been considered 
in the respective ranges. 

(a) Nodal loads equivalent to the uniformly distributed 
load of intensity range from 1000 to 3000 kg/m. 

(b) Uniformly distributed load of intensity range from 500 
to 3000 kg/m. 

(c) Linearly varying load of intensity range varying from 
0-2500 kg/m at one end to 500-3000 kg/m at the other. 

(d) Uniformly distributed load of intensity range varying 
from 1000 to 3000 kg/m, having a uniform eccentricity 
• e 1 with respect to the web plane. The range of this 
eccentricity is from 0.1 to 0.3 times the top flange 


width . 
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3.5.6 Effective Length of the Compression Flange 

In the present problem, it is assumed that no 
lateral restraint is provided for the compression flange but 
the beam is -partially restrained against torsion. Referring 
to the clauses 21.2.1 and 21.2.3 of IS ; 800 (1962), the 
coefficient for calculating the effective length of the comp- 
ression flange is taken as 0.85. The slenderness ratio 1 SR^ 1 
of the compression flange thus calculated with respect to its 
appropriate axis is used in Equation 3.14 (q). 

3.5.7 Permissible Bending and shearing Stresses 

Using the clause 10.2.2.1 and Table VII of IS : 800 

(1962), the critical stress 'C 1 in the compression element 

s 

is calculated. The allowable stress in bending 'P^ ' as 

2 

given by Table III of IS : 800 (1962) is 1650 kg/cm . The 

value of the permissible bending stress c to be used in 

Equation 3.14 (m) is the lesser of C g and P^ c above. 

Referring to the clause 10.3.1 and Table VIII of 

IS : 800 (1962), the value of the permissible shearing stress 

2 

x in Equation 3.14(n) is taken as 1100 kg/cm . 
per 

3.5.8 Permissible Deflection in the Beam 

Referring to the clause 21.4.1 of IS j 800 (1962), 
the limiting value of the deflection * per to be used in 
Equation 3.14(p) is of the span, based on the strength 

or efficiency of the structure or safety against damage to 


finishing. 
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3.5.9 Lower Bounds on the Design variables 

Referring to the range of available rolled steel 
I -sections given in Table I of SP-6(1) (1964), the values of 

bf lmin' bf 2min ' hW min and tw m i n to be used in Equations 
3.14(a), 3.14(c), 3.14(e) and 3.14(f) are prescribed as 5.0, 
5.0, 7.5 and 0.3 cm respectively. 

According to the Clause 17.4 of IS : 800 (1962), the 
controlling thickness for rolled steel being the mean flange 
thickness, the Clause 17.2 of IS j 800 (1962) fixes the values 

of tf l m in and tf 2min t0 be used in Equations 3.14(b) and 
3.14(d), as 0.6 cm. 


3.5.10 Upper Bounds on the Design variables 

Referring to the range of available rolled steel 

I-sections given in Table I of SP-6(1) (1964), the values of 

hf 1 , tf. , bf , tf„ , hw and tw are taken as 
lmax lmax 2max 2max max max 

25.0, 2.5, 25.0, 2.5, 60.0 and 1.25 cm respectively, in 

addition, following conditions also govern the values of 

bf. , bf n and hw 
lmax 2max max 


3.5.10,1 Maximum outstand of flange 

As per the Clause 21.5.1.3 and Table XVIII of IS ; 
800 (1962), the maximum outstand of the compression and the 
tension flange is sixteen and twenty times the respective 
flange thickness. This indirectly put an upper bound on the 
values of bf^ and bf 2max . The values of bf lmax and 
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k^2max use< 3 in Equations 3.14(g) and 3.14(i) are the 

lesser of those given by Articles 3.5.10 ana 3.5.10.1 respect- 
ively. 

3.5.10.2 Maximum depth of the web plate 

Referring to the clause 21.6.1, the minimum thickness 

of an unstiffened web plate should not be less than hw/35 or 

in other words for a given value of tw, hw i 85 tw. This 

indirectly puts an upper bound on hw. The value of hw to 

max 

be used in Equation 3.14(k) is the lesser of that given by 
Articles 3.5.10 and 3.5.10.2 respectively. 

3.6 OBJECTIVE FUNCTION 

Weight being the objective function, for a thin-walled 
beam member of uniform cross-section, the objective function 
reduces to the expression for the cross-sectional area of the 
member. Hence 

f(X) = X-jX 2 + x 3 x 4 + x 5 x & ... 3.15 

3.7 OPTIMIZATION PROBLEM 

Thus, the present optimization problem is of 
constrained type, of the form; 


to find 
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so as to minimize the function 

f (X) = X-.X- + ... 4- x r x, ... 3.15 

1 z b 6 

and subject to the inequality constraints 

9j (X) .< 0 , j = 1, 2, . .., 16 ... 3.14 

The solution is attempted by the interior penalty 
function method, that is widely used for the constrained 
minimization problems. 

3.8 INTERIOR PENALTY FUNCTION METHOD 

In this method, a new function ' 0' , called the 
penalty function, is constructed by augmenting a penalty term 
to the objective function. The penalty term, involving a 
parameter r called the penalty parameter 3 is chosen such 
that its value will be small at points away from the const- 
raint boundaries and will tend to infinity as the constraint 
boundaries are approached. Hence the value of the 0~f unction 
also blows up as the constraint boundaries are approached. 
Thus, once the constrained minimization of 0(x, r) is 
started from any feasible point x 1# the subsequent points 
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generated will always lie within the feasible domain , 
since the constraint boundaries act as barriers during the 
minimization process. 

The 0-f unction used is that defined originally by 
Carroll (1961) 


0(X, r) = f(x) 


? -t- 

j=i g j (x) 


.. 3.16 


It can be seen that the unconstrained minimum of the 
penalty function approaches the constrained, minimum of the 
objective function, as r -* 0. The iterative process to 
obtain this minimum is as follows. 

(i) Start with an initial feasible point X^/ satisfying 

all the constraints with strict i?iequality sign, that 
is, g . (X^ ) < 0 for j = 1, 2, . m. set the counter 
k = 1 and choose the initial value of the penalty 
parameter r^ > 0. 

(ii) Minimize {zf(x, r^) using any of the unconstrained 

minimization methods and obtain the solution X k+ -j_ • 

(iii) Test whether X k+1 is the optimum solution of the 

original problem. If X k+1 is found to be optimum, 
terminate the process. Otherwise go to the next step, 

(iv) Find the value of the next penalty parameter 


'k+1 


= c r. 


where c 1 
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(v) Set the counter k = k + 1, take the new starting point 
and go to step (ii). 

Steps (i) through (v) are collectively known as the 
sequential unconstrained minimization tec unique (SUMT). 

In the present case/ the value of the initial 
penalty parameter r^ is taken as 

f (x x ) 

r i " — i : 

[- i 

j=l g j (x 1 ) 

and c = 0.1 

these values being found satisfactory. 

3.9 UNCONSTRAINED MINIMIZATION 

For obtaining the unconstrained minimum of the 
penalty function, the variable metric method by Davidon- 
p letcher^Powell (19.59, 1963) is used. This method is the 
best general purpose and currently available unconstrained 
optimization technique making use of the function derivatives. 
It is very powerful, stable and continues to progress towards 
the minimum even while minimizing highly distorted and 
eccentric functions. It is basically a descent method that 
makes use of the gradient vector in finding, at & point j 
a better search direction given by the vector 
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S . 


X 



{s. } = 

1 



V 

( 



3 17 

• • ♦ ' 


The iterative procedure of this method can be stated 
by the following algorithm: 

(i) Start with an initial point a nd an (n x n) positive 
definite symmetric matrix [H^] called the Hessian 
matrix. Usually [h.] is taken as the identity matrix 
[i ] . Set iteration number as i = 1 . 

(ii) compute the gradient of the function vf^, at hhe point 
X^ and obtain 

Si = - h 

•k — . 

(iii) Find the optimal step length A^ in the direction S^ 
and set 


X 


i+1 


X. 

1 


* _ 
+ 1 , s . 

x 1 


(iv) Test the new point X^ + ^ for optimality. If X^ + ^ ^- s 
optimal, terminate the iterative process. Otherwise 
go to the next step. 

(v) Update the [h] matrix as 


[H 1+ il = [ H ± ] + + l>iJ 
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where [m^] = 


* _ -u-m 

X , S S 
. i i 

-T — 

s i Q i 


nT 


[n ± ] 


(_r_H.] Q.)( [H.] Q. ) 

? [Hi] Q. 


and q ± = v f(x ±+1 ) - 7f(x.) 


Vf 


i+1 



(vi) Set the new iteration number i = i + 1 and go to 
step (ii). 

As the method makes use of the gradient of the 

function, it is preferable to use the more accurate oubic 

★ 

interpolation method for obtaining -h in step (iii) above. 


3 . 10 COMPUTER PROGRAMME 

The computer programme prepared for obtaining the 
results consists of two programmes: 

(i) STTANA - which deals with the static analysis part, 
and (ii) 0PT2 - which deals with the optimization part, 
based on the results supplied by STTANA. 

The functions performed by various subroutines of the 
programme STTANA are as follows. 

(i) CSPROP ~ calculates the cross-sectional properties 
of a thin -walled I -section. 

(ii) ELSTIF — Generates the stiffness matrix for a thin— 


walled beam element 
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(iii) ASSBLY “■= Assembles the stiffness matrix fk. .] for 

L ij 

the structure from the individual element 


matrices, it further takes out the submatrix 


corresponding to the free nodal displacements 
in the structure . 

( iv ) MATINV -- Performs the inversion of a matrix. 

(v) STAANA — - For a thin-walled beam of I -section of known 

sectional dimensions/ boundary conditions 
and loading, it calculates the maximum values 
of the normal* stress, the shear stress and 
displacement in the structure, after perfor- 
ming the discrete element analysis of the 


structure . 

The optimal solution is achieved by the programme 
0PT2 using following subroutines; 

(i) MAIN ~ Reads the data for the problem, performs 

various steps in the optimization process and 
finally returns the optimal sectional dimen- 
sions and the corresponding values of the 
objective function and the constraints. 

(ii) SUMT2 — Obtains the constrained minimum of a given 

objective function using the sequential 
unconstrained minimization technique (SUMT). 

(iii) PELTY2 -- Generates the objective function and the 

penalty function for a design point X. 
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(iv) TABLEl Gives the critical stress in the compression 

element as per Table VII of IS ; 800 (1962). 

(v) DFP2 — Obtains the unconstrained minimum of a 

function by the D av id on *»F le t che r =P owe 1 1 method, 

(vi) UNIDR 2 — Minimizes a function in one direction using 

the method of cubic interpolation . 

(vii) GRADF2 Gives the gradient vector and its modulus at 

a given point X in the function, by using the 
finite difference technique. 

3.10.1 Characteristics of the Programme 

The programmes STTANA and 0PT2 are quite general and 
possess the following features; 

(a) The programme STTANA can yield the two-, four- and 
eight -element analysis of thin-walled beams having an 

I -section , with different spans and boundary conditions, 
and subjected to a wide variety of loads as detailed in 

Article 3.5. 

(b) By numbering the nodes and generalized displacements at 
each node, and assigning the respective code number to 
each displacement, any set of boundary conditions can 

be generated. Thus, in addition to applying the 

. . , . „„ ,>, 0 ,,«,7 nodes, the intermediate 

boundary conditions at the noaes, 

nodes can also be subjected to such conditions so that 
the continuous beams can also be handled by the programme 
STTANA, using this versatile technique. 
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(c) In each subroutine of the programme 0PT2, care is taken 
to see that during the optimization process the design 
point does not enter the infeasible domain. Although 
theoretically impossible, such type of situation is not 
uncommon due to the numerical errors associated with a 
digital computer. 

(d) At every important stage in the programme 0PT2, 
precaution is taken to cater for a bound design point. 

In the vicinity of such a point, the constraints natu- 
rally ‘blow up' and to get rid of the consequent numerical 
instability, necessary logical provision is made in the 
programme . 

Thus, the programme has been made numerically stable 
to counter all possible difficult situations created during 
its execution. 

The results for all the cases above have been 
discussed and the conclusions drawn in Chapter 6. 


= o “ 
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4. OPTIMAL DESIGN OF THIN -WALLED SECTIONS SUBJECT 

TO STABILITY OR DYNAMIC CONSTRAINTS 

4.1 BRIEF LITERATURE SURVEY 

4.1.1 Analysis of Thin-Walled Sections 

Gallagher (1963) suggested a matrix displacement 
formulation of the beam-column problem and presented a general 
procedure including instability effects in element force- 
displacement relationships applicable to both beams and plates. 
Hartz (1965) presented rapid formulation of stability problems 
for column and frame type structures and of the beam-column 
type problem for determining the effects of axial loads on 
laterally loaded structures. 

Gallagher (1967) extended the concepts of the matrix 
displacement approach to discrete element structural analysis 
to predict general instability. The influence of element 
membrane forces on element effective flexure stiffnesses was 
accounted for and critical applied load intensities were 
predicted for plate, arch and spherical cap structures. 

Dawe (1969) applied the method to the buckling analysis of 
rectangular plates under arbitrary member loading. 

Chajes (1965) presented a simple method of 
accounting for torsional -flexural buckling in case of thin- 
walled open sections. Based on the use of an interaction 
type of equation, he suggested a relatively uncomplicated 
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and straightforward procedure for determining the torsional- 
flexural buckling load in case of singly-symmetric sections/ 
a grouping that includes most of the commonly used shapes. 
Kennedy (1972)/ with the help of experimental results on steel 
angle and tee struts/ advocated the use of a rational buckling 
analysis in stipulating the maximum width-thickness ratios. 

Barsoum (1970) formulated the stiffness matrices for 
the torsional and lateral stability analysis of structures 
composed of flexural members by the matrix displacement method. 

Kawai (1971) performed finite element analysis of 
thin-walled sections based on modern engineering theory of 
beams. Rajasekharan (1971) did the same for thin-walled 
members of open section. Nethercot (1971) demonstrated the 
accuracy of the finite element method in dealing with problems 
of buckling of columns and beams. 

Wittrick (1968) and Williams (1969) developed the 
stiffness matrices and the subsequent computational procedure 
for the buckling or vibration analysis of thin, flat-walled 
structures. It was shown that buckling and vibration 
phenomena for any structure of this type are closely inter- 
related . 

Carnegie (1969) presented a method of matrix 
displacement analysis, with strong convergence characteristics 
for the determination of eigenvalues and eigenvectors in 
vibration problems of continuous systems. 
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Popelar (1969) investigated the stability of the 
free transverse vibrations of a simply supported, thin-walled 
elastic beam with double axeS 1 of symmetry, for arbitrary 
initial velocity distributions. 

Mei (1970) derived the stiffness and consistent 
mass matrices for a thin -walled beam element of open cross*-’ 
section with non-collinear shear centre and centroid. 

4.1.2 optimization of Thin-Walled Sections 

zarghamee (1970) developed a mathematical programming 
technique for the minimum weight design of such' structures 
wherein the general instability may be the controlling 
phenomenon in the design. The method was claimed to be appli- 
cable to all structures for which a finite element model can 
be built. 

Khot (1976) presented an optimization method based c*: 
on optimality criterion for minimum weight of structures with 
stability constraints. The method was developed in the 
displacement method of finite element analysis aid was pro- 
grammed for trusses and frames. 

Taylor (1968) introduced a method for the optimum 

design of a vibrating bar with specified minimum cross- 
section . 

Kamat (1973)' presented a finite element displacement 
the first mode natural frequency of 


formulation to maximize 



113 


a vibrating beam of specified volume with elastically - rest- 
rained ends and resting on a continuous elastic foundation. 

Sheu (1968) developed a method for the elastic 
minimum weight design of a one ‘’dimensional structure that is 
to have a prescribed fundamental frequency. Similar methods 
were also suggested by Turner (1967) and McCart (1969). 

Fried (1969) proposed gradient methods for finite 
element eigen problems. 

However, all these papers above have not considered 
the effect of warping of the section, although this is a 
characteristic feature of a thin-walled beam, in the present 
analysis, this aspect is considered by giving due regard to 
the degree of freedom at a section . 

4.2 STABILITY ANALYSIS 

Consider a prismatic member of length L, having a 
thin -walled, I -shaped cross-section and subjected to a pair 
of compressive longitudinal loads 1 P 1 acting along the 
centroidal axis, as shown in Figure 4.1(a). 

As is understood, from the we 11 -known-theory of 
elastic stability given by say, Timoshenko (1961), the member 
becomes structurally unstable at a particular value of load 
P, called the critical or buckling load 'P cr ' . Adopting a 
two-element discretization of the member [Figure 4.1(a)] , the 
generalized nodal displacements {q ± } and the corresponding 
internal nodal forces {Q ± } can be seen to be as shown in 
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Figures 3.3(b) and 3.3(c) respectively, it may be noted that 

{q } can be expressed in the form of subvectors {q, } , {q ) 
x A B 

and {q >/ as indicated by Equations 3.4 and 3.5, and {Q.} by 

^ i 

subvectors {Q^}, {Qg) an< 3 } , as shown by Equations 3.6 and 

3.7. Obviously (q^) and {q^} are governed by the prescribed 

boundary conditions. 

Recalling from Article 2.11.2, for a given set of 
boundary conditions, the governing equation to determine the 
unknown displacements is the eigenvalue problem 

(Ojj] - x [g ± j]) = 0 ... 2.42 

[k. .] and [g. .] being the stiffness and the stability matrix 
i j i j 

respectively of the member. The non -trivial solution for 

{q^} can be obtained by determining A as an eigenvalue of 

the matrix fg. .I” fk. .1. The value of A is directly 
XJ L 1J 

related to the critical load ? cr , while the vector (q.^ } 
obtained corresponding to a value of 1 represents the parti- 
cular buckling mode of the member. 

4.2.1 Buckling Modes 

Three possible modes exist in which centrally 
loaded columns can buckle : - 

(a) they can bend in the plane of one of the principal 
axes (flexural mode ) , 

(b) they can twist about the shear centre axis 
(torsional mode) s 



or, (c) they can bend an.d twist simultaneously 
( tor sional- ■ flexural r.icJe) . 

For any given member, depending on its length and 
the geometry of its cross-section, one of these three modes 
will be critical. 

Referring to Chajes (1965), for an axially loaded 
column of thin-walled cross-section and having any boundary 
conditions, the critical load ? cr is given by the roots of 
the general equation in 'P', for the torsional-flexural 
buckling 

r 2 (P - P )(P - P )(P - P m ) - P 2 a 2 (P - p ) - p 2 a 2 (p - P ) = 0 

x y <p y y x x 

... 4.1 


where the parameters P , P and P have the form of the two 
c x y <p 

pure buckling loads in planes of the principal axes and a 
purely torsional buckling load about the shear centre axis. 
They are '-given by 


t 2 EI 


XX 


... 4.2(a) 


2 

n El. 


'JDL 


... 4.2(b) 


and 


9 


(GK + 


EI fiP * 


2 

2 } 


... 4.2(c) 


where L is the effective length of the member, 1^. 3nd Iyy 


are given by Equation 2.3(a), K by 2.3(b), I 


S2S2 




by 2.3(c), 
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In the present case, the section is symmetric with 
respect to one axis (x=axis) so that setting a y = 0 reduces 
Equation 4.1 to 

(P - P x ) [r 2 (P - p y )(P - P ) - P 2 a 2 ] = 0 ... 4.3 

There are again three solutions, one of which is P = P which 

represents -purely flexural buckling i n x-z plane [Figure 

4.1(a)]. The other two are the roots of the quadratic term 

inside the square brackets equated to zero and give two 

torsional- flexural buckling loads. The lowest torsional- 

flexural load will always be below p y or P . It may, however, 

be above or below p depending upon the sectional dimensions. 

X 

As mentioned by chajes (1965), for sections having 

T > i , the lowest torsional-flexural buckling load is 
xx yy 

always less than P . Referring to Table I of SP ; 6(1) (1964), 
it may be seen that the entire range of available rolled 
steel I -sections belongs to the above category and hence the 
governing buckling mode will be that in torsional flexure. 

This also emphasises the need to consider the sectional 
warping in the analysis, that arises as a result of the 
flexural torsion. 

In order to determine the buckling load P = P cr 

from Equation 2.42, the elements of matrices [k ±j ] and [g^], 

for an element, can be written in the non-dimensional form 

El -n 

talcing out terms — ^ end | , outside these matrices respect- 
ively. It may be noted that elements of [k^] correspond to 
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flexure in the y-z planes {Figure 4.1(a)] and the term 
El 

— ensures to give the lowest buckling load, 1^ being less 
than I xx * Thus, for an element in Figure 4.4(a), Equation 
2.42 can be written as 
EI 

(" — 3 ^ C^. . ] = Y [g . . ]) {q. ) = 0 ... 4 .4 

i-J - 1 - ij i 

This can be expressed as the general eigenvalue 
problem, as suggested by Equation 2.42 


[k. .] (q.) = X [g. .] (q.} ... 4.5 

1J i lj x 

■■ - ’ pl_2 

where X = — — ... 4.6 

EI 

yy 

It may be noted that the matrices [k. .] and [g . . ] above, 

1 J X J 

are square, -positive definite and symmetric . 

The general eigenvalue problem is solved after 
converting it into the standard form, given by 


[b] (q^ } = X {q 1 } ... 4.7 

where [g^j] = [ L ] M ... 4.8(a) 

obtained after the standard cholesky decomposition, and then 

[B] = fL" 1 ] r^j] [L- 1 ] T ... 4.8(b) 

and (q-j ) = [b] T (qp 


. . 4.8(c) 
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4. 2. 1.1 Number of critical loads 


If ndf is the number of degrees of freedom in the 
system, as governed by the boundary conditions, the elements 
of matrices fkn and [g^j] in Equation 4.5 correspond to 
the free unknown nodal displacements. Hence these matrices 
are of the order (ndf x ndf), while the vector {q }has size 
(ndf x 1). The eigenvalue problem has 'ndf 1 solutions, each 
corresponding to a critical load on the column, it may be 
noted that all solutions are positive by virtue of the 
positive definiteness of matrices fk. .] and fa .1* 

^min sma H- es t eigenvalue, the critical 

load to be considered for the design is given by 


P 


cr 


X 

min 



4.9 


and the corresponding value of buckling stress being critical, 
is 


^bcr 


P 

_cr 

A 


. . . 4.10 


4.3 CONSTRAINTS IN THE OPTIMIZATION PROBLEM 

Out of the fourteen constraints prescribed in the 
problem, twelve are the side constraints, the remaining two 
being the behaviour constraints. 

Recalling that the design vector is 
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the constraint relations can be obtained as follows. 


... 3.1 


4.3.1 side Constraints 


4. 3.1.1 Lower and upper bounds on the design variables 

The twelve relations representing the lower and 
upper bounds on the design variables are the same as those 
given by Equations 3.14(a) through 3.14(1). 


4.3.2 Behaviour Constraints 

4. 3. 2.1 Maximum slenderness ratio in the member 

If ' SR 2 ‘ is the slenderness ratio of the column 
corresponding to the appropriate radius of gyration of the 
section, in accordance with the Clause 19.2.1 and Table XVI 
of IS : 800 (1962), this ratio should not exceed 180. So 


sr 2 

Q 13 ( x ) = 180 ” 1 


... 4.11(a) 


4. 3. 2. 2 Maximum compressive stress in the column 

Inspection of Equation 4.10 shows that, the maximum 
compressive stress that one may expect at the working loads 


is 
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4.4.4 Height of the column 

The range of column heights chosen is from 3.0 to 
6.0 m for the fixed supports, 3.0 to 4.5 m for the simple 
supports and from 3.0 to 4.0 m for the cantilever. 

The effective height can be calculated by multi™ 
plying the actual height H by a coefficient, whose value in 
the respective cases above is 0.67, l.o and 2.0, as given by 
Table XV of IS : 800 (1962). The appropriate value of 
slenderness ratio ' SR^ 1 / thus calculated, is used in 
Equation 4.11(a). 


4.4.5 Lower and Upper Bounds on the Design variables 

The values o£ tf^ In 

Equations 3.14(a) through 3.14(f), and bf^ m , tf^ , ..., 
tW max Ec J uat i ons 3.14(g) through 3.14(1) are the same as 
given in Articles 3.5.9 and 3.5.10 respectively. 


4.4.6 Maximum Outstand of Flanges 


As per the Clause 19.4.1 of IS ; 800 (1962), the 
maximum outstand of the compression and the tension flange 
is sixteen times the respective flange thickness. This 


condition indirectly sets an upper bound on the values of 

b£ lmax 314 bf 2max- The values ° £ “lmax 314 bf 2max to ^ 
used in Equations 3.14(g) and 3.14 (i) are the lesser of 

those given by Articles 4.4.5 and 4.4.6 respectively. 
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4.4.7 Maximum Depth of the web Plate 

As per the clause 19.4.2.1 of IS : 800 (1962), in 

computing the effective area and radius of gyration, the 

depth of the plate should be reckoned as not more than fifty 

times the web plate thickness, in a minimum weight design, 

the total depth of web plate hw may be subjected to this 

limit. The value of hw _ to be used in Equation 3.14 (k) is 

max 

the lesser of that given by Articles 4.4,5 and 4.4,7 respect- 
ively. 


4.4.8 Permissible value of the compressive Stress at Buckling 

The value of 'FS' to be used in Equation 4.11(b) is 
4 and the value of cr b p er in Equation 4.11(c) is given by the 
value of *P ' obtained by using Table II of IS : 800 (1962). 


4.5 OBJECTIVE FUNCTION 


The objective function is given by 


f (X) = x 1 x 2 + x 3 x 4 + x 5 x 6 


... 3.15 


4.6 OPTIMIZATION PROBLEM 

The present constrained optimization problem is 

given as , 


find 
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s'o as to minimize the function 


f (X) = 


x l x 2 


+ 

5 6 


subject to the inequality constraints 
bf . 


9l (X > “ 


' lmin 


x. 


1x0 


tf 


g 2 (x) = 


lmin 
x^ ~ 


lxo 


3 1 

» « » • * j - 


... 3.15 


. .. 3.14(a) 


... 3.14(b) 


and , 


13 


J 14 


» * 


(X) 

6 

tw 


max 


SR 2 

(X) 

180 " 1 


^bw 

(X) 

Cf 


bper 


1^0 


0 


0 


... 3.14(1) 

... 4.11(a) 

... 4.11(c) 


The solution of this problem is sought by using the 
interior penalty function method, as discussed in Articles 
3.8 and 3.9. 


4.7 COMPUTER PROGRAMME 

The computer programme prepared for obtaining the 
following two programmes; 


results comprises 
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(i) STBANA -- which deals with the stability analysis 

part, 

and, (ii) 0PT3 — which deals with the optimization part 

based on the results supplied by 
STBANA. 

The functions performed by various subroutines of 
the programme STBANA are as follows; 

(i) CSPROP — calculates the cross-sectional properties 

of a thin -walled I -section 

(ii) ELSTIE — Generates the stiffness matrix for a thin- 

walled beam element 

(iii) ELSTAB — Generates the stability matrix for a thin- 

walled beam element 

(iv) ASSBLY — Assembles the stiffness matrix [kjj] and 

the stability matrix [g^] from the 
individual element matrices. It further 
isolates submatrices from [k.^] and [g^j], 
corresponding to the free displacements 
in the structure 

(v) EIGSTD — — Converts the general eigenvalue problem 

into the standard one 

(vi ) CHOSKY — — performs the cholesky decomposition of a 

symmetric, positive definite, square matrix. 

(vii) MATMLT — Obtains the product of two matrices 

performs the matrix inversion 


(viii ) MATIN V 
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(ix) EIGEN 


Finds the eigenvalues of a real, symmetric 
matrix by the Givens -Householder method 


and , 


(x) STBANA 


For a thin -iv ailed beam of I— section of 
known sectional dimensions and boundary 
conditions, it calculates the least value 
of the critical load, and the corresponding 
critical stress at buckling, by using a 
discrete element analysis of the structure. 

The optimal solution is obtained by the programme 
0PT3 using following subroutines; 

(i) MAIN — Reads the data for the problem, performs 

various steps in the optimization process 
and finally returns the optimal sectional 
dimensions and the corresponding values of 
the objective function and the constraints 

(ii) SUMT3 — Obtains the constrained minimum of a given 

objective function using the sequential 
unconstrained minimization technique (SUMT) 

(iii) PELTY3 — » Generates the objective function and the 

penalty function for a design point X 
Gives the allowable stress in axial compre- 
ssion as per Table II of IS ; 800 (1962) 
Obtains the uncon strained minimum of a 
function by the Davidon-Fletcher-Powell 


(iv) TABLE 2 


(v) DFP3 


method 
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(vi) UNIDR 3 


(vii) GRADF3 


Minimizes a function in one direction using 
the method of cubic interpolation 
Gives the gradient vector and its modulus 
a t a given point X in the function/ by using 
the finite difference technique. 


4.3 DYNAMIC ANALYSIS 


Consider a prismatic member of length L having a 
thin -walled , l~shaped cross-section and subjected to a given 
set of boundary conditions. Assuming that the member is 
displaced laterally by a small magnitude from its initial 
equilibrium position, it will vibrate in transverse direction. 

Adopting a two -element discretization of the member 
[Figure 3.3(a)] , the generalized nodal displacements {q^ } 
[Figure 3.3(b)] are expressible as sub-vectors {q^ - }, {q B } and 
{q C } by Equations 3.4 and 3.5, where (q A } and {q B } are 
governed by the prescribed boundary conditions. 

Recalling from Article 2.11.3, for a given set of 
boundary conditions and with no axial load 'P‘ on the member, 
the governing Equation 2.44 is reduced to 




(4JD 


[m, .] ) tq, > = 0 


ij 


... 4.12 


where fk. .] and Tm. .] are the stiffness and the mass matrix 
1 aj J L ij J 

respectively. Equation 4.12 represents a general eigenvalue 
problem that gives the natural transverse vibrational 
frequencies '*»' and the corresponding nodal displacements 
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{<1^}/ assumed as unknowns in the problem, one may note the 

similarity between Equations 2.42 and 4.12 and conclude that 

the buckling and vibration phenomena are closely interrelated . 

The non “trivial solution for {q^} can be obtained 
2 

by first determining a as the eigenvalue of the matrix 
r m i jl [kij] * The vector {q^} obtained corresponding to a 

value to, represents a particular vibrating mode of the 
structure . 

4.8.1 Vibrating Modes 

By virtue of the similarity between Equations 2.42 
and 4.12, one can conclude with the help of Article 4.2.1 
that the section being symmetric about x-axis [Figure 4.1(b)] , 
the three possible modes of transverse vibrations are: 

(a) a purely flexural mode in the x-z plane [Figure 
3.3(a)] 

and (b) two torsional-flexural vibrating modes. 

Since the cross-section has I xx > I / the 
governing vibrating mode is expected to be a torsional— 
f lexural one . 

In order to determine the lowest vibrating 
frequency, using I (being less than l xx ) to express the 
elements of matrices [k^j ] and [rru j] in the non-dimensional 
form, for an element of the discretized model. Equation 4.12 


takes the form 
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= O ... 4.13 

of the beam material, 
eigenvalue problem 

... 4.14 

... 4.15 

It may be noted that the matrices [k. .1 and fm. .1 

ij ij 

above are square > positive definite and symmetric. 

As discussed in Article 4.2.1, the general eigen- 
value problem given by Equation 4.13 is solved after converting 
it into the standard form. 

As mentioned in Article 4. 2 .1.1, Equation 4.14 has 
1 ndf ' eigenvalues, all being positive by virtue of the 
positive definiteness of matrices [k. .] and [m. .], where 'ndf' 

1 J i J 

is the degree of freedom of the system. Each solution gives 
a natural vibrational frequency and the corresponding mode of 
transverse vibration. Let be the smallest value of the 
frequency (corresponding to the first mode) of transverse 
vibrations . 

4.9 CONSTRAINTS IN THE OPTIMIZATION PROBLEM 

Out of the fourteen constraints prescribed in the 
problem, twelve are the side constraints, the remaining two 
being the behaviour constraints. Rewriting the design vector 


( ”73 y - [k ij ] " “ 2 PA1 Cm.j]) {q.} 

where ’p 1 is the mass per unit volume 
This can be expressed as the general 

t k ij ] <v = » • Kj ] {q i } 

where \ = . 2 . HA— 
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the constraint relations can be obtained as follows. 


3.1 


4.9.1 Side constraints 

4. 9. 1.1 Lower and Upper Bounds on the Design variables 

The twelve relations representing the lower and 
upper bounds on the design variables are the same as those 
given by Equations 3.14(a) through 3.14(1). 


4.9.2 Behaviour Constraints 

4. 9. 2.1 Lower and Upper Bounds on the First Natural Frequency 


If m . and u represent the lower and upper bounds 
mm max 

on the first natural frequency/ 


g 13 00 



• • * 4 * 16 (a) 


g i4 (x) 



m 

max 


... 4.16(b) 


4.10 DATA FOR THE ANALYSIS 

The number of elements, the boundary conditions in 
the problem and the values of elastic constants are those 
given in Articles 4.4.1, 4.4.2 and 4.4.3 respectively. 
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4.10.1 Span of the Beam 

The range of beam spans chosen is from 3.0 to 6.0 m 
for the fixed and the simple supports and 1.0 to 2.0 m for 
the cantilever. 


4.10.2 Lower and Upper Bounds on the Design variables 

The values of M lmln , in Equations 

3.14(a) through 3.14(f), and bfj^, tf^ tw^ in 

Equations 3.14(g) through 3.14(1) are the same as given in 
Articles 3.5.9 and 3.5.10 respectively. 


4.10.3 Maximum Outstand of Flanges 

The values of maximum outstand of flanges which set 
upper bounds on bf^ and bf^max' are those given by Article 

4.4.6. 


4.10.4 Maximum Depth of the Web Plate 

The value of maximum depth of the web plate which 

sets an upper bound on hw is that given by Article 4.4.7. 
r • max 

4.10.5 Lower and Upper Bounds on the First Natural Frequency 

The range of spans chosen (Article 4.10.1) and the 
available rolled steel sections as given by Table I of 
SP : 6(1) (1964), indirectly decide the lower and upper 
bounds on the first natural frequency. Two different ranges - 
a narrow and a broad one, are chosen for each type of boundary 


conditions . 
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For fixed supports, the values of to . and ® 

mm max 

are (100, 200) and (50, 300) respectively for the two 
different ranges. Similar values for the simple supports 
(75, 125) and (30, 150) and for the cantilever these are 
(100, 200) and (50, 450) respectively. These values have 
used in Equations 4.16. 

4.11 OBJECTIVE FUNCTION 

The objective function is given by 

f(X) = x i x 2 + X 3 X 4 + x 5 x e 


4.12 OPTIMISATION PROBLEM 

The present constrained optimization problem is 

given by; 

to find 



s o as to minimize the function 

f (X) = x t x 2 + . .. + X 5 X 6 


are 


been 


3.15 


3.1 


3.15 


subject to the inequality constraints 
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®l (x) 

_ ° f lmin 

*1 

- 1 v< 

0 

• • • 

3.14(a) 

g 2 (x) 

tf 

lmin 

x 2 

- 1 <5 

0 

• • • 

3.14(b) 

gi 2 <x) 

tw 

~ 1 *< 

• m 

0 

• • • 

3.14(1) 

g 13 S) 

max 

( %in 

(0)i 

* 1 < 

0 

* • « 

4 .16(a) 

and, (X) 

0)0 1 

- 1 < 

0 

m • * 

4 .16(b) 

ffli max 





The solution 

of this problem 

is sought by 

using 

the 


interior penalty function method as discussed in Articles 3.8 
and 3.9. 


4.13 COMPUTER PROGRAMME 

The computer programme prepared for obtaining the 
results comprises the following two programmes: 

(i) FRQANA -= which deals with the dynamic analysis 

part, 

and (ii) 0PT4 — which deals with the optimization part 

based on the results supplied by FRQANA. 
The programme FRQANA comprises subroutines CSPROP , ELSTIF , 
ELMASS , ASSBLY, EIGSTD , CHOSKY, MATMLT , MATINV , EIGEN and 


FRQANA 
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The functions of these subroutines, except those 
given below, are described in Article 4.7. 

ELMASS — Generates the mass matrix for a thin-walled beam 
element 

ASSBLY — Assembles the stiffness matrix [k. .1 and the mass 
matrix [m^.J from the individual element matrices. 
It further isolates submatrices from [k^] and 
[m. .] corresponding to the free displacements in 

X. J 

the structure . 

FRQANA -- For a thin-walled beam of I-section of known 

sectional dimensions and boundary conditions, it 
calculates the first natural frequency of trans- 
verse vibration, using the discrete element 
analysis of the structure. 

The optimal solution is obtained by the programme 
0PT4 using subroutines MAIN, SUMT4 , PELTY4, DFP4 , UNIDR4 and 
GRADF4. The functions of these subroutines are the same as 
those of MAIN, SUMT3, PELTY3 ,..and GRADF3 respectively and 
described in Article 4.7. 

4 . 14 CHARACTERISTICS OF THE PROGRAMMES 

(a) The programmes STBANA and FRQANA are quite general and 
can yield the two—, four— and eight— element analysis 
of the thin— walled beams having an I— section, with 
different spans and boundary conditions . 
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(b) As discussed in para (b) of Article 3.10.1/ the code 
numbering used to generate any boundary conditions/ 
enables one to handle the continuous beams as well. 

(c) Programmes 0PT3 and 0PT4 incorporate logical provisions 
mentioned in paras (c) and (d) of Article 3.10.1/ so 
that they are made numerically stable. 

The results for both these problems have been 
discussed and the conclusions drawn in chapter 6. 


- o 
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5. OPTIMAL DESIGN OP A CANTILEVER RETAINING WALL 


5 . 1* GENERAL 

A wall designed to maintain a difference in tlie 
elevations of the ground surfaces on each side of the wall is 
called a detaining watt. The earth , whose ground surface is 
at the higher elevation is commonly called the back. fit l and 
the wall is said to retain this backfill, in general, the 
term backfill refers to the material behind a wall, whether 
undisturbed ground or fill, that contributes to the pressure 
against the wall. 

The side of the wall adjacent to the backfill is 
called the back and the side which is exposed for most of 
its height is called the face . The bottom surface of the 
wall is called the base or the footing. The line of inter- 
section of the back and base is called the heels a ^d the line 
of intersection of the front of the wall and the base is 
called the toe. Also, a projecting portion on the front of 
the wall is called the toe projection and one on the back, 
the heel projection . The batter of the face or back of a 
retaining wall is its inclination with the vertical. 

Earth, which lies above a horizontal plane at the 
elevation of the top of a wall, is called a surcharge . if 
the ground surface slopes upward from the top of a wall as 

* Articles 5.1 through 5.7, although of elementary nature, 
are presented here for the sake of completeness. 
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shown in Figure 5.1(a), the wall -i c oa-ta 4 . -u, . 

v ' ' uie Wd -i-L is said to have an inclined 

or sloping surcharge . 

Retaining walls are extensively used in connection 
with railways,, highways, bridges, canals, and many other civil 
engineering works . 

5.2 TYPES OF RETAINING WALLS 

The shape of the cross-section, or profile of a 
retaining wall may be determined by the prevailing conditions, 
but usually retaining walls are of four general types; the 
solid gravity wall, the semigravity wall., the cantilever or 
T—wall, and the counterfort wall. The buttressed wall is 
used occasionally under special conditions. Retaining walls 
are mostly constructed of plain or reinforced concrete, 
although plain or brick masonry is sometimes used. 

Amongst all these categories, a cantilever type of 
wall is more widely used than any other type, in spite of the 
durability advantages of the solid gravity and semigravity 
walls, because it is usually more economical for ordinary 
heights, say from 3 to 6 m. A counterfort type of wall may 
prove to be more economical for greater heights. 

5.3 TYPES OF CANTILEVER WALLS 

Common types of cantilever walls are illustrated 
in Figures 5..1(b), (c) and (d). This type of wall is seen 
to be made up of three cantilever beams: the vertical stem. 
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the toe projection and the heel projection, it is also 
called a T-wall because it resembles an inverted *T'. The 
simplest form shown in Figure 5.1(b) is suitable for small 
heights. For higher walls, it may prove economical to provide 
fillets in the comers as shown in Figure 5.1(c). The 
sliding resistance may sometimes, but not always, be increased 
by a downward projection below the base as shown in Figure 
5.1(d). This projection, called a key , is usually placed 
near the centre of the base so as to provide end embedment 
for the vertical reinforcement in the stem. It is more 
effective in resisting sliding if placed at the heel; it is 
sometimes placed at the toe as well. A construction joint 
is provided at the junction of the vertical stem and the top 
of the footing. 

The present optimization problem deals with canti- 
lever walls of ordinary range of heights (3 to 6 m) and so, 
the type shown in Figure 5.1(b) is chosen for this purpose. 

5.4 FORCES ACTING ON RETAINING WALLS 

• The forces acting on a simple form of retaining 
wall, where the cross -sect ion remains constant, are shown in 
Figure 5.2(a). For the sake of simplicity, a trapezoidal 
wall is shown in the illustration. 

The various forces that can be considered for the 
analysis of a retaining wall are as follows; 
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(a) The weight of the wall. The resultant of this weight/ 

* 

of course, acts through the centre of gravity of* the 
portion above the horizontal section considered, as 
indicated by 'W' in Figure 5.2(a). usually, it is 
simpler to divide a cross -section into triangles and 
rectangles, whose areas and centroids can be determined 
easily and to deal with the corresponding forces 
separately. 


(b) Th-~ earth pressure against the back of the wall. The 

pressures usually considered are those due to the lateral 
thrust of the retained earth, as indicated by the 
pressure diagram and the resultant force P acting 
against the back A.B with an obliquity 6 as shown in 
Figure 5.2(a). This is called active pressure. 


(c) The normal components of the foundation pressures . 

These are usually considered as being distributed linearly 
across the base CB giving a trapezoidal pressure diagram, 
as shown in Figure 5.2(a). The resultant of these forces 
acts through the centroid of the pressure diagram and 
is indicated by XV. The base of a wall is usually 
■ horizontal, so that the normal pressures, called the 
foundation or bearing pressures s are usually vertical. 


(d) The horizontal components of tht, feundatzon pr„ss 

The resultant of these components is indicated by XH f 
in Figure 5.2(a). The distribution of these pressures J 
is usually not considered but their resultant is considered. j 
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(e) The earth pressure against the fees of the fall. The 
bottom ox a retaining wall must be placed at a depth 
wherein the soil with the required bearing capacity can 
be reached. The earth in front of the wall may exert 
resistance to forward movement against the buried portion 
of the wall. This is known as passive pressure or 
passive resistance as indicated by the resultant force 
P' acting at an obliquity 6', as shown in Figure 5.2(a). 
The force P' adds to the stability of the wall but is 
omitted from the computations because of uncertainties 
concerning its magnitude. 

The present analysis considers the foregoing forces 
(a) through (e). However, mention may be made of the following 
additional forces that arise depending on the situation and 
the purpose for which a retaining wall is constructed; 

(f) Bridge reactions, tractive effect and centrifugal 
force 

(g) Surcharge loads 

(h) Water and seepage pressures 

(i) uplift 

( j ) Vibration 

( k ) Impact 

and (1) Earthquake effects. 
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5.5 PROPERTIES OF SOIL 

So far as the earth pressure computations are 
concerned, the most significant properties of soil are as 
follows . 

(i) Unit weight of soil (w ) : With various moisture 
contents and the degree of compaction. 

(ii) Sheaving strength of soil (S): The shearing strength 
or resistance of soils is related to the design of 
retaining walls in two ways: A wall is subjected to 
lateral earth pressures exerted by a backfill of 
disturbed soil. The magnitude of these pressures depends 
upon the shearing resistance. Moreover, the wall is 
usually supported by a foundation of undisturbed soil. 

The stability of this support against sliding over the 
foundation, depends on the shearing resistance. It is 
therefore necessary to have information concerning the 
shearing resistance of soils in both the disturbed 

and undisturbed conditions. 

Based on a series of direct shear tests for most 
soils, the shearing strength may be expressed by the coulomb- s 

equation 

S = c + n ta n0 **• 

where s is the unit shearing strength of the soil, n is the 
unit normal pressure on the surface of failure, c and 0 are 
the cohesion and the angle of internal friction respectively 


of the soil 



144 


Clean granular materials such as gravel, sand etc., 
whether saturated or dry, are designated as cohesionless and 
for these, the value of c is taken as zero. Silts and clays 
are categorized under cohesive soils. In the present analysis, 
while investigating the stability of the wall against sliding, 
appropriate values of c and 0 are used in Equation 5.1, taken 
as the basis . 

There are other soil properties like stress strain 
relationship « volume change due to moisture contact, 
permeability etc. but these are not related to the present 
analysis . 

5.6 CLASSES OF EARTH PRESSURES 

The lateral pressures in a soil mass are greatly 
affected by any lateral deformation which the mass may have 
experienced. Such pressures are divided into three classes: 
earth pressure at rest, or natural earth pressure; active 
earth pressure; and passive earth pressure. 

5.6.1 Natural Earth Pressures 

The lateral earth pressures in a natural soil 
deposit which has remained undisturbed are referred to as 
natural earth pressures or earth pressures at rest . 

Lateral earth pressures in deposits formed artificially by 
some construction operation and in which no lateral defor- 
mation has occurred subsequent to deposition are called 
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earth pressures at rest. The magnitude of these pressures 
depends upon the manner in which a soil is deposited as well 
as upon the physical properties of the soil. 


5.6.2 Active Earth Pressures 

These exist in a soil deposit which has been extended 
laterally, such as the backfill of a retaining wall, founded 
on soil which has tilted or moved horizontally a small but 
sufficient amount due to the pressures exerted by the backfill 
and the yielding of the foundations. 

5.6.3 passive Earth Pressures 

These exist in a soil deposit, which has been 
compressed laterally by the movement of a structure which 
abuts against it, or in some other manner. The magnitude of 
the lateral compression of the soil required to develop the 
passive state of stress in a given soil is somewhat greater 
than that of the lateral extension required to produce the 
active state of stress in the soil. Passive pressure is 
often called Passive resistance . 

5.7 EARTH PRESSURE PHENOMENA 

All theories and methods for determining lateral 
earth pressure due to cohesionless or cohesive soil can be 
grouped under Rankine 1 s or coulomb 1 s methods according to 
their basic assumptions. As given by Huntington (1957), 
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whether the soil is cohesive or cohesionless, the conditions 
which govern the application of Rankine's or Coulomb's method 
may be summarized as follows; 

(a) if the back of the wall is a plane surface or a surface 
which can be assumed to be plane without introducing 
significant errors, either method may apply. 

(b) If the back of the wall cannot be assumed to be a plane 
surface, coulomb's condition cannot prevail but Rankine's 
conditions prevail. 

The backs of cantilever and counterfort walls 
depart so widely from plane surfaces that Rankine's method 
can be usually applied without significant error, for the 
cohesionless and the cohesive soils. As mentioned in Article 
5.4, the backfill on the retaining wall is under Rankine ' s 
active state of stress while the soil mass abutting against 
the face is under Rankine's passive state of stress [Figure 
5.2(a)], This passive resistance adds to the stability of 
the wall but is omitted hereafter in the computations because 
of uncertainties concerning its magnitude. 

5.7.1 Active Pressures - Rankine's Conditions 

Referring to Figure 5.2(b), the lateral earth 
pressure p e in the Rankine ' s active state of stress and on 
a vertical plane AV at depth h, is given by Terzaghi (1960) 
as 
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P 


e 


w h [cos i 

o 


cos i — /cos^i «• cos^# 


/ 2 '2 
cos 1 + y cos i - cos 0 


■J 


5.2 


= w g h K say, 


where is the unit weight of soil and i, the inclination 
of surcharge to the horizontal . k represents the quantity 
in the parantheses in Equation 5.2. 

The resultant active pressure on the vertical plane 

AV is 


„ 1.2 
P = - w h K 

2 s 


5 3 

• • # —> • 


acting at the lower third point of the vertical plane and 

acting -parallel to the ground surface. 

If the ground surface is horizontal, i = 0 and the 

1 s inti 

factor K in Equation 5.3 reduces to ± lT~sin$ 


1 u 2 1 - sin# 

p = — W h zr— r~-% 

2 s 1 + sm0 


5.4 


5.8 PRELIMINARIES IN THE OPTIMIZATION PROBLEM 

Consider a typical section of the cantilever 
retaining wall chosen for the optimization problem, as shown 
in Figure 5.3(a). The optimization problem is handled on 
the basis of the following assumptions: 

(i) The face of the wall is provided with a batter, 
keeping the back vertical. 
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(ii) The backfill is assumed to be dry and having a 

horizontal top level; in other words, no surcharge is 
assumed in the problem. 

(iii) The cantilever retaining wall is of reinforced concrete 
whose design is governed by the Indian Code of Practice 
for plain and reinforced concrete, IS : 456 (1964). 

(iv) No tension is to be developed in the foundation. 

(v) 12 mm diameter plain, mild steel bars are provided for 

the main reinforcement in the stem and toe slab and 
16 mm diameter bars in the heel slab. The distribution 
steel is of 10 mm diameter bars. Nominal temperature 
reinforcement is provided in the wall. 

The preassigned parameters in the optimization 
problem are the top width b and the overall height h of the 
wall [Figure 5.3(a)], The design variables are: the stem 
thickness at the bottom, x 1? the thickness of the base slab, 
x 0 ; the total base width, x- and the toe projection, x 4 so 
that the design vector X is expressible as 



Recalling Article 5.4, the various forces that 
would form a critical set in the analysis and subsequent 
design of the present optimization problem, ere shown in 
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Figure 5.3(a). it can. be easily noticed that for a 'no- 
tens ion “-condition ' in the foundation, the critical pressure 
distribution on the base, assumed to be linear, has to be 
triangular in shape with the maximum ordinate value 'p' at 
the toe . 


5.9 CONSTRAINTS IN THE PROBLEM 


The present optimization problem is governed by 
eighteen constraints, ten of which are the behaviour const- 
raints and remaining eight, the side constraints. 

Expressing the design variables x^, . . . , x^ and the 
preassigned parameters b and h in metres, the forces and 
moments are evaluated as kg and kg-m respectively and per 

unit length of the wall, taken as one metre. The areas 

2 3 

and volumes are in m and m respectively, unless otherwise 
mentioned . 

Following quantities can be evaluated first: 


the cross-sectional area A 


soil parameter K 


|(b + x 1 ) (h - x 2 ) + x 2 x 3 

... 5.6 

1 ~ sin 0 c -7 

TTT&& * * * 


0 being the angle of internal friction for the soil retained. 


The resultant lateral earth pressure over height h 

1 2 
= — w .h K 

2 s 


P 


» • • 


5.8 
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acting at a height | from the base, where w is the unit 

s 

weight of soil retained. 


5.9.1 Behaviour Constraints 

5. 9. 1.1 Minimum stem thickness at the bottom to resist 
shearing force 

The external horizontal shearing force on the stem 
bottom section 

F 1 = I K • » s (h ' * 2 > 2 

The internal resisting force in reinforced concrete 


rl 


= j(x- 


c s } 


where j is a constant associated with the permissible stresses 
in the concrete and steel used, c s is the cover provided for 
main steel in the stem and t c the permissible shear stress in 
concrete . Hence 


g x 00 



... 5.9(a) 


5. 9. 1.2 Minimum stem thickness at the bottom to resist 
bending moment 


External bending moment at the stem bottom section 


«i = r <h ' x 2> 

The internal moment of resistance for a balanced section is 
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M 


rl 


= Q (x = c )' 


where Q is a constant associated with the permissible stresses 
in the concrete and steel used. So 


g^(x) 


M rl 


... 5.9(b) 


5. 9. 1.3 Safety against overturning about the toe 

Subdividing the cross-section into rectangles and 
triangles [Figure 5.3(a)], the weights of different portions 
are 


W. = w b(h - x ) . 
1 C 2 


W, 


2 w^ (x.^ - b) (h - x 2 ) 


Vi- = w x^ x_ . 
3 c 3 Z 


being the unit weight of reinforced concrete. 

The weight of soil over the heel portion is 

W 4 = w g (h - x 2 ^ X 3 ” X 4 " X l^ * 

At the time of overturning, the resultant base 
pressure passes through the toe and therefore has no moment 
about the toe. The overturning moment is due to P only and 


M, 


iph 


is given by 
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The balancing moment by the other forces is 

k ' 2 , 


M 


r2 


- W 1 (x 4 + x 1 - --) + W 2 [x 4 + |(x 1 - b)] + 


Ha, 

2 3 


W 4 , 

+ T <X 3 


X 3 " X 1 “ x 4, 


If 1 FS q ' is the factor of safety against overturning/ 


(FS ).M 

g (X) = ■--- A - 1 ...5.9(c) 

M r2 


5. 9. 1.4 Safety against sliding 

Let W = w 1 + W 2 + w 3 + w 

The external sliding force is F 2 = P. 

The frictional resisting force by virtue of Equation 

5 . 1 is 

F r2 = c . x 3 + W tan0 

If 'FS 1 is the factor of safety against sliding, 
s 

then 

(FS )F 2 

g A (X) = — 1 ••• 5 * 9 ( d ) 

4 F r2 


5. 9. 1.5 No-tension condition in the foundation 

If R is the resultant foundation pressure on the 
base , then 

R = I • X 3 
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acting through the centroid of the pressure triangle as 
shown in Figure 5.3(a). considering the equilibrium of the 
given set of forces, one can note that R = w, and the eccen- 
tricity Z of the resultant R with respect to the base centre 
can be expressed by the equation 


M 


r2 “ + which yields 


Z = 


*3 “r2 ;JJ 2 

2 W 


x. 


In the 'no-tension condition 1 , Z is limited to a 


value -g— as seen from Figure 5.3(b) and so 


gc(x) 


x 3 76 


... 5.9(e) 


5. 9. 1.6 Maximum pressure on the foundation 

Referring to Article 5. 9. 1.5, the maximum pressure 
p at the toe is given by 


P = 


2R 

x„ 


2W 

X-p 


If Pp 0r is the permissible bearing capacity of the 


f ound at ion , then 


g 6 (x) = - 


per 


... 5.9(f) 
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5. 9. 1.7 Minimum thickness of the toe slab to resist shearing 
force 

The external shearing force at section A-A [Figure 
5.3(a)] is 


px. x, - X , 

F 3 = — (It — “) - V4 X 2 


The internal resisting force in reinforced concrete 
is 

F r3 * J(x 2 - c b > T c 

where c^ is the cover provided for the main steel in the base 
slab. Therefore 


g 7 (x) = 


... 5.9(g) 


ri 


5 . 9 . 1.8 Minimum thickness of the toe slab to resist bending 
moment 

The external bending moment at section A--A [Figure 
5.3(a)] is 


(x, -x,) x; n w c „2. 


M. 


- fe -H-- 1 - • r + 1 • P ] " t- -4-2 


x.x. 


.3 

: ± 

‘3 “ : 3 ' 

The internal moment of resistance of the balanced section is 


M r3 = Q^ x 2 “ C b r 


Hence , 


g 8 (x) = 


M. 

M* 


r 3 


5.9(h) 
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5.9.2 Side Constraints 

5. 9. 2.1 Lower and upper bounds on the design variables 

Using subscripts 'min' and 'max’ to denote the 
prescribed minimum and maximum values of variables x^/ . . . / x^ 
the constraint relations giving the lower and the upper bounds 
on the design variables are 


g u (x) 

= 

X lmin ^ 

X 1 

... 5.9(k) 

g 12 (x) 

= 

X 2min ^ 

X 2 

... 5.9(1) 

913® 

= 

X 3min = ^ 
x 3 

... 5 .9 (m) 

Si4® 

= 

X . • 

4mm ^ 

x 4 

... 5 .9 (n ) 

g 15 (X) 


Xl - 1 

... 5 .9 ( p) 


x lmax 

916* 


*2 -1 

... 5.9(g) 


x 2max 

g 12 (X) 


X 3 

___iL „ i 

... 


x 3max 

g 18® 


X 4 _ _ 1 
x 4max 

... 5 .9 ( s) 
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5.10 VOLUME OF STEEL PER METRE LENGTH OF THE WA LL 
5.10.1 Stem 

Main steel required in the stem portion , 


M, 


si 


‘ I sV T*] 




1n 4 2 

r . 10 cm 


... 5.10(a) 


where a is the permissible stress in tension for steel. 

S"t 

As per the clause 9.2.1 of IS : 456 (1964), the 
minimum steel should be 0.15% of the gross-sectional area of 
concrete. So 


A sl 


0.15 

100 


x. 


10* cm' 


... 5.10(b) 


The steel required is the greater of the two 

values given by Equations 5.10(a) and 5.10(b). 

As 12 mm diameter bars are used (sectional area 

1.13 cm 2 ), the spacing required is = * (1.13) cm which 

sl 

will be suitably rounded up to a value, say a nd A gl now 

represents the steel actually provided. 

Assume that the main steel in the stem is suitably 

curtailed such that, for every four bars, two alternate bars 
are stepped at a height § V one bar at a height 3 ^ and 
the remaining bar taken for the full height h,, where 
h 1 = h - x^. 

The volume of main steel in the stem can be 


written approximately as 
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. The volume of steel in the toe and the heel portions 
can be similarly worked out and is only indicated in the 
following articles. 


5.10.2 Toe Slab 


Main steel, A g3 


°st j( *2 


cj 


in 4 2 
10 cm 


... 5.16(a) 


or 


0.15 

100 


x. 


10* cm" 


. . . 5.16(b) 


whichever is greater. 

The spacing of 12 mm diameter bars provided (sectional 

area 1.13 cm 2 ) is = (1.13), and after suitably rounding it 

A s3 

to a value y 3 , the actual area provided is called A g3 

0.15 1n 2 2 

Distribution steel, A g4 - j*oo~ * x 2 * 1U 011 

... 5.17 


The spacing of 10 mm diameter bars provided is 

= 122 (0.785) and after suitably rounding it to a value y 4 , 
A s4 

the actual area provided is called A g4 * 


Volume of main steel, v sm 2 A s3 * ^ 

2 

Volume of distribution steel, V gd2 = A s4 ‘ 10 ^ 


... 5.18 
. . . 5.19 


Therefore, 


= v 


sm2 


+ V 


sd2 


Volume of steel in toe slab, V gb 


5.20 
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5.10.3 Heel slab 


Main steel, a _ 
s5 


M, 


°st j (x ; 


c b )_ ' 


in 4 2 

10 cm 


... 5.21(a) 


__ 0.15 ^2 2 

° r " Too" * x 2 • 1° ccn ... 5.21(b) 

whichever is greater. 

The spacing of 16 mm diameter bars provided (sectional 

area 2.01 cm ) is = (2.01) and after suitably rounding,, 

s 5 

it to a value y,-, the actual area provided is called A 

J Sb 

Distribution steel, A g6 is the same as A provided 
at spacing . 


Volume 

of 

main steel. 

V 0 = 

sm3 

A s5 

. lo 2 

3 

cm 

... 5.22 

Volume 

of 

distribution 

steel, 

V sd3 

1! 

> 

CO 

0^ 

.lo 2 cm 3 

... 5.23 

Therefore 







Volume 

of 

steel in heel slab. 

V sh = 

: V O 

sm3 

+ V sd3 

... 5.24 


3 

If w is the density of steel in kg/cm , the total weight 
s t 

I 

of steel may be obtained after increasing it by 10%, to 
account for bar bending, overlapping, hooks etc. Thus 


W st = 1 ‘ 1(V ss + V st + V sh } W st •** 5 * 25 

Also, volume of concrete per running metre of the 
wall may be obtained on the basis of gross sectional area 
from Equation 5.6 as 

t r 2 A ... 5 .26 
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5.12.3 Properties of Concrete and steel 

The concrete mix used for the retaining wall is of 

grade M'*150 and the reinforcement is of plain mild steel oars. 

Referring to Tables VI and VII of IS ; 456 (1964), the 

permissible values of stresses for this concrete grade and 

steel are: a cb = 50 kg'/ era 2 , t c = 5 kg/cm 2 and = 1400 kg/cm 2 . 

As per the Clause 6. 3. 2- (a) of the above reference, the 

modular ratio 'm' has value 19, so that the design constants 

are: j = 0,865 and Q = 8.737 kg/cm 2 . The unit weight of 

3 

reinforced concrete is taken as w c = 2400 kg/m . 

These values, appropriate units 3 have been used 
in Equations 5.9, 5.10(a), 5.16(a) and 5.21(a). 

5.12.4 steel Cover 

The cover provided for the main reinforcement in 

the stem is c =4 cm and that in the base slab, cv = 5 cm. 
s 

5.12.5 Factors of Safety 

The factors of safety against overturning (PS Q ) and 
sliding (FS g ) have been taken as 1.5. 

5.12.6 Unit costs of Steel and concrete 

The material costs of steel and concrete have been 

taken as: r = Rs. 3 per kg and r c = *. 300 per m 3 respect- 
s 

ively . 
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5.12.7 Lower and Upper Bounds on the Design Variables 

As suggested by Huntington (1967), the minimum 
batter for the face should be 1 in 24. on this basis, the 
value of the lower bound x lmin is taken as (b + ~ ) • Based 
on the section requirements, as suggested by Huntington (1967), 
Fergusson (1965) and Tschebotarioff (1951), the values of the 

remaining lower bounds have been fixed as; x 2rn;Ln = 0.30 m, 

h _ , h 

*3min “ 2 and x 4 = 6 ' 

Similarly, the values of the upper bounds on the 

design variables have been decided as; x 2 max ~ 

, . h 

x_ = h and x . „ = r • 

3max 4max 4 


5.13 OPTIMIZATION PROBLEM 

Thus the present optimization problem is of const* 
rained type, in the form; 
to find 


X = 



so as to minimize the objective function 


f(X) = r s W gt + r c V c 


5 5 

4 • • 


... 5.27 


and subject to the inequality constraints 


SjOO ^ 0 


j =1, 2/ • « « / 18 • • * 5*9 
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The solution is attempted by the interior penalty 
function method which has been explained m Article 3.8. In 
this case, while obtaining the unconstrained minimum of the 
penalty function, the method suggested by Powell (1964) has 
been used, as it needs evaluation of the function values ally. 

The iterative procedure of the Powell's method can 
be stated by the following algorithm: =■ 

(i) Choose a starting point 5E ] _. Set the direction vectors 
equal to the co-ordinate unit vectors, i = l,2,...,n 
such that 

f 1 


S 1 = 




(ii) Find the minimizing step lengths X , i = 1, 2,..., 


n 


»»» 'k 

so that £(X. + X. s.) is a minimum along the direction 
1 11 


S^, and define 


_ _ * _ 

X. 1 = X. + X . s. 

l+l l ii 


In minimizing along a direction S^, the method of 
quadratic interpolation can be preferably used, as it 
needs the function values only. 

(iii) Find the integer m, 1 ^ n+1 so that the quantity 


4 = f(x, 4 ) - £<x m > 


is 


maximum. 
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(iv) Check for the convergence; if satisfied, terminate 
the process otherwise go to the next step. 

(v) calculate f 3 _ = f(2j^ - x ) 

f ! = f(3c 1 ) 


and , 


^5 ^ (-^v. , i ) • 


(vi) If either and/or 


+ f 3 ) (f ^ 


use the old set of directions Sw S 9 , ... , S for the 

next cycle of minimi zations, and use X n+1 for the next 

X x and go to step (ii). Otherwise go to the next step. 

(vii) Define the new search direction s = X - X and 

n+1 n 

* — * 

calculate X so that f(X + X S) is a minimum along 

the direction S. Then start the next cycle of minimi- 

* M 

zations by taking (X n+1 + 7 s) as the starting point 

V S 2 S m _ r S m+1 , S m+2 s n , s as the 

new set of search directions and go back to Step (ii). 


5.14 COMPUTER PROGRAMME' 

The computer programme prepared under the name OPTl 
consists of the following subroutines performing functions 
as indicated . 

(i) MAIN — Reads the data, performs various functions 

in the optimization procedure using ot^her 
subroutines and finally gives the optimal 
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dimensions of the section , the minimum 
value of the objective function and the 
corresponding values of constraints in the 
problem. 

(ii) SUMTl Obtains the constrained minimum of a 

function using the sequential unconstrained 
minimization technique (SIJMT). 

(iii) PENLTY — Lists the given set of constraints and 

generates therefrom the objective and the 
penalty functions at a design point X. I 

(iv) POWELL Finds the unconstrained minimum of a 

l 

function by the conjugate directions method. 

(v) UNIDIR «— » Minimizes a function in one direction by , 

the method of quadratic interpolation. ; 

The results for all the cases above have been 
discussed and the conclusions drawn, in chapter 6. [ 

| 

I 

« O - f 

r 

l 



6» DISCUSSION OF RESULTS AND CONCLUSIONS 


G.l COMPUTATIONAL FEATURES OF THE PROGRAMMES 

The computational work is carried out on IBM 7044 
and DEC 1090 computer systems at the Indian Institute of 
Technology, Kanpur and DEC 1077 system at the Tata institute 
of Fundamental Research, Bombay. 

By using the optimization programmes OPTl, 0PT2, 
0PT3 and 0PT4 , the optimal solution in every case has been 
obtained using different starting points. However, during 
their execution, these programmes exhibit certain features. 

6.1.1 Programmes 0PT2, 0PT3 and 0PT4 

The starting design point for these programmes is 
chosen from the available range of rolled steel sections 
given in SP : 6(1) (1964). Each of such points represents 
a practically feasible solution of the problem. 

The characteristic feature of these programmes is 
that, the optimum value of the objective function is 
decreased as the starting point successively approaches the 
optimal design point. During this process, if one starts 
from a design point corresponding to a very heavy section 
of equal flange widths and thicknesses, the optimal section, 
thus obtained, has unequal flange widths and thicknesses. 
This is obvious because the programmes 0PT2, 0PT3 and 0PT4 
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handle each of the optimization problems as a general 
Six “design-variable problem. However/ as the starting point 
approaches the optimum one/ it is observed that while the 
optimal solution is arrived at, the penalty function is 
minimized with an increase in the value of the objective 
function. As reported by Fox (1971)/ such type of behaviour 
is not uncommon with the present form of the penalty function 
chosen for the constrained minimization problems. While this 
is unfortunate, no other course is usually possible. 

The requisite logical provision incorporated in 
the programme switches the final iterative solution back to 
the starting point, which is a section with equal flange 
widths and thicknesses. By successively decreasing the 
sectional dimensions, the starting point can be taken as 
close to the optimum solution as possible. The optimal 
solution, thus obtained, is invariably a section with equal 
flange widths and thicknesses. 

A typical case in the programme OPT2, which 
illustrates all these foregoing features, is given m the 
form of Table 6.1. It can be seen that the optimal solution 
is obtained for the starting point no. 8 which is nearest 
to the optimal point. This is further explained by the 
fact that the next point no. 9 j'ust falls into the infea- 
sible region. 



variation of Optimal Solution With Different Starting Points 
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active constraint is on the normal stress for a span upto 

4.0 m, while for higher spans it is on the deflection. The 
optimum weight increase is from 21.29 to 50.00 kg/m for a 
loading intensity of 3000 kg/m, the active constraint in each 
case being that on the normal stress. 

For a cantilever, the optimum weight range is from 
8.30-18.88 kg/m and 11.34-25.44 kg/m for the loading inten- 
sities of 1000 and 2000 kg/m respectively, the span ranging 
from 1.0 to 2.0 m. In each case, the constraint on the 
deflection remains active. For a loading intensity of 3000 
kg/m, the constraint on the normal stress is seen to be active 
except for a span of 2.0 m where the constraints on the normal 
stress and the deflection are both active (Table 6.7). The 
optimum weight varies from 13.45 to 30.52 kg/m in this case. 

6. 2. 1.2 Linearly varying load 

Tables 6.8 through 6.19 show the results in this 

case . 

In the case of fixed supports, for a span range of 

3.0 to 6.0 m, the optimum weight increases from 8.16 to 19.28 
kg/m and 15.52 to 36.07 kg/m as the loading intensity ranges 
change from 0-1000 kg/m to 2000-3000 kg/m respectively, in 
each case the active constraint happens to be the normal 
stress in the beam, excepting when the span is 3.0 m and the 
loading intensity 1000 kg/m. in this case, the prescribed 
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lower bounds on the top and the bottom flange widths act as 
the active constraints, a s indicated in Table 6.8. 

When the supports are simple, for a span increase 
from 3.0 to 6.0 m, the variation in the optimum weight is 
from 9.66 to 21.74 kg/m and 19.17 to 44.56 kg/m, as the loading 
intensity range changes from 0=1000 kg/m to 2000=3000 kg/m 
respectively. For loading intensities of 0=1000 kg/m, 0=2000 
kg/m, 0-3000 kg/m and 1000=2000 kg/m, the constraint on the 
deflection is seen to bo active, excepting for a span of 6.0 
m (loading intensity being 0=1000 kg/m) and a span of 5.0 m 
(loading intensity being 0=2000 kg/m) where the constraints 
on the stress and deflection are both active, as indicated in 
Tables 6.8 and 6.10 respectively. When the loading intensity 
is 1000-3000 kg/m, the constraint on the stress is active 
upto a span of 4.0 m after which the constraint on the 
deflection is active (Table 6.16). Lastly, for the loading 
intensity of 2000-3000 kg/m, the constraint on the stress is 
found to be active for all spans, except when it is 6.0 m 

A 

in which case the constraints on the stress and on the 
deflection are both found to bo active (Table 6.18). 

In the case of a cantilever, for a loading 
intensity of 1000=0 kg/m, the lower bounds on the top and 
bottom flange widths are seen to be the active constraints 
for spans of 1.0 and 1.5 m (Table 6.9). This pair of active 

constraints is also observed when the span is 1.0 m and 

2000=0 kg/m (Table 6.11) and 


loading intensities are 
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3000-0 kg/m (Table 6.13) respectively. The variation of the 
optimum weight is from G .16 to 11.66 kg/m and 11.97 to 27,12 
kg/m as the loading intensity changes from 1000-0 kg/m to 
3000-2000 kg/m . FOr a span of 1.5 m and loading intensity 
3000—2000 kg/m, the constraints on the stress and the 
deflection are active as seen from Table 6.19. in the 
remaining cases, except those mentioned above, the constraint 
on the deflection is active. 

6. 2. 1.3 Nodal loads 

The results are presented m Tables 6.20 through 

6 . 22 . 

For fixed supports, the variation in the optimum 
weight is from 8.16 to 21.74 kg/m and from 14.39 to 33.87 
kg/m for a span increase of 3.0 to 6.0 m, when ‘the equivalent 
intensity of the nodal loads changes from 1000 to 3000 kg/m. 

The corresponding figures for the optimum weight 
variation in the case of simple supports are 11.97 to 20*84 
kg/m and 21,29 to 50.00 kg/m respectively. 

For a span of 3.0 m under the fixed support 
conditions, when the equivalent loading intensity is 1000 
kg/m, the optimal solution is governed by the constraints 
on the lower bounds on the top and the bottom flange widths. 
On the other hand, the constraint on the deflection is active 
when the span is 4.0 m under the simple support conditions. 

In all the rest of the cases, the constraint on the normal 


stress is active. 
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6. 2. 1.4 Uniformly distributed load with a lateral 
eccentricity 

Tables 6.23 through 6.37 give the results for this 

case . 

In the case of fixed supports , for a loading 
intensity of 1000 kg/m at a lateral eccentricity of 0.1 bf^ 

(bf being the top flange width of the optimal section for 
the corresponding case of the in-plane loading)/ the optimum 
weight increases from 11.05 to 28.84 kg/m as the span changes 
from 3.0 to 6.0 m. When the lateral eccentricity is increasec 
to values 0.2 bf 1 and 0.3 bf^ the corresponding ranges of 
the optimum weight are 12.57-31.82 kg/m and 14.11-34.76 kg/m 
respectively. 

For these values of the lateral eccentricity/ the 
corresponding increases in the ranges of the optimum weight 
are 17.13-39.60 kg/m/ 19.94-47.19 kg/m, 22.97-54.39 kg/m and 
22.05-51.49 kg/m/ 26.06-62.61 kg/m and 29.04-71.67 kg/m when 
the loading intensity assumes values of 2000 kg/m and 30 
kg/m respectively. 

in all these cases, the constraint on the normal 

stress is seen to be active. 

However, In the case of simple supports, when the 

loading intensity is 1000 kg/m at a lateral eccentricity of 
0.1 b fl , the increase in the optimum weight is from 12.84 
to 31.19 kg/m as the span changes frcm 3.0 to 6.0 ». For 
spans upto 4.0 m, the constraint cn the deflection is active 
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while for higher span values , the constraint on the normal 
stress takes over as the active one, as can be seen from 
Table 6.23. Both the norncl stress and the deflection 
constraints are ' imultancously active for spans ranging from 
4.0 tc 4-5 m approximately. 

To r the lateral eccentricity values of 0.2 bf-^ a nd 
0.3 b£ ^ / the corresponding increases m the optimum weight 
are from 13.56 to 33.70 kg/m and from 14.80 to 36.07 kg/m 
respectively. 

As the loading intensity is increased to 2000 and 
3000 kg/m, the respective ranges for these three values of 
the eccentricity are 19.68-46.13 kg/m, 22.34-52.25 kg/m, 

25 ,01-58.87 kg/m and 25.64-59.35 kg/m, 29.85-68.67 kg/m and 
33 , 50-79 . 29 kg/m. In all these cases, the constraint on the 


normal stress is the active one. 

In the case of a cantilever, when the loading 

intensity is 1000 kg/m at an eccentricity of 0.1 To£ y the 

optimum weight varies from 8.31 to 19.17 kg/m as the sp 

_ £ i n +-n 9 0m This variation for the loading 

changes from 1.0 to 2.o m. 

• -X. * n o VvF !S from 9.36 to 22.34 kg/m. 

eccentricity of 0.2 bf^ is 

For the increase in loading to the values of 2000 
kg/m and 3000 kg/m, the corresponding variation in the 
optimum weight for these two values of the lateral eccent 
ricity are 13.11-29.55 kg/m, 15.11-34.21 kg/m and 16.76 
38.55 kg/m and 19.40-45.40 kg/m respectively. 
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In all these cases, the constraint on the normal 
stress is found to be active. 

6.^.2 Optimization of Thin “Walled Sections under stability 
Constraints 

The results of this optimization problem are shown 
graphically in Figure 6.1. 

It may be noticed that for all types of boundary 
conditions, the optimum weight of the column increases non-- 
linearly with the increase in the height of the column. For 
fixed end conditions, the optimum weight increases from 3.16 
to 24.16 kg/m as the column height increases from 3.0 to 
6.0 m. For a column height of 3.0 m, the lower bounds on 
the top and the bottom flange widths act as the active cons- 
traints. For the remaining column heights, the constraint 
on the maximum slenderness ratio of the column is active. 

When the ends are simply supported, the variation 
in the optimum weight is from 11.50 to 26.50 kg/m as the 
column height increases from 3.0 to 4.5 m. The corresponding 
figures for a cantilever column are 24.40 to 32.83 kg/m for 
a height increase from 3.0 to 4.0 m. In all these cases, 
the constraint on the maximum slenderness ratio of the 
column is seen to be active. 

6.2.3 Optimization of Thin-walled sections under Dynamic 
Constraints 

The results for this optimization problem are 
presented graphically in Figures 6.2 through 6.4. 
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For fixed supports, the optimum weight increases 
from 8.16 to 10.88 kg/m, a s the span changes from 3.0 to 
6.0 m and the upper and lower bounds on the first value of 
transverse vibrational frequency have a broad range . Upto 
the span of 5.0 m, the active constraints are the lower 
bounds on the top and the bottom flange widths of the 
section. For the remaining spans, the lower bound on the 
first value of natural frequency is seen to be active. 

The corresponding variation m the optimum weight 
is from 8.16 to 33.50 kg/m when the bounds on the first 
frequency have a narrow range. In this case, the lower 
bounds on flange widths govern the optimal design for a 
span upto 3.50 m, after which the lower bound on the first 
frequency behaves as the active constraint. 

When the supports are fixed, for the increase in 
the span from 3.0 to 6,0 m, the optimum weight increases 
from 8.16 to 19.03 kg/m and from 9.50 to 41.28 kg/m when the 
range of frequency bounds is broad and narrow respectively. 
In the case of a broad range, the lower bounds on flange 
widths govern the optimal design for spans upto 4,0 m. 
Otherwise, the active constraint is always the lower oound 

on the first natural frequency. 

in the case of a cantilever, for a variation of 
span from 1.0 to 2.0 m, the optimum weight enhances from 
8.16 to 37.84 kg/m and from 27.55 to 48.63 kg/m when the 
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bounds on the first natural frequency have a broad and a 
narrow range respectively. In the case of a broad range, 
when the span is 1.0 m, the lower bounds on the two flange 
widths form the active constraints while in the remaining 
cases the lower bound on the first natural frequency always 
governs the optimal design. 

6.2.4 Optimization of a cantilever Retaining Wall 

The results for this optimization problem are 
presented graphically in Figures 6.5 through 6.7. 

As can be naturally expected, for any type of soil, 
the optimal cost rises with the increase in height and top 
width. However, the salient features of results are discussed 
according to the type of soil used for the backfill. 

6. 2. 4.1 Gravel 

For a top width of 0.20 m and the wall height upto 
3.5 m, the optimal cost per running metre of the wall varies 
from Rs. 454.44 to Rs. 578,13. The constraint governing the 
optimal design is the lower bound on the stem root thickness. 
However, for a height of 4.0 m, the active constraint changes 
to the lower bound on the base slab thickness the corres~ 
ponding optimum cost being Rs. 696.16. A further increase in 
the wall height from 4.5 to 6.0 m enhances the optimum cost 
from Rs. 866.94 to Rs. 1530.29 and the active constraint is 
the lower bound on the base slab width. 
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6. 2. 4. 5 Clay 

For a top width of 0.20 m, the increase in the 
optimum cost is from Rs. 518.33 to Rs. 2012.39 as the wall 
height increases from 3.0 to 6.0 m. For a height upto 3.50 m, 
the lower bound on the stem root thickness governs the design. 
When it is above 4.0 m, the active constraint is either the 
lower bound on stem root thickness to resist bending moment 
or the no tension condition in the foundation. Similar trend 
of the active constraints is observed when the top width has 
values of 0.25 and 0.30 m. In these cases, as the height 
increases from 3.0 to 6.0 m, the optimum cost increase is 
from Rs. 562.49-2029.66 and Rs. 601.92-207 3.37 respectively. 

6.3 CONCLUSIONS 

The conclusions are drawn separately for every 
optimization problem handled so far. 

6.3.1 Optimization of Thin -Walled Sections under Static 
Constraints 

This problem has been solved for four different 
cases of the loading and as such it is found logical to draw 
conclusions for each loading case and then deduce overall 
conclusions . 

6. 3. 1.1 Uniformly distributed load 

As the loading intensity varies frcm 1000 to 3000 
kg/m, the optimum weight increase in case of fixed support 
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conditions is 86.40% for a span of 3.0 m, but it reduces to 
55.02% for a span of 6.0 m. The corresponding values in 
respect of simple supports and a cantilever are 65.80%- 
68.35% and 62 .05% -61 .65%, for a span variation of 3. 0-6.0 m 
and 1. 0-2.0 m respectively. 

The effect of boundary conditions may also be noted 
with interest. For a span of 3.0 m, as the supports change 
from fixed to simple, the increase in optimum weight falls 
from 41.87% to 26.20% when the loading intensity increases 
from 1000 to 3000 kg/m. However, for a span of 6.0 m the 
increase in optimum weight has a rise from 18.75% to 28.96%. 

Hence, for smaller load intensities, fixed beams 
are very much economical than simple beams, provided the 
spans are small. For higher intensities also, significant 
economy can be achieved by using fixed beams for small or 
large spans . 

It is worthwhile to note the variation in the type 
of active constraint in the optimal design. For fixed 
supports, the constraint on the normal stress is active under 
all spans and loading intensities. In case of simple supports, 
the deflection constraint governs for small loading intensities. 
For medium loading intensities, the stress constraint is 
active for smaller spans, and the deflection constraint for 
larger spans. When the loading intensities are large, the 
stress constraint dictates the optimal design for all spans. 
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In case of cantilever beams , for all spans, the deflection 
constraint is governing for small and medium loading inten- 
sities but for larger intensities, the stress constraint 
dominates . 


6. 3. 1.2 Linearly varying load 

Generally speaking, for this type of loading, the 
constraint on the normal stress remains active for all spans 
and loading intensity ranges, when the supports are fixed, 
in the case of simple supports, the deflection constraint 
governs the optimal design for all spans when the ranges of 
loading intensities are small or medium. If these are large, 
the stress constraint is active for small spans while the 
deflection constraint dominates m case of large spans. In 
the case of a cantilever beam, the deflection constraint is 
active under all spans and loading intensity ranges. 

In the case of fixed supports, for an increase m 
loading intensity range from 0-1000 kg/m (average value 
500 kg/m) to 2000=3000 kg/m (average value 2000 kg/m), the 
optimum weight enhances by 90.19% and 87.09% for spans of 
3.0 m and 6.0 m respectively, corresponding values in case 
of simple supports and a cantilever are 98.44% =104 .96% and 
46.69% -132.59%, when the span variation in the respective 
case is 3. 0=6.0 m and 1. 0=2.0 m. 

This clearly shows that with the increase in 
loading intensity, cantilever beams in respect of the 
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optimum weight are very much sensitive to the increase in 
span, than the fixed or simply supported beams. Hence, for 
a linearly varying load, it is desirable to cut down the 
spans of cantilever beams to reduce the structural cost, 
ultimately. 

One may also investigate the possibility for a 
choice between simple and fixed beams for equal spans, in 
the former type of beam, for a range of loading intensity from 
0=1000 kg/m, the increase in the optimum weight is 18.38%, 
when the supports are made simple for a span of 3.0 m. 

Similar increase for a span of 6.0 m is 12.76%. corresponding 
values for a loading intensity range of 2000-3000 kg/m are 
23.52% and 23.54% respectively. 

This clearly suggests that for larger load 
intensity ranges, fixed beams are economical than simple 
beams for smaller or larger spans. For smaller load inten- 
sity ranges, fixed beams with smaller spans may be somewhat 
economical; but for large spans the economy achieved is only 
marginal. 

6. 3. 1.3 Nodal loads 

The constraint on the normal stress is always 
active excepting the case of simple supports under small 
spans and loads, when the deflection constraint prevails. 

When the nodal loads have equivalent intensity of 
1000 kg/m, for spans of 3.0 and 6.0 m, the increase in the 
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optimum weight by changing the support conditions from fixed 
to simple, is 46.69% and 32.66% respectively, corresponding 
figures are 47.95% and 47.62% respectively, when the equi- 
valent intensity has a value of 3000 kg/m. 

In practice, nodal loads can assume the form of 
reactions transferred on a main girder, from the cross-girder 
systems. When the load magnitudes are higher, the fixed 
beams are expected to be economical under all spans. Even 
for smaller load magnitudes, the economy achieved is 
comparatively more for smaller spans. 

6. 3, 1.4 Uniformly distributed load with a lateral 

eccentricity 

Under all loading intensities with different 
eccentricities, the constraint on the normal stress is seen 
to govern the optimal design, for almost all spans and 
boundary conditions; except, in the case of simple beams 
carrying loads of smaller intensities and eccentricities , 
the deflection constraint is active, when the spans are 
smaller. 

When the loading intensity is lOOO kg/m, for a 
smaller eccentricity of 0.1 bf^, the increase in the optimum 
weight due to the change-over of suppprt conditions from 
simple to fixed, is 16.20% and 8.15% respectively as the 
span varies from 3.0 to 6.0 m. For a^ large eccentricity of 
0.3 bf ^ , these increases are 4.89% and 3.77% respectively. 



Hence, for smaller loading intensities, fixed beams will be 
economical than simple beams, if the loading eccentricity 
and the span have smaller values. 

For a larger loading intensity of 3000 kg/m, with 
the increase in span from 3.0 to 6.0 m, corresponding vari- 
ation in the optimum weight for the loading eccentricity value 
of 0.1 bf^ is 16 . 28% -=-15 . 27%; for 0.3 bf^, it is 15.36%-10.63%. 
Hence, for larger loading intensities having smaller eccen- 
tricities, fixed beams offer more economy than simple beams, 
for all spans; for larger eccentricities, significant economy 
can be expected for smaller spans. 

In the case of cantilever beams, for all loading 
intensities and spans, the increase in the optimum weight is 
rather insignificant, being of the order of 15% when the 
loading eccentricity is changed from 0,1-0. 2 bf ^ . Therefore, 
if the situation demands the necessity of larger eccentric ite s , 
these are permissible, without greatly affecting the economy. 

As the last investigation m this type of loading, 
values of the optimum weight are compared with those obtained 
for the in-plane loading under identical spans and loading 
intensities. From this information, one can have an idea 
of the significance of lateral eccentricities in the 
structural design, under various boundary conditions. 

For fixed beams, under a small loading intensity 
with an eccentricity of 0.1 bf^, the .increase in the 
optimum weight is around 15-20%; for a large eccentricity 
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of 0.3 bf ^ , the increase is 40-50%. For large loadings, 
such increases for small and large eccentricities are about 
30% and 75% respectively. Hence barring,, perhaps , the small 
intensity loadings with smaller eccentricities, the laterally 
eccentric loads should be scrupulously considered in the 
design of fixed beams. 

In simple beams, such figures are around 5% and 20% 
for small loads and, 15% and 60% for large loads, when the 
eccentricities are small and large respectively. Hence, in 
the design of simple beams, only large eccentricities become 
significant in the design, for all intensities of loadings. 

In cantilevers, for small and large eccentricities, 
similar values are around 5% and 15%, and 25% and 50% 
respectively. This signifies that large loadings, even at a 
small eccentricity, need special consideration m the design 
of cantilevers. 

6. 3. 1.5 General remarks on optimal sections under static 

constraints 

Generally speaking, for any type of loading described 
above, the fixed beams are economical for all spans, parti- 
cularly under large loading intensities; for small intensities, 
they are economical for small spans only. 

In the design of gridworks, the joints need not 
always be designed as rigid, thus imposing the fixed boundary 
conditions on the grid members . In some loading cases , even 
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simple boundary conditions could lead to an optimal design 
of that grid member whose joints are designed for such 
boundary conditions. 

Tables 6.2 through 6.37 provide optimal sectional 
dimensions for an I -section for various loadings and boundary 
conditions. Given a flexibility in the production/ such 
members can be rolled and used to economize thin-walled 
structures/ particularly when the component members are 
required on a large scale. If the situation does not permit 
this, the nearest available section would provide a practically 
feasible optimal structural design. 

6.3.2 Optimization of Thin-Walled sections under Stability 
Constraints 

For a column height of 3.0 m, the optimum weight 
increases by 40.93% as the support conditions change from 
fixed to simple, thus doubling the effective height; however, 
a further increase of 112.17% results as it becomes a 
cantilever column and the effective height is further doubled. 
Corresponding values for a column height of 4.0 m are 129.94% 
and 37.02% respectively. 

Hence, generally speaking, for small heights, 
columns with simple supports may not be so uneconomical 
with respect to the fixed columns but the cantilever columns 
would be highly uneconomical as compared to simple ones. 
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decreased, it is observed that in these soil types, the lower 
bound on the stem root thickness governs the optimal design 
for small heights, while for larger heights the constraint on 
resistance against sliding is active. This is particularly 
true in case of smaller top widths, which ultimately provide 
smaller base widths . 

The situation is slightly different in the case of 
silt and clay which have poor bearing capacities but possess 
cohesive strength that adds to their resistance against 
sliding. For smaller heights, the lower bound on stem root 
thickness governs the optimum design while for larger heights 
the constraint on the safety against sliding is active, this 
being true for all values of top widths, comparatively, silt, 
having higher bearing capacity but lower cohesion, is expected 
to lead to the optimal design that is economical upto a 
certain height, than that obtained with clay as the backfill. 
For heights greater than this, when the constraint on the 
safety against sliding is bound to govern, optimal design in 
case of clay, which possesses more cohesion, is bound to be 
cheaper. 

Referring to Figures 6,5 through 6,7, it is 
deducible that in case of silty backfills, for different 
top widths, the optimal designs of cantilever retaining walls 
are economical upto a height of 4.5-4.75 m. After this, 
height, clayey soils yield a more economical design, although 
the saving may be marginal, say of the order of 10%. 
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6.3.5 Computational Algorithms 

The computer programme OPTl deals with the optimi- 
zation of a cantilever retaining wall, a structure belonging 
to the thick-walled category. The optimal designs of thin- 
walled rolled steel sections under the Static, the Stability 
and the Dynamic constraints are handled by programmes 0PT2, 
0PT3 and 0PT4 respectively. Both these constrained optimi- 
zation problems are approached using the Interior Penalty 
Function method; but in each problem, different methods are 
adopted for the unconstrained minimization involved in the 
procedure, and these present different features. 

In programmes 0PT2, 0PT3 and 0PT4 , the Davidon- 
F let cher -Powell algorithm is applied which makes use of the 
gradients of the objective function. Although numerically 
stable even under eccentric functions, in the present problems 
it is seen to present difficulties while evaluating gradients 
at points near the constraints. As can be seen from the 
optimal solutions presented in Tables 6.2 through 6.37, all 
optimum design points are bound points. Hence, while using 
the search technique, as the design point approaches the 
active constraint, the penalty function tries to ‘blow off', 
which subsequently presents difficulties in evaluating 
gradients, thus halting the minimization process. The 
requisite logical provision incorporated in the programmes, 
however, takes care of such a situation. 
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The programme OPTl employs the Powell's method, 
that incorporates modifications for handling non -quadratic 
objective functions. This method needs function evaluations 
only and hence in the vicinity of the optimal solution, 
normally a bound one, no numerical instability is experienced. 

Regarding the computational accuracy, both algori- 
thms are seen to be equally satisfactory. 

6.4 SUGGESTIONS FOR FURTHER WORK 

The present thesis is just a beginning of the 
problem of optimization of thin- and thick -walled structures 
m Civil Engineering practice. 

The problem can be expanded to include various 
other types of loadings, analysed using different techniques 
and optimized with due consideration to many other varieties 
of constraints . Thus , the next steps suggested to enhance 
the scope of the present problem are as follows; 

(i) In the optimization of thin -walled sections, the 

stability and/or dynamic constraints can be included 
along with the static constraints, so that optimal 
designs could be compared for different combinations 
of constraints . 

(ii) As the reactions from the cross-girders can, in 

general, be a force and a twisting and/or bending 
moment, such concentrated bending and twisting 
moments can be considered among the nodal loads. 
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increase falls to an insignificant 8.27% for a wall height 
of 6.0 m. Similar trend is observed in cases of other soil 
types considered namely, dense sand ( 22 .4% -4 .63%) , loose 
sand (15.89%-1.93%) , silt (17 .65%~1 .67%) and clay (16.12%= 
3.03%). This clearly shows that for any type of soil 
retained, the effect of top width on the optimum cost of the 
cantilever retaining wall is significant at smaller heights 
only; with larger heights, the effect dies out rapidly. 

6, 3. 4. 2 Significance of active constraints 

Amongst the five different classes of soil considered, 
gravel, dense sand and loose sand belong to the cohesionless 
type while the remaining two, namely silt and clay are 
cohesive. In the cohesionless variety, gravel possesses high 
values of the bearing capacity and the angle of internal 
friction, while these values, get successively reduced for 
dense sand and loose sand. Naturally, in case of gravel, 
even at higher wall heights, sufficient frictional resistance 
against sliding is assured by its high value of the angle of 
internal friction. Hence, the lower bound on the stem root 
thickness governs the design. However, m case of dense sand 
and loose sand, the decreasing values of bearing capacity, 
demand higher base widths and subsequently higher wall 
dimensions. This obviously results in successive increase 
in the cost for dense sand and loose sand, it being maximum 
in the latter case. Again, the frictional resistance getting 



optimal sectional Dimensions 


TO 

0 


P 

G 


P 

P 

a 

•Hi 

U 

P 

a 

0 

a 

o 

M 


fO 

H 

0 

P 

0 

P 


CO 

P 

O 


2 rfl 

0 CD 0 X 

•P (1) C fD +J 

q > o n c 

^ *H CJ P *H 


8 

■H 

Si 

- q P- 

i 

j 


0 

H P 
0 £ 


s 

... c cn 

&-ri x. 


IP 

ft 

G 


3 

P 

p S 

to 


4 


0 


0 

G 


P 

0 

CM 

& 

•H 

MH 

E 

0 

P J 

■H 

G 


to 

0 



E 

CM 


*H 

m 


TD 

P 


P 




0 

I 

ij 


E 

P 


•H 

P 


P 

P J 


ft 

i 


o 

! 

' 


P j 

1 

] 


1 P 


i 

- - ^ J 

! 



! , CO i 


G 

to f \ 


P 0 C 0 P 1 


u t 

o u c 


j_ £*u z 

O 4 -> H ( 

P 

1 

■ 

. « ! 


P E "H P 
P* 3 0 P 


\ 
£ P 


0 

p 

H 

O 

& 

Pa 


o. A A J5 ii £L 


6 

e 

G 

-H 


2 


vO 

o 


* lD 

LO 

! - 

N 

Lf 

* rH 

rH 

CO 

rH 



tH 

tn 



tn 

rH 

O 

CD 


rH 


• 

9 

• 

• 

CM 

!> 

oo 

cn 

\ — 1 


CM 

CM 

tO 

CM 

I> 

j CM 

% 

» 

• 

• 

^ ; 

if) 

if) 

i 

vO 

O 

. 

o I 

0 

O 

• 

• 

l * i 

• 

m 

1 P | 

| P ! 


m 

! CO 1 

1 CM 

if) 

rH 

: H i 

; w 

(M 

to 

° 


O 

9 

• i 

• 

• 

tO 

o S 

00 

00 

o 


o 


• 

• 

• 

• 

r- 

CM 

o 

r- 


cn 

o 

CM 



rH 

rH 

vO 

o 

sh 

O 

• j 

* 

» 1 

9 

vD 

l> 

CO 

CO 

O [ 

p 

! 

O 



• 

• \ 

9 


CM 

o 

t> 

r- 1 

o\ 

o 

CM 

1 


tH ¥ 

rH 

«n | 

CO 

CO ! 

CO 

-fc rH j 

r~i 

rH 

rH 

^ i 
i 

& 


tn 

tn | 

i 

rH } 

CM | 

i rH 

O 

Ch l 

ID 

; o 

• 

* l 

i * j 

i • 

o\ j 

CM | 

CO 1 

i LO 


rH 5 

1 ^ i 

[ CM 

1 — i j 

i vo i 

! CM | 

1 o\ 

* 

| * : 

* 

* 



| LO 

1 ^ 


o 

to 

ro 


O 

o 

cn 


ro 


G 

3 



rH 


rH 

CM 

P 

to 

0 





( 

*H 


G 1 



vO 

rH 

to 1 

U 

TO 

0 

CM 

* 

* 

• 

* J 

P 

0 

•H 

P 

to 


i> 

CO 

0 °N 

X 

0 

P 




* 

•H 

•H 

G 






T5 E 

Ph 

0 


o 

CM 

o 

o 

\ 


E 

CM 


4 

• 

9 

>i D» 


H 

P 

to 

t> 

vO 

vjO 

e ^ 


TO 

rQ 

if) 

r- 

CT> 

O 








p 

M O 


rH 






0 o 


rd 



tO 

P 

tO 

P o 


E ! 

i ^ ! 

• 

• 

l 1 

• 

♦ \ 

•rH t — i 

j 

-H , 


! 

I tO 

o 

co ; 


| 

P ' 

rs 

1 — it. 1 



! i 

s \ 

*« 

J 

| o' 

) 1 

; o 

| tM ; 

1 o | 

o 


i 

i 

rH 

• 


1 « 

« 

OH p 


j 

P 

to 

i> ! 

1 tO 

tO 

tn *h 



P 

LO i 

i l> ! 

Ch ! 

o 

C 0 


j 


| 

1 ! 

! 

1 — t 

*H C 

t 



j 

j i 


3 

TO 0 

. c 

g 



i ^ ! 


j 

rd P 

nJ 



o 

! O i 

o 

o 

O G 

i ft C 


• 

_ ! 

1 * 1 

• 

I 

P H 

) CO 

•H 


CO \ 


lO 

} vo 


I 

0 

P g 

0 0 

£ £ 


G 

•H 

to 

if) 

0 

U 

P 

0 

P 

m 


0 

£ 

p 

c 

*H 

s 

*H 

P 

a 

0 


W P 
0 0 
G T5 


E 

•H 

3 

£ 

* 


19S 


**Maximum 



199 


Table 6 . 3 

Optimal Sectional Dimensions 

i 

Loading? uniformly distributed 

intensity: 1QOO kg/m ; Lateral Eccentricity: 


Span Y 
in m ! 



Cantilever 


Optimal dimensions in mm 



Optimum 
weight 
in kg/m 


1.0 51*0 6.4 51.0 6.4 101.0 4.0 8.30 


1.5 78.8 6.7 78.8 6.7 144.0 4.7 13.65 


2.0 96.8 7.2 96.8 7.2 192.0 5.3 18.88 


**Maximum deflection in the beam. 


nil 


Active 
constra - 
int/ 


®15 


SlS 


9 15 
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Table 6.5 


Optimal sectional Dimensions 


Loading; uniformly distributed 


Intensity: 2000 kg/m ; Lateral Eccentricity: 


Span 
in m 




Cantilever 


Optimal dimensions in mm 




Optimum 
weight 
in kg/m 


1.0 71.0 6.5 71.0 6.5 121.0 4.3 11.34 


1.5 92.8 7.0 92.8 7.0 182.0 5.2 17.61 


2.0 108.0 8.6 108.0 8.6 233.0 5.9 25.44 


**Maximum deflection in the beam. 


nil 


Active 

constra- 

int/s 


** 

9 15 


g 15 


g 15 
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Table 6.9 

Optimal Sectional Dimensions 

Loading; linearly varying 

Intensity; 1000-0 kg/m ; Lateral Eccentricity; nil 


Cantilever 


Span 
in m 


Optimal dimensions in mm 



Optimum 
weight 
in kg/m 


Active 
con st ra- 
in t/s 


1.0 50.0 6.4 50.0 6.4 100.0 4.0 


8.16 



1.5 50.0 6.4 50.0 6.4 100. 0 4.0 0.16 g 1# g 3 


2.0 73.0 6.5 73.0 6.5 123.0 4.4 11.66 g 15 


9-Lower bound on the top flange width. 
©9Lower bound on the bottom flange width. 
**J4aximum deflection in the beam. 
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Table 6.11 

Optimal Sectional Dimensions 


Loading: linearly varying 

Intensity: 2000-0 kg/m ; Lateral Eccentricity; nil 


Span 
in m 

T 

Cantilever 


Optimal dimensions 

in mm 


1 

1 

! Optimum 

1 

1 

J Active 


bfi 

t f i 

I ) 

i b£ 2 1 

1 1 

tf 2 j 
1 

■ 

hw J 

f 

tw 

J weight 

1 in kg/m 

1 

1. .. 

J constra- 
! int/s 

l 

1.0 

50.0 

6.4 

50.0 

6.4 

100.0 

4.0 

8.16 

A- 

" “XJ/ 

g l' g 3 

1.5 

68.0 

6.5 

68.0 

6.5 

118.0 

4.3 

10.88 

none 

2.0 

82.4 

6.8 

82. 4 

6.8 

156 .0 

4.9 

14.80 

** 

* 15 


&Lower bound on the top flange width. 


OOLower bound on the bottom flange width. 


**Maximum deflection in the beam 
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Table 6.13 

Optimal Sectional Dimensions 


Loading; linearly varying 

Intensity; 3000=0 kg/m ; Lateral Eccentricity; nil 


I 

1 

\ 


cantilever 



j--/ '-‘ 4 . X | 

in m | Optimal dimensions 

1 

in mm 


Optimum 

Active 

i i r 

j bf i ; tf i j d£ 2 

j t£ 2 

1 

hw j 

i 

^ i 

tw 

weight 
in kg/m 

constra- 

int/s 

1.0 50.0 6.4 50.0 

6.4 

100.0 

4.0 

8.16 

n\ fTv/Tv 

si/ s 3 

1.5 76.6 6.6 76.6 

6.6 

133.0 

4.5 

12.67 

** 

g 15 

2.0 91.6 7.0 91.6 

7.0 

179.0 

5.2 

17.26 

g 15 


®Lower bound on the top flange width. 
WLower bound on the bottom flange width. 
^Maximum deflection m the beam. 
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Table 6.15 

Optimal Sectional Dimensions 

Loading; linearly varying 

Intensity; 2000=1000 kg/m ; Lateral Eccentricity; nil 


Span 


Cantilever 


in m 


Optimal dimensions 

in mm 


Cotimum 

1 

; Active 


bf 1 

i 

I 1 

bf 2 i tf 2 ! 

hw ; 

i 

tw 

weight 
in kg/m 

J constra- 

J int/s 

1 

1 .0 

57.0 

6.4 

57.0 6.4 

107.0 

4.1 

9.20 

** 

915 

1.5 

32.8 

6.8 

82.8 6.8 

157.0 

4.9 

14.89 

915 

2.0 

100.0 

7.6 

100.0 7.6 

206 .0 

5.5 

20.83 

g 15 


**Maximum deflection in the beam. 
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Table 6.19 

Optimal sectional Dimensions 


Loading; linearly varying 

Intensity; 3000-2000 kg/m ; Lateral Eccentricity; nil 


T 


Cantilever 


Span 


in m 


Optimal dimensions 

in mm 

1 

1 

1 

Optimum 

1 

; Active 


b f i 

| t£ l 

i bf 2 1 

! 

tf 2 1 

hw 

1 t 

| I 

; tw ; 

i i 

! _ 

weight 
m kg/m 

[ constra- 

; int/s 

1 

1 .0 

75.0 

6.5 

75.0 

6 . 5 

125.0 

4.4 

11.97 

** 

g 15 

1.5 

95.6 

7.1 

95.6 

7.1 

189.0 

5.3 

18.52 

Z 3 ' g 15 

2.0 

116.0 

8.7 

116.0 

8.7 

241.0 

6.0 

27 .12 

s 15 


^Maximum normal stress ;n the beam. 


** Maximum deflection m the beam 
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flange width. **Maximum deflection in the beam 
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Table 6.31 

Optimal Sectional Dimensions 

Loading: uniformly distributed 

Intensity: 2000 kg/m ; Lateral Eccentricity: 0.1 bf^ 

(from Table 6.5) 


Canti lever 


open! 

in m 


Optimal dimensions 

in mm 


T" 

i 

! optimum 

i 

; Active 


b£ 1 

tf l 

bf 2 

1 

tf 2 | 

hw 

J tw 

t 

! weight 

! in kg/m 

1 

J constra- 

; int/s 

1 

1.0 

77.6 

6.7 

1 77.6 

6.7 

138.0 

4.6 

13.11 

* 

g 13 

1.5 

100.0 

7.7 

100.0 

7.7 

208.0 

5.5 

21.13 

g l 3 

o 

1 cl 

126.8 

8.7 

126.8 

8.7 

253.0 

6.1 

29.55 

g 13 

^Maximum normal 

stress 

in the 

beam . 
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Table 6.32 

Optimal Sectional Dimensions 


Loading: uniformly distributed 

Intensity: 2000 kg/m ; Lateral Eccentricity: 0.2 bf 1 

(from Table 6.5) 
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in m 

.... 
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Table 6.37 

Optimal Sectional Dimensions 

Loading: uniformly distributed 

Intensity: 3000 kg/m ; Lateral Eccentricity: 0.2 bf ^ 

(from Table 6.7) 

Span 
in m 


* 


1.5 

98.4 

7.2 

98.4 

7.2 

196.0 

5.3 

19.40 

g 13 

1.5 

133.4 

8.8 

133.4 

8.8 

264.0 

6.3 

31.34 

g l 3 

2.0 

165.0 

10.3 

165.0 

10.3 

333.0 

7.2 

45.40 

g 13 


* Maximum normal stress in the beam. 




Optimum weight in kg/m 


o - Broad range 
A - Narrow range 



Fig* 6 • 3 


4 5 

Span of beam L in metres 

Optimum weight of simple beams for 
various spans and frequency ranges 



Optimum weight in kg / m 


o - Broad range 
a - Narrow range 


' to ^ 20 

Span of beam L in metres 

Fig . 6.4 Optimum weight of cantilever 
beams for various spans and 
frequency ranges* 



Optimum cost per metre length in rupees 



heights and soils- 



Optimum cost per metre length in rupees 



Fig. 6-6 Optimum cost of reinforced concrete 
cantilever retaining walls for various 
heights and soils. 


Optimum cost per metre length in rupees 



Height h in metres 

Fig. 6.7 Optimum cost of reinforced concrete 
cantilever retaining walls for various 
heights and soils. 
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Table 5.1 

Representative Values of Soil Properties 


1 1 

No. J Type of soil ; 

i i 

i i 

i i 

i i 

i i 

w 

s 

(kg/m 3 ) 

i 0 

1 

1 

1 (degrees) 

I 

P 

per 

(kg/m 2 ) 

~ f 

1 

I C 

1 

! (kg/m 2 ) 

I 

__ _L ^ _ 

1. Gravel 

2000 

38 

55000 

0 

2* Dense sand 

1920 

33 

45000 

0 

3. Loose sand 

1880 

28 

30000 

0 

4. Silt 

1720 

22 

20000 

2000 


5 . Clay 


17 20 


18 


16000 


5000 
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6.2.1 optimization of Thin -walled Sections under static 
Constraints 

For four different loading conditions considered, 
the results are presented for different spans and boundary 
conditions in Tables 6.2 through 6.37. 


6. 2. 1.1 Uniformly distributed load 

The results are presented in Tables 6.2 through 

6.7. 

For fixed support conditions, the optimum weight 
changes from 9.05 to 25.01 kg/m as the span changes from 3.0 
to 6.0 m, with the loading intensity remaining at 1000 kg/m. 
As the loading intensity is increased from 2000 to 3000 kg/m, 
the corresponding range of the optimum weights is 13.45- 
32.32 kg/m and 16.87-38.77 kg/m respectively. In each case 
the constraint on the normal stress in the beam is seen to 


be active . 

However, in the case of simple support conditions, 
for a loading intensity of 1000 kg/m, a change in the span 
fran 3.0 to 6.0 m enhances the optimum weight frctn 12.84 to 
29.70 kg/m. In each case, the constraint on the deflection 
in the beam is found to be active, except for a span of 
5.0 m wherein constraints on the normal stress rnd the 
deflection are both found to be active, as indie t 
Table 6.2. For the loading intensity of 2000 kg/m, the 
optimum weight ranges from 16.87 to 40.01 kg/m an 
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failure of the wall by sliding on the foundation. When the 
height is 6.0 m , the active constraint is on the minimum stem 
root thickness to resist bending moment. 

When the top width has values of 0.25 and 0.30 m, 
the respective range of the optimum cost is Rs. 577.87-2171.79 
and Rs. 630.70-2197.80, as the wall height varies from 3.0 to 

6.0 m. As before, the active constraint is the lower bound 
on the stem root thickness for smaller wall heights and that 
on the safety against sliding, for larger heights. 

6. 2. 4. 4 Silt 

For a top width of 0.20 m and the range of heights 
from 3.0 to 6.0 m, the optimum cost varies from Rs. 496.81 to 
Rs. 2157,44. For heights upto 4.0 m, the active constraint is 
the lower bound on the stem root thickness; between 4.0 to 

5.0 m, it is on the safety against sliding while for larger 
heights, it is the minimum stem root thickness to resist 
bending moment. 

For top widths of 0.25 and 0.30 m, the optimum cost 
variation is from Rs. 523.92=2185.26 and Rs. 584.50=2193.39 
respectively, as the wall height changes from 3.0 to 6.0 m. 
The active constraint is the lower bound on the stem root 
thickness for smaller heights and that on the safety against 
sliding, for larger heights of the wall. 
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Table 6.7 

Optimal sectional Dimensions 


Loading; uniformly distributed 

Intensity; 3000 kg/m ; Lateral Eccentricity; nil 



cantilever 

. 

Q JvJctll 

in in 


Optimal dimensions 

m mm 


1 

1 

Optimum J 

Active 


bfi 

j t£ l 

i bf 2 ! 

1 » 

tf 2 ; 

i 

hw ; 

i 

tw 

weight J 

in kg/m ; 

1 

1 

constra- 

int/s 

1.0 

73.4 

6.7 

78.4 

6.7 

142.0 

4.7 

13.45 

* 

g 13 

1.5 

100.0 

7.3 

100.0 

7.8 

209.0 

5.5 

21.29 

g 13 

2.0 

130.4 

8.7 

130.4 

8.7 

259.0 

6.2 

30.52 

** 

g 13' g 15 


*Maximum normal stress m the beam. 
**Maximum deflection in the beam. 
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Table 6.17 

Optimal Sectional Dimensions 

Loading; linearly varying 

Intensity; 3000-1000 kg/m ; Lateral Eccentricity; nil 

| Cantilever 

Active 
constra- 
in t/ s 


1.0 

63.0 

6 . 5 

63.0 

6.5 

113.0 

4.2 

10.11 

none 

1.5 

86 . 3 

6.9 

86.8 

6.9 

167 .0 

5.0 

15.92 

** 

g is 

2.0 

100.0 

8.2 

100.0 

8.2 

218. C 

5.7 

22.66 

g 15 



**Maximum deflection in the beam 
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Table 6.27 


Optimal Sectional Dimensions 


Loading; uniformly distributed 


Intensity; lOOO kg/m ; Lateral Eccentricity; 0.2 bf. 


(from Table 6.3) 


Span 
in m 


Cantilever 


Optimal dimensions m mm 



Optimum 
weight 
in kg/m 


Active 
constra- 
int /s 


1.0 58.0 6.4 58.0 6.4 108.0 4.1 9.36 


1.5 86.4 6.9 86.4 6.9 166.0 5.0 15.82 g. 


2.0 100.0 8.1 100.0 8.1 216.0 5.6 22.34 


*tiaximum normal stress m the beam, 
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